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SECTION  I 


ADAPTION  OF  THE  SAP  IV1  COMPUTER  CODE 
TO  THE  AIRCRAFT  SHELTER  ANALYSIS  PROGRAM 


INTRODUCTION 

1.  GENERAL  BACKGROUND 

The  United  States  Air  Force  must  continually  develop  and  improve 
shelters  and  shelter  systems  for  the  protection  of  aircraft.  These 
protective  facilities  must  function  satisfactorily  under  a number  of 
external  environments,  the  worst  of  which  is  weapon-induced  structural 
loading.  Efficient  design  and  evaluation  of  protective  shelters  for  a 
variety  of  static  and  dynamic  loads  necessitates  use  of  advanced  com- 
puter techniques.  Such  computer  analyses  must  be  capable  of  determining 
the  linear  response  of  structures  with  geometrical  configurations  typical 
of  aircraft  shelter  components. 

2.  OBJECTIVE  AND  SCOPE  OF  STUDY 

This  study  was  to  select  an  existing  computer  code  based  upon  the 
finite  element  method  of  structural  analysis  and  to  revise  this  code  to 
enable  the  Air  Force  to  effectively  analyze  shelter  development.  The 
ability  of  the  code  to  predict  the  linear  response  of  shell  and  folded 
plate  structures  subjected  to  arbitrary  static  and  dynamic  loading 
conditions  was  of  primary  importance. 

The  scope  of  this  effort  as  outlined  in  the  technical  directive  was 
limited  to  the  modification  of  an  existing  computer  program  to  obtain  a 
code  with  the  following  characteristics: 

Elements  to  be  included  in  the  finite  element  approach: 

a.  A plate  finite  element. 

b.  A shallow  or  semi-thick  shell  finite  element. 

c.  Beam-columns  and  plate  or  shell  stiffeners. 


1.  Bathe,  K.  J.;  Wilson,  E.  D. ; Peterson,  F.  E. , SAP  IV  - A STRUCTURAL 
ANALYSIS  PROGRAM  FOR  STATIC  AND  DYNAMIC  RESPONSE  OF  LINEAR  SYSTEMS, 
Earthquake  Engineering  Research  Center,  Report  No.  EERC-73- 11,  June 
1973,  Revised  April  1974,  College  of  Engineering,  Berkeley,  California. 
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Materials  to  be  modeled  (treated  as  linearly  elastic): 

a.  Reinforced  concrete  (cracked  or  uncracked). 

b.  Plain  concrete. 

c.  Structural  steel  (including  cold  rolled  corrugated  steel). 

The  nodal  boundary  restraint  options: 

a.  Axial. 

b.  Lateral. 

c.  Moment  resisting. 

d.  Variable  torsional  stiffness  (i.e.,  end  wall  simulation). 

Typical  input  parameters  (free  format  preferred): 

a.  Geometrical  shape. 

b.  Material  properties. 

c.  Static  loading  conditions. 

d.  Dynamic  load  versus  histories. 

e.  Output  control  indicators. 

i 

Output  data: 

a.  Displacements. 

b.  Reactions  at  boundaries. 

c.  Stress/strain. 

d.  Principal  stress/strain. 

e.  Displacement/stress  versus  time  histories. 

A plot  capability  of  input  and  output  data  was  required. 

Although  some  portions  of  the  code  were  rewritten,  the  purpose  was 
not  to  develop  a completely  new  algorithm  or  program,  or  to  test  physical 
models.  The  effort  was  limited  to  analytical  modeling. 
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The  approach  to  the  modification  of  the  code  was  as  follows: 

a.  A hierarchy  of  problems  to  be  sequentially  addressed  and  solved 
by  the  computer  code  while  the  final  code  was  being  developed  was  formu- 
lated and  submitted  for  approval. 

b.  Concurrently,  existing  codes  were  reviewed  and  a code  was  selected 
based  on  technical  requirements  of  this  project. 

c.  The  selected  code  was  modified  as  necessary  and  the  selected 
problems  were  solved  by  using  the  revised  code. 

This  report  describes  the  formulation  of  the  program  and  validation 
of  the  analytical  technique. 
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SECTION  II 

SELECTION  AND  USE  OF  A COMPUTER  CODE 


1.  SELECTION  OF  A COMPUTER  CODE 

Several  computer  codes  incorporating  many  of  the  features  required 
for  this  project  were  considered.  Codes  such  as  ANSYS  and  MARC  were 
quickly  eliminated  because  they  contained  the  ability  to  handle  nonlinear 
material  and  geometric  relations  which  was  not  required.  In  addition, 
these  codes  have  proprietary  restrictions  which  render  them  unavailable 
for  scrutiny  and  modifications.  The  proprietary  restraint  also  applies 
to  STARDYNE,  which  otherwise  would  have  been  an  excellent  possibility. 
Other  codes  were  rejected  because  of  insufficient  documentation  or  lack 
of  verification. 

2 

Two  codes,  NASTRAN  and  SAP  IV,  were  reviewed  in  some  detail.  It 
was  recommended  that  SAP  IV  be  adopted  as  the  basic  source  code  for  the 
following  reasons: 

a.  The  code  contained  a large  number  of  features  which  fulfilled 
technical  requirements  of  the  project. 

b.  The  code  had  been  used  for  a sufficient  period  of  time  to  detect 
and  correct  the  major  errors.  On  the  other  hand,  the  code  does  contain 
some  of  the  latest  advances  in  finite  element  theory  and  time  integration 
schemes. 


c.  The  code  is  documented,  open,  and  available  for  modifications. 

d.  NASTRAN  also  contains  many  of  the  same  capabilities  as  SAP  IV, 
but  it  permits  nonlinear  analyses.  This  latter  feature  greatly  adds  to 
the  complexity  of  the  code,  making  it  more  difficult  to  modify.  Also, 
the  generalized  nature  of  NASTRAN  results  in  the  need  for  more  computer 
time  and  storage  than  that  required  for  the  same  problem  run  with  the 
use  of  SAP  IV. 

2.  THE  USE  OF  SAP  IV 

A version  of  SAP  IV  prepared  by  Wright  Patterson  AFB  was  obtained 
and  installed,  and  initial  problems  were  run  with  very  little  diffi- 
culty. Other  than  the  usual  problems  associated  with  interpreting 


2.  NASA  SP-223,THE  NASTRAN  PROGRAMMER'S  MANUAL,  Douglas,  F. , Editor, 
National  Aeronautics  and  Space  Administration,  Washington,  D.  C.  1970. 
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statements  in  a user's  manual,  no  difficulties  were  encountered.  The 
code  performed  efficiently  and  the  results  compared  quite  favorably  with 
classical  results  for  certain  specific  test  problems.  Furthermore,  on 
the  basis  of  numerical  experimentation,  it  was  found  that  the  modal 
solution  technique  for  dynamic  problems  appeared  to  be  much  more  effi- 
cient than  the  direct  integration  method. 

A set  of  illustrative  problems  was  chosen  to  exercise  the  various 
features  of  SAP  IV  that  were  of  primary  interest  in  this  study.  In 
particular,  the  static  and  dynamic  integration  routines  were  to  be  used 
on  problems  utilizing  beam,  plate,  and  thick  shell  elements  only.  The 
complexity  of  these  problems  ranged  from  simple  beams  and  plates  to 
structures  such  as  aircraft  shelters  and  doors.  The  description  of  the 
problems,  together  with  selected  results,  is  given  in  a later  section. 

When  SAP  IV  was  chosen,  it  was  clear  that  certain  modifications  and 
additions  would  be  necessary  to  meet  the  technical  requirements  of  this 
project.  Additional  items  became  apparent  with  the  effort  associated 
with  obtaining  solutions  to  the  set  of  illustrative  problems.  The  tasks 
that  had  to  be  addressed  consisted  of  the  following: 

a.  Since  SAP  IV  is  based  entirely  on  linear  theory,  a systematic 
approach  for  analyzing  reinforced  or  plain  concrete  members  was  not  con- 
tained in  the  SAP  IV  User's  Manual.  A procedure  for  computing  equivalent 
properties  for  concrete  beams,  plates,  or  shells  was  required. 

b.  The  input  format  for  SAP  TV  is  very  rigid  with  such  anomalies  as 
blank  cards  from  time  to  time  to  denote  zero  values,  the  end  of  a section 
of  input  data,  and  problem  termination.  For  the  beginning  user,  such 
rules  can  be  very  time  consuming.  Consequently,  the  need  for  a free 
format  input  program  was  quite  evident. 

c.  SAP  IV  did  not  contain  a subroutine  for  computing  principal 
stresses  and  strains.  This  subroutine  is  a desirable  feature  for  deter- 
mining the  possibility  of  failure. 

d.  SAP  IV  did  not  contain  a mesh  plot  routine,  making  the  verifica- 
tion of  input  geometrical  data  extremely  difficult.  Frequently,  for 
very  large  problems,  the  only  method  is  to  use  engineering  judgment  on 
the  nature  of  the  results  to  a specific  loading  function.  For  dynamic 
problems  this  is  a very  risky  approach.  Thus,  the  addition  of  the 
capability  for  plotting  mesh  points  was  assigned  a very  high  priority. 

e.  The  output  of  SAP  IV  consists  entirely  of  numerical  values  and 
plots  constructed  on  the  printer.  For  dynamic  problems,  such  data  can 
be  overwhelming  and  difficult  to  interpret.  The  need  for  a time  history 
plotting  package  then  became  evident. 
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f.  The  use  of  the  SAP  IV  modal  solution  technique  for  dynamic 
problems  frequently  resulted  in  the  diagnosis  that  the  storage  limit 
was  exceeded.  Apparent  reason  for  this  is  that  the  generalized  force 
vector  is  computed  for  all  time  steps  before  the  modal  equations  are 
integrated.  Consequently,  any  limitation  on  computer  memory  allocation 
also  limits  the  number  of  time  steps,  or  the  duration  of  time  for  which 
a solution  can  be  obtained. 

g.  It  was  discovered  that  a limited  number  of  generalized  displace- 
ments could  be  requested  in  output.  This  appears  to  be  a flaw  within 
SAP  IV  and  should  not  be  a major  one  since  there  is  no  similar  limita- 
tion on  generalized  stresses. 

Some  of  these  limitations  were  expected;  others  were  not  and  signi- 
ficant effort  was  required  to  verify  that  these  problems  were  associated 
with  the  code  and  not  the  input  data.  Because  of  the  limited  time 
available  to  complete  the  project,  the  highest  priority  was  placed  on 
items  a,  b,  d,  and  e with  the  other  items  to  be  considered  only  if  time 
became  available.  The  results  of  this  part  of  the  project  are  presented 
in  the  next  section. 


SECTION  III 


ADDITIONS  TO  THE  PROGRAM 


1.  EQUIVALENT  STIFFNESS  PARAMETERS  FOR  CONCRETE  SECTIONS 

Many  aircraft  shelters  and  other  structures  of  interest  to  the  Air 
Force  are  made  of  concrete  or  reinforced  concrete.  Since  such  structural 
members  are  designed  to  crack  under  highly  probable  loading  conditions, 
an  analysis  must  reflect  this  nonlinear  behavior  in  some  fashion.  If  a 
code  (such  as  SAP  IV)  based  on  linear  theory  is  to  be  used,  properties 
for  these  structural  members  must  be  chosen  so  that  the  actual  nonlinear 
behavior  is  approximated  or  averaged  in  some  sense. 

If  a procedure  can  be  established  for  beams,  then  the  extension  to 
plates  and  thick  shells  is  relatively  straightforward.  The  procedure 
used  for  obtaining  an  equivalent  bending  stiffness  of  a beam  can  be 
applied  directly  to  orthotropic  plates  by  using  the  beam  formulas  in 
each  of  the  principal  material  directions.  There  is  a factor  l-v2(v  is 
Poisson's  ratio)  that  is  not  present  in  beam  theory.  However,  in  light 
of  the  assumptions  made  in  arriving  at  an  equivalent  linear  stiffness, 
the  incorporation  of  such  a term  is  not  difficult. 

The  use  of  the  thick  or  thin  shell  element  in  SAP  IV  introduces 
additional  complications  in  that  actual  elastic  constants  rather  than 
bending  stiffnesses  must  be  satisfied.  Because  bending  is  the  dominant 
mode  of  deformation  for  most  engineering  structures,  it  is  imperative 
that  the  resultant  bending  stiffness  correspond  closely  to  that  of  the 
actual  structural  element.  A convenient  way  of  achieving  this  corres- 
pondence is  to  choose  the  appropriate  elastic  constants  for  concrete 
and  then  select  values  for  the  thickness  and  mass  density  so  that  the 
resultant  bending  stiffness  and  mass  per  unit  shell  reference  area  are 
equal  to  the  experimentally  or  analytically  determined  values  for  that 
particular  section.  If  a different  bending  stiffness  exists  in  the 
direction  normal  to  that  already  modeled,  the  selection  of  appropriate 
orthotropic  elastic  material  properties  in  the  element  can  be  used  to 
approximate  the  different  stiffness. 

The  basis  for  equivalent  plate  and  thick  shell  element  properties 
is  contained  in  the  assumptions  used  for  beams.  Thus,  the  remainder  of 
this  subsection  is  devoted  to  beams  with  several  illustrative  examples. 
Extensions  to  plates  and  thick  shells  are  considered  to  be  rather  straight- 
forward and  are  only  included  in  the  sample  problems. 

When  reinforced  concrete  members  deform  under  the  action  of  some  load, 
cracking  of  the  concrete  occurs  on  the  tensile  side  of  the  member.  This 
highly  nonlinear  behavior  must  be  approximated  with  an  appropriate  linear 
property  if  a linear  elastic  analysis  is  to  be  used  with  any  degree  of 
accuracy. 
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Unsymmetrical  sections  introduce  another  complicating  factor  in 
that  the  cracked  stiffness  for  both  directions  of  bending  may  be  dif- 
ferent. As  a reasonable  approximation  for  problems  in  which  any  one 
section  may  experience  both  positive  and  negative  bending  moments,  it  is 
proposed  that  an  additional  averaging  procedure  be  used.  First,  obtain 
the  average  stiffnesses  of  the  cracked  and  uncracked  sections  for  bend- 
ing in  both  directions.  Then  use  the  average  of  these  two  stiffnesses 
to  obtain  a resultant  linear  equivalent  stiffness  to  be  used  in  a 
dynamic  analysis. 

3 

The  approximate  formula  suggested  by  Biggs  for  the  average  bending 
stiffness  of  rectangular  sections  is 

E bd3 

EcIb  = (5.5ps  +0.083) 


in  which 


Ec  = Young's  modulus  for  concrete, 

IB  = average  second  area  moment, 
b = width  of  the  member, 
d = effective  depth  of  the  member,  and 


p$  = steel  ratio. 


Since  it  is  not  clear  just  what  assumptions  were  used  in  deriving 
this  equation,  it  is  proposed  that  first  principles  be  used  in  obtaining 
average  stiffness  parameters.  The  following  examples  illustrate  the 
procedure,  and  Biggs'  formula  is  used  for  comparison. 

Example  1.  Symmetric  Section. 

Consider  a rectangular  beam  (Figure  1)  with  the  same  reinforcing 
top  and  bottom  so  that  the  average  stiffness  is  the  same  for  both  posi- 
tive and  negative  bending. 


T.  Biggs,  J.  M. , Introduction  to  Structural  Dynamics,  McGraw  Hill  Book 
Co.,  New  York,  N.  Y.,  1964. 


As  =0.75  in2, 


Figure  1.  Symmetric  Rectangular  Section 

Area  of  steel,  top  and  bottom,  A^  * A^  = 0.75  in2  with  distances 
dp  = dn  = 12  in.  The  steel  ratio  is,  by  definition, 

As+As 

ps  = (8  i n ) ( 1 2 in)  = °'0156 
Young's  modulus  for  concrete  is  given  by 

Ec  = (w)1’5  33/FJ" 

in  which  w is  the  weight  density  in  lb/ft3  and  f'  is  the  concrete  strength 
in  psi.  For  4000  psi  concrete 

E = (145)1,5(33)^000  = 3.64  x 106  psi 


The  ratio  of  steel  to  concrete  Young's  moduli  is 


n = Es/E 


c 


29  x 10 6 psi  _ 
3.64  x 10b  psi 


Thus  the  steel  area  is  transformed  to  an  effective  concrete  area  8 times 
as  large.  For  the  uncracked  section  (Figure  2)  the  moment  of  inertia 
becomes 


Figure  2.  Transformed  Uncracked  Section 
A = TT  bh3+2nAs(5  in^  = TI  (8)(14)P+2(6  in^(25  in2)  = 2129  in4 


For  the  cracked  section  (Figure  3),  the  location  of_the  neutral 
axis  must  be  determined  first.  Let  the  neutral  axis  be  y from  the  upper 
surface. 


Figure  3.  Transformed  Cracked  Section 


The  total  equivalent  concrete  area  is 


Then  first  moments  about  the  top  surface  yield 


Now,  the  second  area  moment  about  the  neutral  axis  for  the  cracked 
section  is 


m 

1 //  E 

m 

Ijj  A = (1/3)8  in  (3.32  in)3+6  in2 (3.32  in  -2  in) 

+6  in2 (12  in-3.32)2  = 560  in4 

The  average  second  area  moment  of  the  cracked  and  the  uncracked  section  is 
^ = 2129  jn^.560  in4  = 1345  in4 


The  value  obtained  from  Biggs  is 


!B  = T”  (5-5ps  +0.083)  = 


For  the  purpose  of  comparison 


_ (8  in) (12  in) 


{5.5(0. 0156)+0. 083}  = 1168  in4 


= 1.15 


Example  2.  Unsymmetric  Reinforcing 

Consider  a rectangular  beam  with  different  reinforcing  at  the  top 
and  bottom. 


In  Figure  1,  let 


A*  =0.5  in2 


A°  = 1.25  in2 


d*  = db  = 12  in 


The  steel  ratio  for  the  cross  section  is 


ps  * (Qiii?yi25)  ■ °-°182 


I 


The  transformed  uncracked  section  is  illustrated  in  Figure  4. 


nA$  = 4.00  in2 


7.  in-y  1 <>> 


.£>  •>» 


‘nA,.  = 10.00  in2 
s 


Figure  4.  Transformed  Uncracked  Section 


Total  equivalent  concrete  area 


A = (8)(14)+4. 0+10.0  = 126  in2 


The  distance  from  the  bottom  face  to  the  neutral  axis  is  y where 


yA  = (112)(7)+(4.0)(12)+(10)(2)  = 852  in3 
y = 6.76  in 


The  second  area  moment  about  the  neutral  axis  is 

jJa  = y^-{8) (14) 3+(8) (14) (7-6. 76)2+4.00(12-6.76)2+10(4.76)2=  2172  in' 


and 


The  average  value  for  the  second  area  moment  for  positive  bending  becomes 
,p  . 1 (2172+825).  U99 

cl  ^ 

The  transformed  cracked  section  for  negative  bending  is  shown  in 
Figure  6. 


Figure  6.  Transformed  Cracked  Section  for  Negative  Bending 

The  total  equivalent  concrete  area  is 

A = 8y+(4)+10  = (8y+14)  in2 

With  reference  to  the  bottom  face,  the  location  of  the  neutral  axis  is 
given  by 

yA  = | (8y)+(2) (10)+(12) (4)  = 4y2+68 
or 

4y2+14y-68  = 0 
y = 2.73  in 
and 

!N  A = I (8)  (2.  73)3  +10(2.  73-2)2+4(12-2. 73)2=  403  in'* 


ft 
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The  average  value  for  the  second  area  moment  for  negative  bending  is 


ij  = 1(403+21 72)  = 1288  in“ 


The  average  value  for  both  positive  and  negative  bending  is  then 

'a  ' ’ 1394 


Biggs'  formula  yields 

!B  = 5. 5P$+0 .083  = (M1 2I3(5 . 5 (0 . 01 82 ) +0 . 083 1 = 1266  in'* 


and 


-p-  = 1.10 

*B 


Example  3.  Unsymmetrical  Cross-section  with  Reinforcing 

The  following  cross-section  (Figure  7)  has  been  suggested  as  a 
suitable  model  for  a section  of  corrugated  material  from  an  aircraft 
shelter. 


f'  = 5000 

Ec  = 4 . 07  x 106  psi 
n = Es/Ec  = 29/4.07  = 7.13 

The  transformed  uncracked  section  is  shown  in  Figure  8.  The  equivalent 
concrete  area  is 


Figure  8.  Transformed  Uncracked  Section 

A = (6) (8)+(5) (4. 5)+(2) (1 .07)+0. 36  = 73.0  in2 

and  taking  moments  about  the  top  face,  the  location  of  the  neutral  axis 
is  given  by 

yA  = (0. 36) (1 .25) +(48) (3)+(l . 07) (5. 75)+(22. 5) (3.25)+( 1 .07) (9.625)  = 346.5  in3 

y = 4.75  in 


With 

it  follows  that 


The  second  area  moment  is 

!N  A = T^(6)3+48(4.75-3)2-Hfl(5)(4.  5)3+(22.5)(2.25+1.25)2 
+(0. 36)  (3. 5)2  +(1 .07)  (5. 75-4. 75)2  +( 1 .07)(  5. 75-0.875)2 
= 636  in4 
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Figure  9.  Transformed  Cracked  Section  for  Positive  Bending 
A = 8y+0. 36+1. 07+1. 07  = (8y+2. 5G)  in2 

and  taking  moments  about  the  top  face  the  neutral  axis  is  obtained  from 

yA  » (8#|)+(0.36)(1.25)+(1,07)(5.75)+(1.07)(9.625)  = 4y2+16.9 
or 

4y2+2. 50y- 16.9  = 0 


and 


The  second  area  moment  is 


The  average  value  for  positive  bending  is 


^(636+98)  = 367  in 


The  transformed  cracked  section  for  negative  bending  is  shown  in 
Figure  10. 


Figure  10.  Transformed  Cracked  Section  for  Negative  Bending 


The  equivalent  concrete  area  is 


and  the  location  of  the  neutral  axis  is  obtained  by  taking  moments  about 
the  bottom  face 

yA  = 5y(^)+(l .07) (0.875 )+(l .07) (4. 75) +(0.36) (9.25)  = 2.5y2+9.35 


2. 5y2+2. 50y-9.35  * 0 


y = 1.50  in 

The  second  area  moment  is 

!N  A = (1 . 50)3+(l .07) (1 .50-0. 875 )2+(l .07 )( 4.75-1 .50 )2 

+(0. 36) (9.25-1 . 50)2=  39  in1* 

The  average  value  for  negative  bending  is 

I1?  = %(636+39)  = 337  in4 

a 

For  both  positive  and  negative  bending,  the  recommended  value  is 

Ta  = Js(Ia+Ia)  = 352  in<* 

Biggs  does  not  suggest  a procedure  for  a section  as  complicated  as 
this  one.  However,  his  approach  could  be  adapted  as  follows:  The  steel 
cross  section  area  is  A,  = 0.05+0. 1 5<+0. 15  = 0.35  in2.  For  positive 
bending  use  b=8  in  and  0=9.625  in.  Then 


rB  = {5<5ps  +0-083}  ■ (8) -*H5. 5(fe|5'^+0. 083)  = 385  in4 


For  negative  bending  use  b = 5 in  and  d = 9.25  in 


l|j  - (5)-1^25l3{5.5  (0^35^+0.o83}  = 247  1n* 
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The  average  of  the  two  values  Is 


IB  = ^(385+247)  = 316  in' 


T-  - 1.11 


For  the  three  beam  sections  presented,  the  beam  stiffnesses  calculated 
by  averaging  four  stiffnesses  (cracked  and  uncracked  for  positive  bending 
and  cracked  and  uncracked  for  negative  bending)  are  from  10  percent  to 
15  percent  greater  than  those  determined  from  the  Biggs  formula.  The 
Biggs  formula  is  not  directly  applicable  to  non-rectangular  sections,  and 
for  complicated  sections  its  use  becomes  highly  questionable.  For  simple 
cross-section,  however,  its  use  in  determining  dynamic  bending  stiffnesses 
can  be  justified  from  the  foregoing  test  calculations.  For  complicated 
beam  cross-sections,  although  parameters  may  be  chosen  so  that  the  Biggs 
formula  stiffness  compares  well  with  the  four  stiffness  average  method, 
it  is  suggested  that  only  the  latter  method  be  used  in  determining  the 
dynamic  stiffness. 

2.  FREE  FORMAT  INPUT  PROGRAM 

To  alleviate  the  time  consuming  and  confusing  format  rules  for  intro- 
ducing data  into  the  SAP  IV  code,  a program  known  as  FFIP  was  developed. 
FFIP,  short  for  Free-Format  Input  Program,  enables  the  user  to  prepare 
input  data  for  SAP  IV  with  considerable  ease  and  speed  and  to  avoid  the 
rigidity  in  form  necessary  for  creating  input  data  directly  applicable 
to  the  SAP  IV  program.  Ultimately,  input  data,  when  properly  prepared, 
can  be  completely  shuffled,  and  the  resultant  output  of  FFIP  will  be 
identical  every  time.  FFIP  is  a separate,  independent  program  to  be  run 
just  before  executing  the  SAP  IV  program  in  a single  job  run,  and  it 
creates  a resultant  output  tile.  Tape  1,  which  is  to  be  used  as  input 
to  the  SAP  IV  program.  To  utilize  FFIP  proficiently,  the  user  is  required 
only  to  be  familiar  with  FFIP  input  syntax  rules,  which  are  simple  and 
straightforward,  and  FFIP  keyword  parameters.  Explanations  of  FFIP 
input  requirements  and  further  elaborations  of  FFIP  characteristics  are 
available  in  Appendix  A. 

3.  GRAPHICS  PLOT  PACKAGES 

Two  plotting  packages  were  developed  for  the  SAP  IV  code  to  simplify 
both  data  input  and  output.  For  dynamic  problems  the  time  history 
package  provides  the  user  with  plots  of  displacement  or  stress  versus 
time.  The  mesh  plot  package  permits  rapid  verification  of  input  geometrical 
data. 


3.1  Time  History  Plot  Package 

A time  history  plot  package  is  available  for  the  user  who  wishes 
to  analyze  voluminous  SAP  IV  output  with  ease  and  efficiency.  The  plot 
package  has  been  included  in  SAP  IV  without  major  modification  to  the 
code  and  the  standard  software  package  is  required  for  Its  execution. 
Plots  can  be  created  by  requesting  printer  plot  from  SAP  IV  ( i . e . , KKK-2) 
with  data  for  specific  component  members.  No  other  input  data  are  re- 
quired and  each  resultant  plot  contains  a maximum  of  three  curves  which 
are  automatically  grouped  depending  on  component  members  specified.  A 
normalization  factor  is  provided  for  minima  and  maxima  of  all  three 
curves  when  applicable.  Also,  the  beginning  and  ending  time  data  points 
are  automatically  plotted  to  the  extrema.  Each  curve  permits  up  to  100 
points.  If  there  are  more  than  100  points,  the  data  will  be  incremented. 
Further  elaborations,  including  samples  and  functional  descriptions  of 
the  program  are  given  in  Appendix  B. 

3.2  Mesh  Plot  Package 

A mesh  plot  package  is  available  for  rapid  and  simple  verifi- 
cation of  input  geometric  data.  It  is  a separate  program  which  can  only 
be  run  after  a successful  execution  of  SAP  IV  in  a single  job  run.  In 
addition  to  the  normal  input  data  to  the  SAP  IV  program,  a variable 
named  MODEX  in  the  master  control  card  section  of  the  SAP  IV  manual 
must  be  set  to  1.  SAP  IV,  subsequently,  produces  an  output  file.  Tape  8, 
which  provides  one  of  the  required  Inputs  to  the  mesh  program.  Besides 
TAPE  8,  card  input  must  include  data,  such  as  visual  orientations  of 
generated  object  defined  in  3-D  space  and  selections  of  optional  node 
and  element  labelings  and  various  or  all  elements  to  be  plotted.  Card 
input  is  again  in  free-format  and  is  easy  to  use.  Explanations  of 
syntax  rules,  a list  of  keyword  parameters  and  further  details  of  mesh 
plot  characteristics  are  available  in  Appendix  C. 


PROBLEMS  THAT  DEMONSTRATE  THE  APPLICABILITY  OF  SAP  IV 
TO  AIRCRAFT  SHELTER  ANALYSIS 


1 . BACKGROUND 

A hierarchy  of  problems  was  selected  to  illustrate  the  use  of  SAP  IV 
for  obtaining  solutions  in  cases  that  ranged  from  the  most  elementary 
type  to  those  that  were  fairly  complex  and  representative  of  the  type 
that  might  be  encountered  in  a detailed  analysis  of  an  aircraft  shelter. 
The  simpler  problems  were  chosen  to  demonstrate  basic  features  of  the 
code  and  to  allow  comparison  with  solutions  obtained  from  strength  of 
materials  formulas.  The  more  complex  problems  were  designed  to  represent 
the  structures  and  load  environments  actually  encountered  in  Air  Force 
design  and  experiments.  Since  the  objective  of  this  part  of  the  study 
was  to  determine  the  appropriateness  of  SAP  IV  for  this  class  of  problems, 
no  attempt  was  made  to  perform  the  kind  of  detailed  analysis  that  might 
be  required  under  a given  situation. 

Of  the  large  number  of  elements  available  in  SAP  IV,  the  following 
problems  were  solved  using  only  beam,  thin  shell  or  plate,  and  thick 
shell  or  plate  elements.  Input  for  each  problem  was  introduced  into 
the  code  using  the  Free  Format  Input  Program  (see  Appendix  A).  The  mesh 
plots  and  time  history  plots  have  been  selected  so  that  over  the  range 
of  problems  considered,  most  of  the  options  available  with  the  plot 
package  are  exercised. 

2.  CANTILEVER  BEAM 

2.1  Static  Analysis  Under  Five  Different  Loads  (Figure  11) 


2.1.1  Problem  Description 

A 25-foot  long  cantilevered  W 12  x 40  beam  was  subjected 
to  the  following  concentrated  loads  acting  on  the  free  end  of  the  beam: 

a.  1000- lb  axial  load 

b.  1000-lb  vertical  load 

c.  300, 000- in- lb  moment 

d.  5000- in- lb  torque 

e.  All  of  the  above 


i 


I 
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THEORY  Ax  = 9.4044  x 10'4in  THEORY  Ay  = -1.0011  in 

SAP  IV  A = 9.4044  x 10"4in  SAP  IV  A = _i.0088  in 

NODE  6 NODE  6 y 
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2.1.2  Five  3-D  beam  elements  were  used  in  the  analysis.  The 
assignment  of  nodes  and  beam  elements  is  shown  in  Figure  12. 


Figure  12.  Cantilever  Beam  Element  Assignment 


2.1.3  Input  Data  - Five  Load  Cases 

HEADER=*STATIC  ANAL.  CONC.  LOADS  CANT.  BEAM 
NUMNP=7 ,NELTYP=1 ,LL=5 
IX(1,7)=1 ,1,1,1 ,1,1 
I X ( 2 6)=0,0,1 ,0,1,0, 

XYZT ( 1 )=0. ,0.  ,0.  ,0. , 

XYZT(2)=60. ,0.  ,0.  ,0.  , 

XYZT(6)=300. ,0.  ,0.  ,0. , 

XYZT(7)=0. ,60. , 

NBEAM=5,BNEPC=1 ,BNMPC=1 , 

BMPC(1 )=2.9E7,0.3,7.339E-4,.289 

BEPC ( 1 ) = 1 1 . ,3. 51  ,4. 128,. 956,44. 1 ,310.  , 

BEAM(1)=1 ,2,7,1 ,1, 

BEAM(5)=5,6,7,1 ,1 , 

CLMD(6,1 )=1000. , 

CLMD(6,2)=0. ,-1000.  , 

CLMD(6,3)= , , , , ,3. 0E5, 

CLMD(6,4)=, , ,5. 0E3 

CLMD(6,5)=1000. ,-1000. ,,5.0E3,,3.0E5 

I 


2.1.4  Results 


I 


A comparison  of  deflections  as  computed  by  SAP  IV  and  a standard 
strength  of  material  approach  was  made  for  each  of  the  five  load  cases. 
The  maximum  difference  in  computed  deflections  was  1.56  percent,  and 
this  can  probably  be  attributed  to  the  fact  that  deflections  computed 
using  strength  of  materials  did  not  take  shear  deformations  into  account, 
whereas,  SAP  IV  did.  Because  the  problem  analyzed  was  statically  deter- 
minate, no  appreciable  differences  in  shears  and  moments  were  computed 
by  either  approach. 

Load  Case  1 (1000  lb  axial)  X-translation 

SAP  IV  - node  6 Ay  = 9.4044  x 10"H  In 

Strength  of  Materials  A = Pi/AE  = (1000  lb)^  Q 106psi-j 

= 9.4044  x 10-l*in 


Load  Case  2 (1000  lb  vertical)  y- translation 

SAP  IV  - node  6 A^  = -1.0088  in 

Strength  of  Materials  4y  = ^ = (1000  l^lzglffipsilSlO  in 

= 1.0011  in 

Difference  in  predicted  deflections  is  0.8  percent. 

Load  Case  3 (300,000  in-lb  concentrated  moment)  y-translation 
SAP  IV  Node  6 Ay  = +1.5017  in 


Load  Case  4 (5000  in-lb  torque)  Angle  of  Twist 

SAP  IV  x rot.  = 0.14067  radians 
Node  6 

Strpnath  of  Materia*  n - 11 (5QQQ  Jn-lb)  (3QQ  in) 

Strength  of  Materials  p JG  (0.956  in1* ) (11 .154  x 103  in4 ) 

= 0.14067  radians 

Load  Case  5 (Combined  loads)  y-translation 
SAP  IV  Ay  = +0.49289  in 
Node  6 

Strength  of  Materials  Ay  = -1.0011+1.5017  = 0.5006  in 
Difference  in  predicted  deflections  is  1.56  percent. 

2.2  Effect  of  Shear  Deformations  on  a Deep  Beam  and  a Standard  Beam 
2.2.1  Problem  Description 

As  seen  in  problem  2.1,  the  inclusion  of  shear  deformations 
in  an  analysis  leads  to  computed  deflections  somewhat  larger  than  de- 
flections computed  by  standard  means.  In  an  attempt  to  study  the  effect 
shear  deformations  have  on  the  deflections  computed  for  a standard  beam 
and  a "deep"  beam,  analyses  were  made  on  a standard  15-foot  cantilever 
and  a deep  15-foot  cantilever.  For  each  beam,  analyses  were  made  in- 
cluding shear  deformations  and  not  including  shear  deformations.  A 
"deep"  beam  was  considered  to  be  a beam  with  a length  to  depth  ratio  of 
5:1  or  less  (see  Figure  13). 

The  following  problems  were  considered: 


1 


CASE 

1 


CASE 

2 


10,111  lb 


50,000  lb 


Figure  13.  Shallow  and  Deep  Cantilever  Beams 


In  both  cases  loads  were  chosen  that  would  yield  a bending  stress  of 
approximately  20.5  ksi  ac  the  fixed  end. 

2.2.2  Three  3-D  beam  elements  were  used  for  both  problems.  The 
arrangements  of  the  nodes  and  beam  elements  used  in  the  analysis  are 
shown  in  Figure  14. 


Node 


Element 


Figure  14.  Cantilever  Beam  Element  Assignment 


2.2.3  Input  Data  - Case  I with  Shear 


HF  A PFR  = *STO • BEAM  SHEAR  W12X65 
NUMNP  = 5 »NFLTYP= 1 »Ll  = l , 

IX(1,5)=1, 1,1, 1,1,1, 

T X ( 2 4 ) =0,0 , 1 , 1 , 1 ,0  , 

XYZT  ( 1 1 =n.  ,0.  ,0.  ,P. , 

xyzt ( ? ) =6o. ,n. ,n. , 

y YZ  T ( 4 ) =] qo,  , 

XYZT ( 5 ) =0.  , 5.  , 

N|PFAM=3  ,BMFPC=1  »BNMPC=1» 

PMPC ( 1 ) =?.QF7,  . 3, IS.? 2 ,490. , 

RFPC<l)=19.1,4.73,4.73,2.19,17S.,b3- 

PFAM(1)=1,2,5,1,1 

RRAM ( 3 ) =3 ,4,5 , 1 ,1 

FL^PI 4, 1>=,— 10111., 


2.2.4  Input  Data  - Case  I without  Shear 


HFAOFR=*STn.  NO  FHFAR  W12X65 
N<IMNP  = S ,Mci.TYP  = 1 ,LL  = 1 , 
IX(1,5)=1, 1,1,1 ,1,1, 

T X ( 2 4)=n,n,i ,1,1 ,n, 

XYZT ( 1 ) = R » ,0, » 0 • » 0 , , 

VY7T ( 2 ) =fin. ,o. ,0. » 

X V7J ( 4 ) =1 RO. , 

XYZT ( 5 ) =n. ,S. , 

MPFAM  = q ,RMFPr = 1 ,PNMPC=1 , 

R MP  c ( 1)=2.PE7,.3,1C.22,4  90., 
RRPO ( 1 ) = 10. 1 , , ,2. 19, 1 7 b, ,b33. , 

RRAM(  1 )=  l ,2,4,1  ,1 

orA^(3)=3,4,S»l ,1 

0L^n(4, 1 j = ,-101 11 . , 


Case  II  with  Shear 


HFAPFR=*nFFP  RF AM  FHFAR  '.-.'36X1  7F 
NIIMMP=F  ,NFLTYP  = 1 ,LL  = 1 
I X ( 1 * 5 ) = 1 , 1,1,  1*1,1, 

ix<  ? 4)=n«o*i»i»i.o. 

XY2T ( 1 ) =n.  , 

XYZT(2)=6n.  , 

XYZT(4)=1R0.. 

XV7T ( 5 ) = n.  ,5.  , 

MRFAM=3  *BMFPC=1  *RNk/PC=1  » 

RMPC( 1)=?. of 7*. 3, IF. 22  » 4 P 0 • , 

RFPC ( 1 )=39«8*21»26*21«26. 7*03*226 
RFAM ( 1 )=1,2*5»1»1 
RFAV(3)=3,4,b,l*l 
r|  vn  ( 4 , 1 ) = ,-sno-OO.  * 


Case  II  without  Shear 


HF A nER  = *PFFP  REAM  NO  SHEAR  W 3 6 X 1 3 b 

NIImNP*5»NELTYP=1  *LL  = 1 
!X(l*b)  = l.l. 1*1*1. 1* 

T X ( 2 4) =0.0.1. 1*1*0. 

XY7T ( 1 ) =0. , 
xv7T(21=Ao. , 

XVZT (4 ) =1«0.  * 

<v7T(F  )=o.  , F . , 

wpr  am  = ^ ,R.vrpr=  1 ,PNf,PC  = 1 * 

nwpfl 1)=2.PF7*.3»1F.22*4Q0., 

OFPC(l)=39.8,*.7.n3.226.,782n. 

RFAM ( 1 ) = 1 »?*5*1*1 
PFAM( 3 ) =3.4 ,5 . 1 .1 
r|  Mp<4.]  ) = ,-FOonn.. 


I 


2.2.7  Notes  on  Input  Data 

The  input  data  for  the  two  cases  involving  a standard 
beam  (W  12  x 65)  are  very  similar  to  the  input  data  discussed  when 
considering  the  analysis  of  a cantilever  beam  under  five  loading  cases. 
One  should  note  that  for  the  analysis  considering  shear  deformations, 
the  element  geometric  properties  include  an  area  of  4.73  square  inches 
that  is  associated  with  shear  forces  in  the  local  two  direction,  whereas, 
for  the  second  case  the  input  does  not  include  this  piece  of  data  and 
the  analysis  will  not  include  shear  deformations. 

For  the  two  cases  involving  a deep  beam  (W  36  x 135) 
it  should  be  noted  that  for  the  first  case,  where  shear  deformations  are 
taken  into  account,  the  element  geometric  properties  include  an  area  of 
21.26  square  inches  associated  with  shear  in  the  local  two  direction. 

In  the  second  case,  where  shear  deformations  are  not  taken  into  account, 
this  shear  area  has  not  been  included. 


2.2.8  Results 


For  the  standard  beam  (W  12  x 65)  the  maximum 
deflection  computed  by  SAP  IV  including  shear  deformations 
was  A = 1.3061  inches.  The  maximum  deflection  computed  by  SAP  IV  not 
Including  shear  deformations  was  A = 1.2716  inches  which  compares  to 
within  five  places  with  the  strength  of  materials  solutions.  The  error 
in  th^^trength  of  materials  solution  was  2.7  percent. 

For  the  deep  beam  (W  36  x 135)  the  maximum  deflection 
computed  by  SAP\£V  including  shear  deformations  was  A = 0.46656  inch. 

The  maximum  defleHjon  computed  by  SAP  IV,  not  considering  shear  defor- 
mations, was  a = 0?>2561  inch,  which  again  compared  to  within  five 
places  with  the  strength  of  materials  solution.  The  error  in  the  strength 
of  materials  solution,  this  time,  was  8.9  percent.  Thus,  for  extremely 
deep,  short  beams,' the  error  encountered  when  computing  deflections 
without  considering  shear  deformations  begins  to  be  significant. 

Therefore,  In  any  analysis  of  an  aircraft  shelter  where  there  is  some 
question  as  to  whether  or  not  the  shelter  section  may  constitute  a 
"deep"  beam,  shear  deformation  should  be  taken  into  account.  As  the 
Input  shows,  this  can  be  done  with  very  little  effort.  All  that  is 
required  Is  to  include  in  the  element  geometry  input  the  areas  associ- 
ated with  shear  forces  in  the  local  two  direction.  (Figure  15.) 
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1 


Case  1 


Case  2 


SAPIV  WITHOUT  SHEAR  DEFORMATIONS  Ay  = 1.2716  in 
SAPIV  WITH  SHEAR  DEFORMATIONS  Ay  = 1.3061  in 


SAPIV  WITHOUT  SHEAR  DEFORMATIONS  Ay  = 0.42861  in 
SAPIV  WITH  SHEAR  DEFORMATIONS  Ay  = 0.46656  in 


ELEMENT  MODEL 


3.  Fixed- Fixed  Beam 


3.1  Static  Analysis  of  Fixed-Fixed  Beam  (Figure  16) 

3.1.1  A box  beam  fixed  at  both  ends  is  subjected  to 

i 

1.  A concentrated  load  at  the  center 

2.  A uniformly  distributed  load 


Load 
Case  (1) 

t= 

Material  = Steel 

Member:  4 in  x 4 in  x 40  in  box  beam  with 
wall  thickness  1/4  inch 

Figure  16.  Fixed-Fixed  Beam  Load  Cases 


100  lb 


Load 

Case  (2) 


25  lb/in 

4T.I  F ****  + **  WiHt 

a ^ 


3.1.2  A group  of  four  3-D  beam  elements  is  used  for  this 
analysis  with  the  location  of  the  nodes  and  the  element  numbers  (square 
boxes)  shown  in  Figure  17. 

+y 

j 1 2,3  4k  Element 


* Node 

1 2 3 4 5 0 

in  ^10  in  ' 10  iJlO  in^ 

Material  = Steel 

Member:  4 in  x 4 in  x 40  in  box  beam 
with  wall  thickness  1/4  in 


4 


17  riomant  flcci'nnmonf  _ Fi  var(_Fi  Pnamc 


■•P-.J-WWWW'  7 H ..  ..  _ _ 


The  input  data  for  these  problems  are  as  follows: 
3.1.3  Input  Data  - Case  1 


HFADER=*STATIC  BEAM  PROR  LEY  * 
NUMNP=ft »NFLTYP=1 »LL=1 * 

TX( 1 »5»6)=1*1*1»1*1»1» 

T x ( ? 4)=o»o»l »i »i *n» 

XYZT  ( 1 ) = 

XYZT  ( 2 ) =1 n.  , 

XYZT ( 5 ) = 4n. * 

XYZT (6 ) =0.  , 1 . , 

MRF AM  = 4 . ^mfpc=1»PNMPC=1» 

PMPf(  1 )=3#0F7».333».OOO7-a46 
PFPC ( 1 ) =1 . 59»  » ». 079 *.766*. 76ft 
RFAM( 1 ) =1 ,2*6*  1 » 1 
RFAM(4)=4,5»6»1 »1 
CLMD  ( 3» 1 ) =,-100. ♦ 


3.1.4  Input  Data  - Case  2 


HEADF R = * STATIC  BEAM  PROBLEM 
MUWNP=6 *NrLTYP=l »LL  = 1 
T X ( I * ft  * ft  ) = I * 1 * 1 * 1 * I * 1 » 

T X ( 2 4 ) =n,Ot  1 » 1 » 1 *0 
XYZT ( 1 ) = 

XYZT ( 2 I =1 o.  » 

XYZT ( 5 ) =4n. » 

XYZT ( 6 ) =0 • »1.  ♦ 

N R F A M = 4 »BMEPC= 1 »BNMPC=  1 

rmpo( 1 )=3.0F7».333». 0007346 *-1^.72 

RFPri  1 )=i  ,6Q»»».079*.7fift*.76fi 

R F A M ( I ) = 1 »?  »6*1 *1  * 

RFAM(4)=4*ft»6»l »1 * 

FLM(1  1 = 1. 

OFrL  FY=  I . * 
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3.1.5  Comparison  of  SAP  IV  Results  with  Strenqth  of 
Materials  Solution 

Case  1.  Concentrated  Load  at  Center. 


100  lb 


Maximum 
deflection 
(at  center) 

V 

Shear  force  50 
diagram  (lb) 


SAP  IV  Results 
- 0.14505  x 10‘*  in 


Strength  of  Materials  Results 
- 0.14505  x 10"2  in 


20  in 


20  in 


20  in 


20  in 


(lb-in) 

Bending  moment 


500 

yiTk 


DUU 


Manual  of  Steel  Construction,"  Seventh  Edition,  AISC,  pp  2-203. 

Figure  18.  Static  Results  - Fixed-Fixed  Beam  Point  Load 

The  results  of  Fiqure  18  Indicate  the  following: 

1.  The  maximum  deflections  agree  to  within  five  significant  figures 

2*  The  shear  force  diagrams  are  the  sc<me,  as  are  the  bending 
moment  diagrams. 


Case  2.  Uniformly  Distributed  Load. 


i 


25  lb/in 


Maximum 
deflection 
(at  center) 


Shear  force 
diagram 


SAP  IV  Results  Strength  of  Materials  Results 

-0.72527xl0"2in  -0. 72527xl0"2in 


10  in  10  in 

10  in  10  in 

20  in 

20  in 

Figure  18a.  Static  Results  - Fixed-Fixed  Beams  - Distributed  Load 


The  results  of  Figure  18a  show  that: 


1.  The  maximum  defections  agree  to  within  five  significant  figures. 


2.  The  difference  between  the  two  shearing  force  diagrams  and  the 
two  bending  moment  diagrams  is  due  to  the  assumptions  of  constant  element 
strain  in  SAP  IV. 
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3.2  Dynamic  Analysis  of  Fixed-Fixed  Beam 

3.2.1  A box  beam  (Figure  19)  fixed  at  both  ends  is  subjected  to 
a time  dependent  load  at  the  center  as  follows: 


F(t) 


40  in 


Material  = Steel 

Member  = 4 x 4 x 40  inches  box  beam 

with  wall  thickness  = 1/4  inch 


Figure  19.  Fixed-Fixed  Beam  - Dynamic  Point  Load 

3.2.2  The  finite  element  model  and  location  of  nodes  are 
the  same  as  for  the  static  analysis  (3.1.2). 

Two  different  analysis  methods,  mode  superposition 
and  direct  step-by-step  integration,  were  used.  The  input  data  for 
these  analyses  are  as  follows: 
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3.2.3  Mode  Superposition 

HFAnFR  = *nvNA^IC  ANALYSIS  OF  A FIXED-FIXED  BE -AM 
NL'MNP  = 6.NFLTYP=1  »NF  = ft  ♦ NT>YN  = 2 . 

I X ( 1 * 5 , 6 ) = 1 . 1 ♦ 1 ♦ 1 . 1 ♦ .1  . 

r y ( ? ,i*o. 

*Y7T  ( 1 ) - 

XY7T ( i ) =1  A.  , 
yY7T(S|=4o,, 

W7T(6)=A*.»I.. 

moc  am  = /,  ,r»Mr  nr  - i ,pM‘.'Pfs  1 

nf‘nr(  ] | ri  ,nf7,  , 'ni  , .noATi/.f, 

orpf ( 1 ) =1 . SO,  , , .079,  . 76b  . .766 

nCAW( 1 1=1 .7,^,1 *1 

BFAM(/»)a/4»S,6tl*l 

1FPP-1  «NFN  = ] ,KT  = 4nnn,NUT  = lon.0T:,no'ii 

MPn«2)*l**l. 

Ml  p ( i » =■?  ,_ion. 

t ( i »i ) =n.  ,t  ( l »? ) = .os,t  ( l .3 1 = s.n 

F(1  »1  ) =0 • « F ( 1»2)=1«»F(  1 » 3 ) = 1 • 

y y Y = 7 

frovP(7,7,*i ,7,f 

v y v c,  - 2 

TM  2*2  l =1  »?  ♦*»  7»8»1  2 
T M 2 » 3 1 *1  . 2 »<S 

• 

3.2.4  Direct  Integration 

H FAnFR=*DYNAwIC  ANALYSIS  OF  A FIXtD-FIXKD  BEAM 
m|  jM»'P  = ft  ,NFL  TYP=1  » NF  = f>  *ND YN  = A 
T Y ( 1 « 5 » A ) = 1 ♦ 1 *1  » 1 * 1 » 1 * 
yx(7  <o  =o»0 » 1 ♦ 1 ♦ 1 *o » 

W7T  ( i ) = 

yv7T(21=1^. * 

XV7TI * ) =4p. * 
yv7T|h)56o,,1  , , 

»inrA*':4  inMrDOil  »PN“Dr= 1 » 

n"Pf ( 1 1 *1 .07, , 333 » ,0no734A» 

pcPr(l)=1.AO,,..n7P*.7AA..7AA 

npAMf  1 ) =1  ,?  *6» 1 » ] 

n F A M { A ) = A*S»f>»l  » 1 

mfm=  i *m  i = ',nno , NOT  = 1 no  * r> T = «o inoi 

N P ( 3 * 2 ) s 1 * * 1 • 

NLP ( I ) =3 ,_i no. 

T(1,1)=0.»T(1.?)=.0S»T(1.*)=S. 
f(l ,1 ) =0. »F( 1*21=1. *F(1»3)=1. 

Y Y Y = 7 » 

I ro»'P  (25=1  ,7,  A 
T rOMP (31  = 1 * 7 » A 
vyv  c = p , 

T S ( ? ♦ 2 > » 1 . ? * fi . 7 * 8 « 1 ? 

TM2*3)*l.2*A 


^ .... 


3.2.5  Results 


Both  the  vertical  and  horizontal  displacements  of  the 
beam  center  (node  3)  are  presented  in  Figure  20.  The  horizontal  com- 
ponent is  zero  for  all  time  while  the  vertical  component  increases 
negatively  to  an  average  value  corresponding  to  the  static  deflection. 


DYNAMIC  ANtAl  Y5 1 5 Zc  A c IXED-FIXEi?  3FAM 
NORMALIZED  DISPLACEMENT "RESPONSE  AT  NODE 


NuRKRt  l?E0  FRCTGR-. 


162F-G2 


COMPONENT  l 


COMPONENT  2 


Figure' 20.’  Midpoint  Displacement  Response  (Fixed-Fixed  Beam) 


4.  PLANE  FRAMES 

4.1  Static  Analysis  of  Rectangular  Frame  Under  Four  Different 
Loads. 

4.1.1  Problem  Description 

A 10-foot  by  10-foot  rectangular  frame  was  subjected 
to  four  different  loading  conditions.  These  four  load  cases  are  illus- 
trated In  Figure  21.  F 


10  Ib/ln  Uniform  Load 


1000  lb  Cone.  Lateral  1000  lb  Cone.  Vertical 
Load  Load 


Combination  of  all  Loads 


Figure  Zl.  Rectangular  Frame-Static  Load  Cases 

The  frame  was  composed  of  W 8 x 15  horizontal  and  vertical 
structural  steel  members  with  the  following  properties: 

A = 4.4321n  ; 1^  = 48.1  in'1,  I = 3.40  In4;  J = 0.136  In4 


4.1.3  Input  Data  - Rectangular  - Frame  - Static 


HFADFR=*STAT  IC  ANAL.  RECT.  FHRAN'E  UNIF.  AND  CONC.  LOADS 
M(iMNP=1  3,NFLTYP=1  ,LL  = 4 
T X ( 1 » 1 3 ) = 1 »1»1»1»1»1» 

XY7T ( 1 ) = 

XYZT ( 2 ) =0. .30. » 

XYZT  ( 5 ) =0. » 120. 

XYZT ( 6 ) =3o. , 120.  » 

XYZT  < 9 ) =1?0. *120.  ♦ 

XY7T ( 1 0 1 =1 20.  ,9o.  , 

XYZT ( 13 ) =1 20. >0. , 

IX ( 2 1?  1 =O,0» 1 * 1 »1  .0, 

N R F A Y = 1 2 ,on|FPC=1  ,RMMPC  = 1 , 

RMPC( 1 )=2.9F7,.3,7.34E-4,.284 
R*rPC(l)=4«43»»».136»3»4»48.1  * 

RRAM( 1 1=1. 2,  6, 1,1 
opaw ( 4 ) =4, 5 ,6 , 1 , 1 
RFAM(5l=5,6»4»l»l 
nTAM( 8 ) =8,9,4, 1 » ] 

RFAM(9)=9, 10,8, 1,1 
RFAM( 12 1 = 1 2 » 13,4, 1 ,1 
CLMP( 1,11=150. 
clmp( 1,41=150. 

01 MP( 2, 1 1=300. 

FLMD( 2,4  1=  300  . 

CLMD  ( 3, 1 1=300. 

CLMD ( 3,4) =300. 

C L M O ( 4,1) = 300. 

'LMP( 4 ,4 ) =300. 

CLMD ( 5 * 1 1 = 150. *-150. » 

CLMD ( 5 ,4 1 =1 50.  ,-150.  , 

CLMD( 6, 1 1 = ,-300. , 

CLMD( 6 ,4 1 = , -300.  , 

CLMOI 7,1)  = ,-300.  , 

CLMD1 7,2  1=1000.  , 

C L M D ( 7,3 ) = ,-1000.  , 

CLMD ( 7,4 1 =1 000. ,-1300.  , 
ri  vp(  R,,  1 ) = ,-300,  , 
ry  M n ( R.4)  = .-700.  , 

CLMD ( 9 ,1 ) =-150. ,-150.  , 

OLM 0(9,4)=-!  50.  ,-150.  , 
rLMD ( 1 o , 1 < =-300.  , 

CLWD ( lo ,4 1 =-300.  , 

CLMD (11  ,1  ) =-300.  , 
n MD(  11  ,4 1 =-300.  , 

CLMD ( 12,1) =-300.  , 

C L M P ( 12,4) =-300.  , 

CLMD  ( 1 7 ,1  i =-]  50.  , 

O.MPt  1 7 ,4  1 =-1  50.  , 


— — 

! 

4.1.4  Results.  A comparison  of  moments  in  in/lb  as  computed  by 
SAP  IV  and  a standard  moment  distribution  technique  for  the  first  three 
load  cases  is  shown  in  Figure  23. 


Load  Case  1 


W = 10  lb  /in 


Moment  Distribution  SAP  IV 


Maximum  Difference  = 13.4  percent 


Load  Case  2 


Moment  Distribution  SAP  IV 


Maximum  Difference  =0.7  percent 


Moment  Distribution  SAP  IV 

Maximum  Difference  =0.4  percent 

Figure  23.  Static  Results  - Rectangular  Frame 

It  should  be  noted  that  the  maximum  difference  in  computed  moments 
occurred  for  load  case  1.  This  is  probably  because  the  uniform  load  was 
approximated  by  equivalent  concentrated  loads  applied  at  the  various  nodes. 
A more  accurate  analysis  could  be  achieved  by  including  more  beam  elements 
in  the  numerical  analysis. 
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4.2  Response  Time  History  of  a Rectangular  Frame  (Figure  24) 
Subjected  to  a Traveling  Pressure  Wave 


4.2.1  Problem  Description 


The  same  10-by  10-foot  rectangular  frame  considered  in 
problem  4.1  was  struck  by  a pressure  wave  of  10  lb/in  traveling  at 
200  ft/sec. 


W * 10  lb/in 
t^  = time  of  arrival 


Figure  24.  Dynamic  Loading  - Rectangular  Frame 

4.2.2  Finite  Element  Model 

The  same  12  3-D  beam  elements  as  were  used  in  problem 
4.1  were  used  in  this  problem.  The  assignment  of  nodes  and  beam  elements 
is  shown  in  Figure  25. 


Figure  25.  Element  Assignment  - Rectangular  Frame 
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4.2.3  Input  Data  - Rectangular  Frame  - Travelling  Pressure  Wave 

HFAnFR=* <FCT.  FRAME  TRAVELING  PRESS.  WAVE 
NUMNP=13»NFLTYP=1 ,NF=lo*NDYN=2 , 

T X ( 1 » 13 ) =1 , 1 , 1 » 1 » 1 . 1 . 

T X ( ? 121=0,0,1 ,1,1 ,n, 

XY7T  < 1 1 = 

XYZT I 2 1 = o. ,70.  , 

XY7T ( 5 ) =o.  , 1 20.  ♦ 

XY7T < 6)=3o. , 1 ?o.  , 

XY7T ( 9 1=170.,  1 ?o.  , 

XV7T( 10 1 = 1 20.  ,90.  » 

XVZT ( 1? 1 = 120.  ,30*  » 

XY7T( 131=170., 0., 

N op AM =12 ,qNFPC=l »BNMPC=1 , 

pmpc(  1)=?.9F7».3,7.34E-4,.28A 
<3FPC(  1 1=4. 43,  ,,.136, 3. 4,48.1, 

°FAM( 1 1=1  .2 ,6 , 1 , 1 
qfaV(4i=4,5,6*1 ,1 
QFA w( S ) =5 , A, 4 , 1 ,1 
RPAV( 8)=8,9,4,1  ,1 
PFAM (91=9, 10,8,1  ,1 
r»FA  A'  (171  = 17, 13, 8, 1,1 

TFPR=1 ,NFN=1 ,NAT  = 3,NT  = 50nO,0T  = .0002  »NOT  = 10 
•MP(  ’ *1  1 =1  ,-l  , lr’0. 

WD(  7,11  = 1,1  ,300. 

N P ( o,i  i=l,l  ,300. 

NP( 4, 1 ) =1  ,1 ,300. 

NP( ^.11=1  ,1  , 1 So . 

NP  ( 5,21=1,1,-150. 

HP ( 6 ,2 ) =1  ,7 ,-7po. 

NP( 7,2 1 =1 ,7,-300. 

NP( R,21=l  ,7,-300. 

A'  D ( O , 1 )=1  .7,-1  40. 

MP( o,7)=l ,7, -ISO. 

MD ( 10,1 )=1 ,3,-300. 

NP( 11*1 1=1 ,3,-300. 

NP( 1 2* 1 1 =1 ,3,-300. 

NP(  13,1  1 = 1 ,3,-1  SO. 

AT  ( 1 1=0., AT  (2)  = .025,AT  (3  ) = • 0 i>  » 

nlp ( i ) = ; ,1 . 

T ( 1 , 1 ) =0. ,T (1,21=2. 

0(1,1  1 = 1.  ,F(  1,71  = 1, 

X K K = 2 » T GO'iP  ( 1 » 7 , S 9,11,131  = 1,2,3,6 
y vy  o = p , 

T S ( ?,1  1=1  ,7,7,6,7,8,9,17 
IS (7,21=1, 7,°, 6, 7, 8, 9,1? 
tS( 7,41=1 .7,3,6,7,8,9,12 
TS(p,Sl=  1,7,7,6,7,8,9,17 
IS (?, 61=  1, 7,7,6, 7, 8, 9,1? 

T S ( ?,7)=  1,7,3,6,7,8,9,17 
IS (2, 81=  1,7,7,6,7,8,9,12 
IS  (7, 9)=i, 7, 7, 6, 7,8,9,]7 
1 S(?,lll=1,7, 7, 6, 7, 8, 9,1? 

TS(7,1?)=1, ?, 3, 6, 7,8,9,17 


4.2.4  Results 


I 

1 


f 


The  first  part  of  the  results  gives  the  eigenvalues  and 
mode  shapes  of  the  frame.  Figure  26  shows  the  first  five  eigenvalues 
and  the  mode  shapes  associated  with  the  first  three  eigenvalues. 
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W]  = 144.9  Rad/sec 


Wz  = 566.1  Rad/sec 


W4  = 979.2  Rad/sec 
W5  = 1817  Rad/sec 


Figure  26.  Rectangular  Frame  - Natural  Frequencies  and  Mode  Shapes 

The  displacement  time  history  for  node  7 is  given  in 
Figure  27.  Components  1 and  2 in  this  figure  correspond  to  the  x and  y 
displacements,  respectively.  To  obtain  absolute  results  the  normalized 
values  given  in  each  figure  must  be  multiplied  by  the  normalized 
values  given  in  each  curve.  For  the  stress  time  history  shown  in  Figure 
28,  component  1 is  the  axial  force  at  node  7,  component  2 is  the  shear 
force  at  node  7 and  component  3 is  the  out  of  plane  shear  force  at  node 
7,  zero  for  all  time.  All  of  these  components  act  on  element  7. 

The  displacement  response  of  node  7 exhibits  a period  of 
0.043  second,  which  corresponds  to  the  first  mode  vibration  period. 
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4.3  Static  Analysis  of  a Polygonal  Frame  Under  a Uniform  Pressure 


4.3.1  Problem  Description 

A five-member  polygonal  frame  was  subjected  to  a 
uniform  load  (Figure  29)  of  10  lb/in.  The  corners  of  the  frame  were 
on  a 24- foot  arc  to  simulate  the  geometry  of  a typical  aircraft  shelter. 
The  members  of  the  frame  were  W 12  x 36  with  properties: 


Figure  29.  Polygonal  Frame  - Static  Loading 


A = 10,6  in2,  Ixx  = 281 .0  in*.  I = 25.5  in4 
J = 0.830  in\  E = 29  x 106lb/in2,  v = 0.3. 


Five  3-D  beam  elements  were  used  in  the  analysis. 
The  assignment  of  nodes  and  beam  elements  is  shown  in  Figure  31. 


Figure  31.  Element  Assignment  - Polygonal  Frame  (Five  Elements) 
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is  shown  in  Figure  32. 


Figure  32.  Local  Coordinate  System  of  Beam  Element 


The  node  k is  any  nodal  point  which  lies  in  the  local  1-2  plane 
but  not  on  the  1 axis  and  is  used  to  orient  the  element.  For  this 
problem  node  7 was  used  as  the  node  K for  all  of  the  beam  elements. 

It  should  be  noted  that  for  elements  4 and  5 the  node  K lies  below 
the  1 axis  and  the  signs  of  the  output  forces  with  regard  to  the  local 
coordinate  system  are  reversed. 

4. 3. 2. 3 Input  Data  - Polygonal  Frame  - Static  - Fixed 

End  Forces 


HFADF  K = *POl.YGONAL  FRAhE  STATIC  UNIF.  PRESS.  FIXED  tl\iU  FuKCEo 
NUMNP  = 7 »NFLTYP=1  ,LL  = 1 ♦ 

TX ( 1 *6,7) =1  , 1 , 1 , 1 , 1 , 1 * 

T X ( 2 8 ) = n,0  * 1 *1  » 1 »o, 

XY7T ( 1 ) = 

XY7T ( 2 ) = 68.  , 169. » 

XY7T ( 3 ) =1P0. *274.  , 

YY7T (4 ' =377. *274.  , 

X V7 T { 5 i = 6?1 . ,160., 

XY7T(6»=576.,0. 

X Y7  T (7 ) =o.  , 300. * 

MPFAM=S ,RNFDC=1 »BNFEFS=? ,RNMPC=1 * 

ROPr  ( 1 l = 2.or7*.  3.7.339F-4.  .7  HO 
R F P C ( 1 )=lo.6»*».83,26.8,281., 
nrcri  | i i=,oo4,8»*** 2 7288.8 
PFPF J( 1 ) = .On4.8  ***,-2  7289.8 
PFPFT ( 2 ) = *-006.8  » » » ,-2  7288.8 
PFFFJ( 2 ) = ,-904.8* ,,  *27288.8 
RFAM( ) )=1  ,2*7,1  *1  ,1 
R F A M ( 2 ) =2 ,7  *7, 1 , ] * 1 
OFAP  ( 3)=3, 6,7,1 *1 *1 
orA v ( 4 ! =4, 8, 7, 1 *1*2 
oFM'l  8 ) =6,6,7, 1 ,1 ,2 
Fl_M(  1 >*1  . , 
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4. 3.2. 4 Results 

The  moments  in  in/lb  acting  on  the  ends  of  the 
various  elements  as  computed  by  SAP  IV  are  shown  in  Figure  33. 


26,890 


26,890 


26,900 


26,900 


28,860 


\ 

1,860 


Figure  33.  Bending  Moments  for  Polygonal  Frame 
(Equivalent  Fixed  End  Forces) 


It  should  be  noted  that  these  values  are  all  within  5 percent  of  what 
would  be  computed  by  a moment  distribution  approach. 

4.3.3  Analysis  by  Equivalent  Concentrated  Forces  at  Six  Nodes 

4. 3. 3.1  Problem  Description 

The  five  member  polygonal  frame  (Figure  34) 
described  in  the  introduction  was  subjected  to  the  uniform  load  of  10 
lb/in.  For  this  analysis  the  uniform  load  was  replaced  by  equivalent 
concentrated  forces  at  the  six  nodes. 

4. 3. 3. 2 Finite  Element  Model 

The  uniform  load  was  represented  by  the  equivalent 
concentrated  forces  acting  on  the  six  nodes.  The  assignment  of  nodes 
and  elements  for  this  problem  Is  the  same  as  that  described  In  problem 
4.3.2. 
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Figure  34.  Element  Assignment  - Polygonal  Frame 


The  concentrated  loads  refer  to  global  coordinates  x and  y,  therefore 
the  si gn  reversal  that  was  irieritioned  in  prublem  4.3.2  does  not  apply. 


4. 3. 3. 3 Input  Data 
trated  Forces  at  Six  Nodes. 


Polygonal  Frame  - Static  - Concen 


HF ADFR=*POL YGONA 1 FRAME-  STATIC  UNIF.  PRESS  • FORCE  6 NGDc. 
MUMMP=7  »NFLTYP= I »LL  = 1 » 

T X ( 1 » 6 * 7 ) = 1 .1,1  *1  ,1  *1* 

IX( ? 5)*n*o*l»l*l *Ci 
•X  Y 7 T ( 1 ) = 

YY7T ( 2 ) = SR  . , 1 69.  * 
yY7T(3)=109.*?7A., 

YY7T ( A 1 =377. ,27A.  , 

YY7T  ( 5 ) =S?1  . ,1  69.  , 

YY7  T (6 ) =576. «n, 

XY7  T ( 7 ) =a.  , 3nn  . » 

MRFAM  = S*RN|F  PC  =1  *BNMPC=  1 * 

RMPCI 1 >=2.9F7*. 3, 7.339F-A ♦ .289 
np  or  I 1 ! r 1 n , A , . . , H 3 • ? *>  .5*281.  1 
R C fl  M ( 1 ) a 1 , ? , 7 , 1 ♦ 1 , 1 

r»r*A««  / O t _ n n ~r  i i ^ 

>•■%£>—  ( * • * • 9 i * i * 1 

o r a r ( 3 l = 3 . A » 7 * 1 * 1 * 1 
nFAMI A ) =4, S * 7, 1 ,1,1 
nP»IM(S):^,P,,7,l  .1  ,1 
CLMFX  1 , 1 ) = 861  .*-780., 

CLMr»(2«l)  = n93.,-lnl?.» 

CLMHI  3*1  )=*)32.  »- 16  37.  ♦ 

CLMD(A,l)=-532.»-1637.* 

cLMPl5*l)a-1393.»-10l2.« 

CLMH( 6* 1 1 a-861 . *-280. * 


4. 3. 3. 4 Results 


The  moments  in  in-lb  acting  of  the  end 
five  beam  elements  as  computed  by  SAP  IV  are  shown  in  Figure  35 
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J 


Figure  35.  Bending  Moments  for  Polygonal  Frame 
(Equivalent  Concentrated  Forces) 

When  compared  with  the  results  obtained  in  problem  4.3.2,  these  moments 
are  grossly  different  from  those  obtained  by  either  a moment  distribution 
approach  or  an  equivalent  fixed  end  force  finite  element  analysis. 

The  reason  that  this  equivalent  nodal  force 
approach  was  attempted  Is  that  in  a dynamic  analysis  loads  must  be  ex- 
pressed in  terms  of  equivalent  nodal  forces.  As  expressed  in  problem 
4.3.4,  a much  more  accurate  analysis  is  achieved  when  the  finite  element 
model  includes  more  nodal  points. 


Nodes. 


4.3.4  Analysis  by  Equivalent  Concentrated  Forces  at  Sixteen 


4. 3. 4.1  Problem  Description 


For  this  analysis  the  uniform  load  of  10  lb/in, 
acting  on  the  polygonal  frame  described  in  the  introduction,  was  replaced 
by  equivalent  concentrated  loads  acting  at  16  nodes. 
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4. 3. 4. 2 Finite  Element  Model 

Fifteen  3-D  beam  elements  were  used  in  the  analysis. 
The  assignment  of  nodes  and  beam  elements  is  shown  in  Figure  36,  along  with 
global  coordinates  of  the  various  nodes. 


As  stated,  the  uniform  load  was  replaced  by  con- 
centrated loads  at  16  nodes.  Figure  37  shows  the  x and  y components  of 
the  loads  that  were  used  in  the  finite  element  model. 


Figure  37.  Polygonal  Frame  - Equivalent  Concentrated  Forces  (16  Nodes) 
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HFADFR  = *POLY(~tONAL  frame  STATIC  UN  I F < 
N'|MNP=1  7»MFLTYP=1  , L L = 1 , 

TX<  I ,16»17)=1»1,1»1  ,1,1 
TX ( 2 15 )=n,0» 1 , 1 ,1 ,0, 

XY7T  ( 1 ) = 

XYZT(2)=18.»56.» 

XYZT ( 3 ) =37. , 113.  , 

XY7T(4)=55.,169., 

XYZT ( 5 ) =ln3. .204.  » 

XY7T(6)=]5i . ,239. * - 

YY7T ( 7) =199. ,274.  , 

XYZT (8) =2 5 8. ,274.,  . 

XY7T (9 ) =31 8. ,274. , 

YY7TI 10 ) =377. ,274.  , 

XV7TI 11 1=425. .239.  , 

XYZ  T ( 1?  ) =4  7 7 . ,204.  , 

XYZT ( 13 1=521 . .169. » 

XYZT ( 14  ) =579.  ,113.. 

XYZT ( 151=558.  ,56.  » 

XYZT ( 16  1=576.  , 

XYZT( 1 7 ) = , 300 • » 

M«F  AM=15,PNFPC=1 »RNMPC=1  , 


LOAD  CONC.  LuAD6  Al  ib  NODtS 


RFPT ( 1 1 =1  a. 6, , 
RFAM( 1 ) =1  , 2,  ] 7 
DFAM(2)=?,3,I7 
PPA^( 3 1=3,4, 17 
RFAP (41=4,5,17 
RFAM( 51=5,6,17 
RFAM(6)=6,7,17 
P F A M ( 71=7.8,17 
PF  AM ( 8 1 =«, 9, ] 7 


PFAM (91=9,10.17, 1,1 

PFAM( ln)=io,ll,17,i,l 

P F A M (111  = 11, 12,17, 1.1 
PFAP( l?)  = i 2,13,17,1  ,1 
PFAM( 131=13,14,17,1,1 
p F A M (14)=14, 15, 17, 1,1 
PC  A w ( 1 5 1 = 1 6 , 16,17,1  ,1 
Cl  mp>(  1 ,1  1=287.  ,-93.  , 
CLMD( 2 , 1 1=573. ,-186. , 
f L M i ) i 3 , 1 1=573.  ,-186., 
CLM 0(4,11=465. ,-337., 
CLMD(*,1  1 = 354 . ,-488.  , 
r|_MO(  6,1  1=354.  ,-488.  , 

O L M o ( 7,1 1 = 1 78.  ,-546., 

CL YD (8,1) = ,-603 . , 

CL YD (9,1 1 = ,-603  . , 

FLYO ( 1 o , i 1 =-178. ,-546.  , 
CLMD( 11 ,1 1 =-354. ,-488. , 
CLMD(  12,1 )=  — 354. ,-488.  , 
CLwD( 17,1 i =-465. ,—347. , 
rLY0( 14,1 i=-573.,-186. 
F|.YD(  15,n  =-573.  86.  , 

OL YO (16,1 i=-?87.,-186.. 


.83,25.5,281  . , 
1 ,1 
1 , 1 

1,1 

1,1 
1 ,1 

1,1 

1,1 
1 , 1 


4. 3. 4. 3 Input  Data  - Polygonal  Frame  - Static  - Concen- 
trated Forces  at  16  Nodes. 


4. 3. 4. 4 Results 


The  moments  in  in- lb  acting  on  the  ends  of 
elements  1,  3,  4,  6,  7,  9,  10,  12,  13  and  15  as  computed  by  SAP  IV  are 
shown  in  Figure  38. 


Figure  38.  Bending  Moments  - Polygonal  Frame  (15  Elements) 


The  maximum  variation  from  moments  that  would 
be  computed  by  a moment  distribution  approach  for  the  uniform  load  is 
16.3  percent.  A more  accurate  analysis  would  be  achieved  if  the  uniform 
load  were  distributed  to  more  nodes.  However,  a simpler  and  more  accurate 
static  analysis  is  achieved  by  replacing  the  uniform  load  with  equivalent 
fixed  end  forces  as  exemplified  in  problem  4.3.2.  This  technique  of 
expressing  a uniform  load  in  terms  of  equivalent  fixed  end  forces  is 
recommended  in  the  SAP  IV  manual  as  a method  by  which  uniform  loads 
may  be  specified  for  a static  analysis.  However,  for  a dynamic  analysis 
(problem  4.3.3),  a uniform  load  must  be  expressed  in  terms  of  nodal  forces 
and  this  16-node  model  will  be  used  In  the  subsequent  dynamic  analysis. 

4.4  Response  Time  History  of  a Polygonal  Frame  Subjected  to  a 
Traveling  Uniform  Load. 

4.4.1  Problem  Description 

The  five-sided  polygonal  frame  discussed  in  problem  4.3 
was  subject  to  a pressure  wave  of  magnitude  10  lb/in  traveling  at  1000 
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ft/sec.  As  in  problem  4.3.4,  the  uniform  load  was  replaced  by  equivalent 
concentrated  force  acting  at  16  nodes  (Figure  39). 


t. .=  Arrival  Time  Value 


Figure  39.  Polygonal  Frame  - Dynamic  Loading 


The  indicated  loads  are  in  units  of  pounds 


4.4.2  Finite  Element  Model 


The  same  15  3-D  beam  elements  used  in  problem  4.3.4  were 
used  for  this  dynamic  case.  Also,  the  assignment  of  nodes  and  beam 
elements  for  this  problem  was  the  same  as  for  problem  4.3.4. 


4.4.3  Input  Data  - Polygonal  Frame  - Traveling  Pressure  Wave 


HFADFR=*POLY.  FRAME  TRAv.  UNIF.  LOAD  lOOOFPS  10  LRS./IN 
mi  |MMP=  l 7 » MrL  T VP=  1 » NF= 1 0 » NOYN  = 7 , 

T X ( 1, 16, 17)=1»1»1,  1,1,1, 

1 X ( 7 1 6 1 = o , 0,1*1*  1,0, 

W7  T ( 1 > = 

X Y7T ( 2 ) =1 8 . , 66 . 

YY7T ( 3 1=77. ,113., 

XY7T(4  ) = 66.  , 169.  , 

XY7T ( 6 ) =1 03. ,704.  , 

Y V7  T ( 6 \ = 1 61  . , ? 3 9 . , 

X Y7  T ( 7) =1 09.  ,774.  , 

YV7T ( 81=768.  ,774.  , 
y V7  T ( 9 ) =31 8 . ,?7A.  , 

V Y 7 T ( 101=377.  ,?7A.  , 

X Y 7 T ( 1 1 )=4?6. ,239.  , 

XY7T ( 12  ) = A 7 3 . ,204. , 

X v 7 j ( 13)=S71. ,169., 

XY7T( I4)=6  3 0.  , i i 5 • , 

W7T(  16  )=66R.  ,66.  , 

XY7T ( 16 1=676.  , 

W7T  ( 17)=o.  .300.  . 

MOFAM=]6,RNMPr=l,RNFPC=l, 

RVDr ( 1 ) =?. OF 7, . 3, 7. 339F-4, .280 

OPPCM  1=  l o . 6 , , ,.83,26.6,731.  , 

poam (11=1,2,17, 1,1 

nr a 2 1*2, 3, 17, 1,1 

p,FAN(  31=3,4,17,1,1 

PPAV (4)=4,6,17,1,1 

PFAM ( 61=5,6,17,1,1 

°FAM (61=6,7, 17,1,] 

n r A M ( 71=7,3,17,1  ,1 

oPAM( 8 1=8,0,17,1,] 

RFAM (9l=Q,10, 17, 1,1 

OPAM (101=10, 11  ,17,1  ,1 

PPAM (I] 1=11,12,17,1,1 

O F A M ( 12  1=12,13,17,1,1 

nrAM(  131  = 13,14,17,1  ,1 

n f a w ( 141  = 14,15,17,1  ,1 

n F A M ( 161  = 15,16,17,1  ,1 

MFM=1 ,N AT =6, NT =6000, NOT =1 o,nT=.0002 

M o{l, 11=1,1 ,2R7. 

ND ( 1,21=1 ♦ 1 ,-93. 

MP( 7,1  1=1  , 1 ,673. 
mo ( 7,21=1.1 ,-186. 
mo (7, 11=1.1, 573. 

MP( 3 ,2 1=1 , 1 ,-186. 

NO ( 4 , 1 )=1  .7  .466. 


Input  Data  - Polygonal  Frame  - Traveling  Pressure  Wave  (concluded) 


MP-(  ft,?  )=1  ,7,-117. 
m°(S1  '=1  .2,354. 

NP<  * ,2 ) =1 ,2 ,-48R. 

VP (6,1 )=i ,?,354. 

MP( 6 .2  ) = 1 ,2 ,-488. 

NP( 7,1) =1 ,3,178. 

NP( 7,2 ) =1 ,3,-546. 

MP( 8 ,2 ) =1 ,3,-603. 

MP( 9,2 1 =1  ,3  ,-603. 

VP(  10,7  1 = 7 ,17,-1  78. 
v D ( 7 0,21  = 7 .0  , -846 . 

NP( 11,1 )=1 ,4,-354. 

ND(  1 1 ,2 1=1  ,4,-488. 

MP ( 12,1  1 = 7 ,4,-354. 

NP ( 1 ? ,2 ) = 1 ,4,-488. 
wp ( j 3 , 1 ) =7  ,4,-465. 

NP i 13,2)-’  ,4,-337. 

\i  P ( 74,1  ) = 7 ,6,-573. 

MP(  ] 4,? ) = 1 ,6,-186. 

MD ( 16,7 1=7 ,6,-673. 

M D ( 16,21=7  ,6,-186. 

MP ( 14,1 )=i ,6,-287. 

M P ( 7 6 , 2 ) = 7 ,6,-03. 

ATI) ) =0 • , AT (2)=. 0071 ,Ar(3)=.02*AT(4)=. 0339, AT (5)=. 0434 
NLP ( 1 ) = ? ,7  . 

T ( 1 , 1 ) =0. ,T ( 1 ,2  )=2. 

F ( 1 , 1 ) = 1. ,F(  1 ,2  1 = 1 . 

K k'K  =2  , 

TOfiMP  (1  ,4,7.10,13,161  = 1,2,3,6 
'4rKS=7 

1 6 ( 2,1  )=1  .7,6,7,3,12 
74(2,41=1,2,6,7,8,12 
7 4 (7, 71  = 1,2,6,7, 8, 12 
74(2,101=7 .7,6,7,8,12 
T 4 ( 2 , 13)  = 7 ,2,6,7^8,72 
14 (?, 161=7,2,6,7, 8, 12 


4.4.4  Results 


The  first  part  of  the  results  give  the  first  10  eigenvalues 
and  eigenvectors  of  the  frame.  Figure  40  shows  the  first  three  mode  shapes 


56.5  Rad/Sec 


123.9  Rad/Sec 


/ W3  = 226.4  Rad/Sec 

Figure  40.  Polygonal  Frame  Eigenvalues  and  Eigenvectors 


Time-history  response  curves  for  node  7 are  shown  in  Figures  41 
and  42.  Components  1 and  2 refer  to  the  x and  y motion,  respectively. 
The  rotation  response  curves  refer  to  z-axis  rotation.  The  expected 
period  of  vibration  for  all  of  these  curves  is  approximately  equal  to 
the  period  of  the  first  mode  of  vibration  of  0.1112  second. 

POLY.  FRAME  1RRV.  UNIF.  LORD  1000FPS  10  LBS. /IN. 

NORMALIZED  DISPLACEMENT  RESPONSE  RT  NODE  7 


NORMRL 1 FED  FRCTOR- 


. 196E+00 


COMPONENT  1 
COMPONENT  2 


Figure  41.  Displacement  Response 
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POLY.  FRAME  TRRV.  UNIF.  LOAD  1000FP5  iO  LBS. /IN. 
NORMALIZED  ROTATION  RESPONSE  AT  NODE  7 


NORMALIZED  FACTOR - 


TIME  (SECOND) 


. 128E-02 


COMPONENT 


S 


Figure  42.  Rotation  Response 
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FORCE  RESULTANTS 


I 


^ y jsr  ^ 't»-  ry 


For  the  time  histories  shown  in  Figure  43,  component  1 is  the 
axial  force  at  node  7 and  component  2 is  the  shear  force  at  node  7. 
All  of  these  components  act  on  element  7. 

POLY.  FRAME  1RAV.  UNIF.  LOAD  1000FPS  10  LBS. /IN. 
NORMAL  1 7ED  STRESS  RESPONSE  IN  ELEMENT  7B 


NORMALIZED  FACTOR- 


TIME (SECOND) 

,nLir  n,  COMPONENT  1 

.SOUE+Qi*  COMPONENT  2 


Figure  43.  Force  Response  Element  7 


4.5  Response-Time  History  of  a Reinforced  Concrete  Polygonal  Frame 
4.5.1  Problem  Description 

A five-sided  polygonal  frame  with  geometry  similar  to 
the  frames  discussed  in  problems  4.3  and  4.4  was  subjected  to  a travel- 
ing pressure  wave.  The  frame  was  made  up  of  members  that  were  2-foot 
sections  of  a reinforced  concrete  covered  corrugated  steel  arch  air- 
craft shelter.  The  cross-section  is  shown  in  Figure  44. 


Figure  44.  Shelter  Section 


For  the  analysis  an  equivalent  reinforced  concrete  section  was  chosen 
by  matching  flexural  rigidities  of  the  uncracked  sections.  This  section 
is  shown  in  Figure  45. 


65 


Figure  45.  Equivalent  Reinforced  Concrete  Section 


Steel  In  the  lower  portion  of  the  model  with  an  area  of  5.22  In2 
was  selected  to  match  the  area  of  steel  In  the  corrugated  steel  panel 
and  was  placed  at  a position  which  coincided  with  the  center  of  gravity 
of  the  steel  panel.  It  should  also  be  noted  that  the  area  of  concrete 
in  the  lower  portion  of  the  model  was  equal  to  the  area  of  concrete  in 
the  steel  panel. 

The  approximate  moment  of  inertia  of  the  shelter  section  cast  in 
the  corrugated  steel  panel  was  determined  by  transforming  the  rein- 
forcing steel  and  panel  steel  into  an  equivalent  amount  of  concrete. 

Two  methods  were  used.  The  width  of  the  cross-section  was  treated  as  24 
l/2-1nch  wide  strips.  Using  this  technique,  the  moment  of  inertia  of  an 
uncracked  2-foot  wide  section  of  the  shelter  was  computed  to  be  56,034 
in'*.  The  moment  of  inertia  of  the  uncracked  equivalent  reinforced  section 
(Figure  45)  is  58,665  in'*.  This  is  a difference  of  only  4.5  percent. 





The  method  for  selecting  the  equivalent  stiffness  of  a concrete 
section  to  use  In  a dynamic  analysis  has  previously  been  discussed. 
Using  these  procedures  an  effective  stiffness  of  36.189  In1*  was  used, 
along  with  an  area  of  711.8  in2,  a torsional  rigidity  of  20,000  in  an 
a weak  axis  flexural  rigidity  of  20,000  in4.  It  should  be  noted  that 
the  latter  three  values  do  not  affect  the  results  of  this  analysis. 

The  pressure  wave  used  for  the  analysis  was  selected  to  model  the 
free  field  overpressure  data  traces  from  the  mixed  company  event  at  a 
range  of  600  feet.  The  wave  selected  had  a maximum  pressure  of  36.7 
psi,  a duration  of  160  milliseconds  and  a velocity  of  2000  ft/sec. 

The  wave  chosen  is  shown  in  Figure  46. 


time  (milliseconds) 


Figure  46.  Pressure  Wave 


For  the  2-foot  wide  section,  the  peak  pressure  of  36.7  psi  was 
equivalent  to  a uniform  load  of  881  1 b/ In,  which  was  replaced  by 


equivalent  concentrated  forces  acting  at  16  nodes.  These  equivalent 
forces  are  shown  in  Figure  47  and  are  also  in  pounds. 


M 50480 


50480 


25240 — + 


Figure  47.  Equivalent  Nodal  Forces  - Concrete  Polygonal  Frame 


The  arrival  times  in  seconds,  of  the  wave  traveling  at  2000  ft/sec  are 
also  Indicated  in  the  sketch. 

4.5.2  Finite  Element  Model 

The  same  16  nodal  points  used  in  problems  4.3  and  4.4 
were  used  for  this  analysis.  Fifteen  3-D  beam  elements  were  used  to 
model  the  five-sided  polygonal  frame. 


r 


i • 


4.5.3  Input  Data  - Concrete  Polygonal  Frame  - Dynamic 

HFAnFR  = *POLY.  FRAME  MIXED  COMPANY  600  FT.  2000  F PS 
NUVNP=  1 7 » NFL  T YP  = 1 * NF  = 1 0 ♦ Nr>YN  = 2 , 

IX(1*16»17)=1, 1*1*1, 1*1 
TX( 2 15 )=o,0» 1 * 1 • 1 *0 
XY7T ( 1 1 = 

XYZT ( 2 ) =18. *56. 

XYZT(3)=37.»113. 

XY7T (4 ) =56.  ,169. 

X Y7  T (6 )=]o3.  ,204. 

XY7T( 61=161 . ,239. 

XY7T ( 7) =1 09. ,774. 

XV7T ( 

XYZT ( 9 ) =31 8. ,274. 

X YZ  T ( 1 O ) =3  77, ,274. 

XY7T( 11 1=426. ,239. 

YY7T i 12  1=473,  ,204 , 

XYZT ( 13 ) =521 .,169. 

XYZT ( 14 ) =630.  ,113. 

XY7T  f 1 6 ) =668. ,66. 

VY7T  ( 1 6 1 =r-76.  , 

XY7T( 171=o. ,80o., 

Mor a m= 1 6 ,nN^PC= 1 »BNFPf= 1 . 

° VP  C( 1 1=3. 644FA,.17,2.1716F-4,. 0839 
OFPri  1 1 =71  1 .8,  , ,20000.  , poopo  . .361 89.6 
OTA  M (11=1.2,17,1*1 
PFAV ( 21=2,3,17,1,1 

RFAY <3)=3, 4, 17*1*1 
OF  AM(  4 )=4, 6, 17, 1,1 
RFAV(5)=5, 6*17, 1,1 
of/im  {61=6,7*17.1  ,1 
oca  w (71=7,8,17, 1,1 
RFAVf 8 1=8,9,17,1,1 
ocaM{ 91=9,1 O, 1 7 ,1  ,1 
n f/\  m ( 10  1 = 10,11,17,1,1 
Q F A M ( 111  = 11,17,17,1,1 
n f a v ( 121  = 12*13*17,1,1 

PFAM(l3)=i3,]4,l7,i,l 

OFAY (141=1 4, 16*  17, 1,1 

PTAM( 16 ) = 1 8* 16,17, 1 , 1 

NFN=1  ,r;AT  = 6 »NT=200  , NOT  = 2 »DT=  .001 

NP (1*11=1.1 .26240. 

NP( 1 ,2 1 * l .1  .-87O0. 

MO| 2, 1 1 = 1 .1 .60480. 

MP( 2,21=1.1 .-16490. 

NP (3, 11=1, 1,6 04 80. 

MB ( 3*21=1 . 1 .- 1 6400. 
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Input  Data  - Concrete  Polygonal  Frame  - Dynamic  (Concluded) 


MP(A»1 )=1 ,2,40840. 
MP( 4*2 ) =1 ,7,-2967 n. 
NP( 8,1 )=1 ,2 .81200. 
MP( 8 »2 ) =1 ,2,-4  2940. 
m P { 4,1  1=1  ,? ,7]  200. 
NP‘(  6,2  1 =1  ,2  ,-42940. 
NP{ 7,1 1=1,3,15600. 
NPC  7,2 ) =1  ,7,-4  8010. 
MD( 8 ,2 ) =1 ,7 ,-87080. 
mp( o,2)=]  ,3  , — 97o8o. 


MP(  10,1  )=1  ,° 
M P ( 10,2)='  .7 

n p ( n ♦ i ) = i » 4 

MP( 1 1 ,7 1 =1  ,4 
NP(  12*1  1=1  ,4 
MP(  1 7 » 7 1 = 1 , 4 
NP<  1 3 » 1 ) = 1 ,4 
M P ( 13*21=1  ,4 
N& (14*11=1*5 
MP( 14,21=1 ,9 
M P ( 18,1  )=1  , * 
N O ( 1 5 » ? ) = 1 ,c 
MR ( 14,1 ) =1 ,8 

P'P(  16,21  = 1 ,8 


-1 8600. 
-4801 0. 
-31200. 

- 4 2 9 4 n , 

-31200. 
-42940. 
-40840. 
-29670. 
-90480. 
-16400. 
-80480 . 
- 1 64no  , 
-7*24'-'. 


P'P(  16,7  1 = 1 ,8,-8700. 

AT(  1 ) =0,AT  ( 7 1 = .002 3, AT  ( 3 1 = .0083  ,AT  (4  ) = .(U57,AT  ( o)  = .021  / 
MLP ( 1 ) = 5 , 1 . 

T(l,l»=0.,T(l,2)=.025,l(l,3)=.08,T(l,4)=.16,T(l,8)=.2b 
F<  1 , 1 > = ] . ,F(  1 *2  ) = . 545  ,FM  ,3  1 = . 135,F<  1 ,4  ) =0.  »F(  1 ,5)  =0. 

FKK=2,lSr>=l 
iroMP(4 ,7 , 9,9 ) =1,2 
*0|0K6  = 2 

TS( 7*1 1=1 ,7,6,7,8,17 
16(7,41=1,7,6,7,8,12  ‘ 

T 55  C 2,7  1=1  ,7,6,7,8,12 
!6( 2,101=1  ,7,6,7,8,12 
16(7,131  = 1 ,7,6,7,8,12 
1. 6 (7, 18)=1, 7, 6, 7, 8, 12 


I 


4.5.4  Results 


The  first  part  of  the  results  gives  the  first  10 
eigenvalues  and  mode  shapes.  Figure  48  shows  the  first  three  mode  shapes 
and  frequencies. 


Figure  48.  Concrete  Polygonal  Frame 
Eigenvalues  and  Eigenvectors 


H 

I 

I 

The  first  natural  frequency  of  the  reinforced  concrete  polygonal 
frame  was  50.9  rad/sec. 

The  first  natural  frequency  of  the  steel  polygonal  frame  was  56.5 
rad/sec.  The  natural  frequency  of  a beam  is  given  by 

fn  = n2  a/ET/AA 

The  members  of  the  steel  frame  had  properties:  E = 29  x 106  psi, 

I = 281.0  in4,  A = 10.6  in2and  A = 0.289  lb/in3. 

The  quantity  under  the  radical  is  equal  to 

SUM  = 5.158  x 10“ 

The  members  of  the  reinforced  concrete  frame  had  properties: 

E = 3.64  x 106  psi,  I = 36,189  in4,  A = 711.8  in2  and 
A = 0.0839  lb/in3. 

The  quantity  under  the  radical  is  equal  to 

v^I/AA  = 4.696  x 104 


The  results  of  problems  4.4  and  4.5  are  shown  to  be  consistent  since 
“st^conc 


and 

yfc- 1 / AA  S t _ i i g 

A\jh\  cone 


As  a further  example  of  the  results  consider  the  displacement 
response  history  at  node  7 as  seen  in  Figure  49.  Components  1 and  2 
correspond  to  the  horizontal  and  vertical  displacements,  respectively. 
The  horizontal  displacement  has  a period  of  0.125  second  which  com- 
pares well  with  the  expected  period  (from  first  node)  of  0.123  second. 


I 
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NORMAL I 7E0  FACTOR- 


TIME  l SECOND) 


.SiCJE+OG 


COMPONENT 

COMPONENT 


Figure  49.  Displacement  Response  History 
Concrete  Polygonal  Frame 


5.  Plates 


5.1  Square  Plates 

5.1.1  Load  Functions 


following  loads 


A simply  supported  square  plate  is  subjected  to  the 


Case  1.  Static  point  load  at  center 


Case  2.  Dynamic  point  load  at  center  with  the  time 
history  shown  in  Figure  50 

Case  3.  Static  uniform  pressure 

Case  4.  Dynamic  uniform  pressure  with  the  time  history 
shown  in  Figure  50 


0.002 


0.01  t,f 


Load  - Time  History 


F ( 1 b ) 


l/4(in) 

_L 


l/4(in) 

P(psi)  | 

Q 32SHJ.1. 


jure  50.  Square  Plate  - Static  and  Dynamic  Loading  (Continued) 


h. 


Figure  50.  Square  Plate  - Static  and  Dynamic  Loading  (Concluded). 


5.1.2  Finite  Element  Model 

Since  this  plate  is  symmetrical  with  respect  to  the 
x-axis  and  the  y-axis  with  origin  at  the  center,  one-quarter  of  the 
plate  is  sufficient  for  this  analysis.  Four  thin  plate  elements  were 
used  with  the  location  of  nodes  and  elements  shown  in  Figure  51. 


Figure  51.  Element  Assignment  - Square  Plate 


75 


HF  A D E R = * S T A T I C ANALYSIS  OF  SIMPLY  SUPPORTED  RECTANGULAR  PLATF 
NUMNP=9 »NFLTYP=1 ,LL  = ] * 

T X ( 1 ) =0  ♦ l » 1 *0  » 1 * 1 

tx(?)=oo.o,o»i,i 
T X ( 3 ) =0  ,0,0  » 1 * 1 , 1 
IX(4)=o,i»l,o,o*l 
rx(^)=n»o,Oin*o*l 

T X ( 6 ) =o,n,o , 1 »n»l 
TX(7)=1»1»1»0»0»1 
I x ( R ) = 1 »o, l ,o ,p » l 
T x ( 9 ) = l , o , l *i » n * l 
XV7T(  1 > = ?o. 

XY7T ( ?|=?n,,]0.! 

X v 7 T ( 3 ) = ?o. ,?n. 

X v 7 T ( 4 ) =1 o.  , 

XY2  T ( 5 ) =lo. , 10.  , 

XY7T ( 6 ) = 1 O • ,70.  , 

XYZT ( 7 ) = 

XY7T ( 8 ) = ,1 0. 

X Y7T ( 9) =,30. 

NPL ATE=4,PMDM= I . 

pmp I (1)  = .n0O7346,  , , ,3.376E7, 1.  125E7,  ,3.37bb7,  , 1 . 12*>E7 
' p I ATF< 1 4 ) = • ? 6 

•TPL  ATF  ( 1 ) =1  .7,6,4, ,1 
TP|  ATP  ( 7 1 = 7 ,3 ,6,6  , ,1 
T PLATE ( 3 ) =4,6 ,8 ,7, ,1 
T PL  AT--  '41=6, 6, 9, 8,  ,1 
OLMP( 3, 1) = , ,-76. 
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HEADER  = *DYNAV,IC  ANALYSIS  OF  A SIMPLY  SUPPORT  LD  PL  A I L 
NUMNP=9  ,NFLTYP=1 , NF = 1 0 »NDYN= 2 , 

T X { 1 )=o,i  .1  *0*1,1,  X YZ  T ( 1 ) =20  • , 

T X < 21=0,0,0,0,1  » 1 > XYZT(2>=20.,10., 

T X ( 3)=0»o,0»l»l»l*  XYZT(3)=20.,20., 

1 X ( 4 1 =0  *1  , 1 ,0 ,0  » 1 , XYZT( 4 ) =10.  , 

tx(si=o»o,o*o»o,i»  xyzt(5)=1o»»io«, 

TX(6)=0, 0,0, 1,0,1,  XYZT(6)=lO.,20.» 

TX< 7 1 =1 » l ,1  *0,0,1  , X Y Z T ( 7 ) = 

TX(8)=l,o,l»o,o,l,  XYZT(8)=0.»lo., 

TX(9)=1,o,l,l,o,l»  XYZT(o)=0.»20,» 

NPLATE=4»PNDM=1 , 

PPP  I ( 1 1 = .0007346, , , »3.375E7 , 1. 125E7, ,3.37bE7» ,1 . 125L7, 
PLATFd  4)  =.25, 

I PLATE ( 1 ) =1  ,2 ,5  ,4  , , 1 , 

I PL ATE ( 2 ) =2 ,3 ,6,5  * , 1 , 

TPLATE(7)=4,5,3,7,  ,1  , 

I PLATE (4  1=  5,6,9,8,  ,1  , 
tfpp=1 , N F N = T , M T = I ooo,NOT  = TO,nT  = .OOnl 
MP( 3,3 ) =1 , ,1 .0, 
w L P ( 1 ) = 4 » - 2 6 . 

T ( 1 ,1 ) =0. , T ( 1 ,2 ) =.002 , T ( 1,3)=. 01,  T ( 1 , 41=1.0 
F(l,  1)=0.,F(1,2)  = 1.0,F(1, 31=0. ,F ( 1 * 41=0., 

TC0MP(1  91=1,2,3,4,5,6 
K K'X  S = 2 

TS( 6,1  1=1  ,2, 3, 4, 6, 6 
(6(6,21=1  .2, 3,  A, 5, 6 
T S ( 6,31  = 1 ,2, 3, 6, 6, 6 
T 6 ( 6,41=1  , 2 , 3 , U , 6 , 6 

I 
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5.1.5  Input  Data  - Case  3 - Static  Uniform  Pressure 


HFAPFR=*5TAT ic  analysis  of  simply  SUPPORTED  RECTANGULAR  platl 

N!!MNP  = R »NFLTYP=1  »LL  = 1 > 

tx( i )=o,i ,i  ,n»i  *1 

T x ( ?)=n,n»n*n*i »i 

rx ( ^ ) =o»n,n»i , l »i 

rx(4)=n»i«l»n,n»i 

I X ( s ) =0  »0»0 ,0  »0  » 1 
T X ( 6) =0*0,0 , 1 ,0*  1 
T X ( 7 ) = 1 *] ,1 .0*9*1 

tx ( fi » = 1 »o» i »o  »n*i 
IX  < 9) =1 »n, 1 * 1 ,o» 1 
XY7T  ( 1 ) -?n, 

XY7T ( 2 1 =2".  * 10.  * 

XY7T( 7 ) = 2 n . *20. 

XY7T ( A ) - Lo. , 

XY7T ( S ) = 1 n. , 1 0. * 

XV7  T ( 6 * =1 n. *20.  * 

X Y7  T ( 7 ) = 

XY7T(8)=*10. 

XY7T ( 9)  = *20. 

NPLATE=4*PNDM= 1 . 

DFL^l  =1  • 

P^c.  I ( 1 > *.nOn7346. , ,,3.37&E7.1.12i>E7,  .3.3  75E7*  ,1.1?£>E7 
DLATFIl  41=.2F,~25. 

I PL  ATE ( 1)=1 »2*5*4**1 
I PL  A T E ( 2 ) = 2 » 3 » 6 » 5 * * 1 
TPLATF(3)=4*5,8»7»,1 
TDL  ATF  (4  1 =«j  ,6,9,8  * , 1 
Fl.M  ( 1 ) = 1.  , 


5.1.6  Input  Data  - Case  4 - Dynamic  - Uniform  Pressure 


i 


. 


i 


I 


HFADFR=*DYNAMIC  analysis  of  a simply  supported  PLATF. 

NUVNP  = 9 ,NFLTYP=1 , NF  = 1 0 » NDYN= 2 * 

TX(l)=0»l,l»n»l,I»  X YZ  T ( 1 ) =20 • * 

TX<2)=0»0,0»0»1»1»  XYZT ( 2 > =20. » 10. » 

TX(3)=0»0,0»1»1*1»  XYZT ( 3 ) =20 • » 20  • » 

TX(4)=0>1»1 ,0*0,1,  XYZT (4 1=10.. 

IX(5)=0»9»0»0»n,l»  XYZT<5>=10.,10.» 

TX(6)=0,n»0»l,0»l,  XY7T(6)=l0.,20.» 

T X ( 7 > =1 ,1 ,1 ,0,0*1 * XYZT ( 7 ) = 

IX(8)=l*n,1.n»n,l,  XYZT ( R ) =0. , 10. » 

TX<9)=1, 0,1, 1,0,1,  XY7.T(  9 ) =0.  * 20.  » 

NPLATE  = 4»PND^= 1 , 

PMP  I ( 1 ) =.^007346,  , , ,3.37  5F  7.1.1 2 E 7*  ,3.37b  E7  , ,1.12  5E  t ♦ 
PLATE  I 1 4)  = .25, 

T PLATE! 11 =1 ,2 ,5,4,  ,1 
TPL ATF ( 2 ) =?  ,3 ,6,5,  , 1 
T n>L  ATE  ( 7 ) =4 ,5 ,8 ,7  , , 1 
IDEATE! 41=5,6,9,8,,! 

7 FPP= 1 ,NFN  = 1. » N T= ] non , NOT  = 1 0,DT  = .nonl 
NP( 1 ,3)=1, ,-62.5 
MP( 2»3)=1, ,-125. 

NP ( 3 » 3 ) =1  , ,-62.5 
NP( 4,3) =1 , ,-125. 

ND( 5 ,3 ) =1  , ,-750. 

NP( 6 , 3 ) = 1 , ,-125. 

MP( 7 , 3 1 = 1 , ,-62. s 
MPI 8,31=1  . ,-125. 

MD(  o,?  t=i , , -62 . 6 

NLP ( ] ) =4, l . 

T(7  ,1  ) =0. ,T ( 1 , 2 1 = .00  2 ,T ( 1 ,3 1 =.01  , T ( 1 ,4) = 1 .0 
Fn,l)=0.,F(t»2>  = 1.0»FH,9)=0.,F(l,4)=0. 

V «K  =2 

ICOMP (1  97=1,2,3,4,5,6 
PK-<S  = 2 

T S ( 6 , 1 ) = 1 , 2 , 3 , 4 , 6 , 6 
!5(6,2)=1 .2, 3, 4, 5, 6 
T S ( 6 , ^l = 1 ,2*3, 4, 5, 6 
T c- (6, 4 1=1  ,2,3, 4, 5, 6 


I 
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5.1.7  Results 


Case  1 - Static 


SAP  IV 


Maximum  deflection 
(at  center) 

Bending  moment  at 

x = 15  in,  y = 15  in 


Thin  Plate  Theory 
-i  (4) 


-0.3990  x 10"1  in  -0.4223  x 10  in 


Mxx  - Myy  - -’-136  "xx  ■ "yy  = -°-9831 


x 10  lb-in 


xx  yy 

x 10  lb-in  ^ 


(Resultant  in  Element  2) 


Case  3 - Static 


Maximum  deflection 
(at  center) 


-5.2135  in 


Maximum  bending  moment  1481  lb- in 

(Resultant  in 
Element  2) 


-6.0484  in 

1916  lb-in 
(at  center) 


Fundamental  Frequency  for  Symmetric  Model  (rad/sec) 


Mode  Number 


SAP  IV 
0.1936  x 103 
0.9469  x 103 
1.358  x 103 


Displacement  - Time  Histories 


Thin  Plate  Theory 
0.1908  x 103 
0.9539  x 103 
1.717  x 103 


The  response  of  the  center  of  the  plate  for  the 
two  dynamic  load  functions  is  given  in  Figures  52  and  53. 


Timoshenko,  S. , and  Woinowsky-Krieger,  S. , Theory  of  Plates  and  Shells, 
Page  143-149. 

3Szilard,  Rudolph,  Theory  and  Analysis  of  Plates,  page  650. 


5Clark,  S.  K. , Dynamics  of  Continuous  Elements,  page  172. 
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DYNAMIC  ANALYSIS  OF  fl  S 
NORMALIZED  DISPLACEMENT 


TIME  (SECOND) 


COMPONENT 


NORMALIZED  FACTOR 


COMPONENT  2 
COMPONENT'  3 


Figure  52.  Nodal  Displacement  of  Plate  Center 
Point  Load 


DISPLACEMENT 


Case  1 Case  3 


Material:  Steel 

Members:  40  In  x 40  in  x 1/4  in  plate 
2 in  x 2 in  x 40  in  box  beams 
with  a wall  thickness  of  1/4  in. 


Figure  54.  Square  Plate  with  Two  Reinforcing  Ribs 
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5.2.2  Finite  Element  Model 


Due  to  symmetry,  a quarter  of  the  plate  is  considered 
for  this  analysis.  A group  of  four  thin  plate  elements  and  four  3-D 
beam  elements  is  used  with  the  location  of  nodes  and  the  elements 
shown  in  Figure  55. 


Numeral  - node  number 

(3  - plate  ele- 
ment number 

A-  beam  ele- 
ment number 


Figure  55.  Element  Assignment  Two-Ribbed  Plate 
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Case  1 - Static  Point  Load  - at  Center 


HFADFR=*STATIC  PLATE  WITH  REINFORCING  RIbS 
NUMNP  = 9 »NFLTYP=2»LL=1 

t x ( l ) =o , i »i *o  * l *1* 

t x ( 2 1 =0 ,o « 0 ,0  » 1,1, 

TX( 3 ) =o ,o,o » 1 » 1 » 1 » 

I X ( 4 ) =o  ,1 ,1 »o»o , 1 , 

T x ( *> ) =n ,o ,o ,o ,n » l * 

T X ( 6 ) =o  ,n , o , 1 , o , ] , 

T X ( 7 ) = T *1,1  *0,0,1, 
fX( 91=1,0,1 ,0,0,1, 

T X ( O ) = 1 , 0 , 1 , 1 , 0 , ] , 

XY7T  ( 1 1 =?n. 

XYZT  ( 2 1 =2o. ,10. , 

XYZT  ( 3 1 =2o. ,20, , 

XYZT(4)=ln. , 

XY2T ( 5 ) =lo. , 10. , 

XY7T(6)=ln.  ,20.  , 

X Y7  T ( 7 ) = , 

XY7T  ( 8 1 = , 1 0.  , 

V V7  T ( 9 ) = ,70.  , 

MRFAM=W  ,Rr>irpr=  ] , RNMPC-  1 , 
d M P c ( 1 )=3.r7,. 333, .00073 6 6 
RFPCll  )=. 799,,, .0395, .RR-3, .383 
RF  AM ( 1 1 =1 ,2  ,4 , 1 ,1  » 

RFAM (21=2,3,4,1,1, 

R F A M ( 3 1=3,6,4,1,1, 

R F A M ( 4 ) = ',,9,4,1  ,1  , 

N|PL  A T E = 4 , P NDM=  1 , 

PMP  I ( l)=.ooo7346»,  , » 3 • 37 5 F t ,1  • 1 2 ‘J  E f , »3.37oF7  ,,1,12pE7 
PLATFU  4 1 = . 2 5 
T PL  A TF ( 1 1 = 1 ,2,6,4,  ,] 

IP|  ATE( 21=2,3 ,6,5,  ,1 
JPLATF (31=4,6,8,7,  ,1 
TPLATE(4)=5,6,9,8,  ,1 
CLMD  ( 3,1  ) = , ,-26. 


5.2.4  Input  Data  - Case  2 - Dynamic  Point  Load  at  Center 


HF  A D E R = Y N A M I C ANALYSIS  OF  A SIMPLY-SUPPORTED  PLATE  w 2 REIN 
NUMN,P  = Q ,NELTYP  = 2»NF  = ln  »NDYN=? 

TY( 1 )=o.l.l *0*1 »1»  XY7T ( 1 ) =20 , 

TX (?) =n  ,0.0*0  « 1 * 1 , XYZT ( 7 ) =70.  ♦ lo.  , 

rx(3)=n*n,o»i*i»i»  xyzt ( 3 ) =20.  *20. » 

TX(A)=0»1,1,0»0»1*  XYZT (A ) =10. > 

TX(A)=n»o,o,o,o»i*  XYZT(5>=10.»10.* 

IX(fi)=0*n.0»]»n.l,  XYZT(6)=10.»20.» 

TX  ( 7) =1 ,1 , 1 ,n *0»1 » X Y7  T ( 7 ) = 

TX(‘U=1*n.1»0»n,l»  XYZT(8>=0.,10.» 

TX{9)=1»0.1»1*0»1»  XYZT ( 9 ) =0. *20. * 

r>!RP  AM  = A ,RMF PC=  1 *PNMPC=  1 » 

pupr (1  ) = 3.E7» .333  » .0OO73A6 

"PPC ( 1 ) =.705*  » ».039S».383».383 

3 F A M ( 1 ) = l ,7  ,A»  1 »1 

RFAM( 2 )=2.3,A»1 »1 

RFAf/(  3)=3*6»A.  1 * 1 

RFAf/  ( A ) =6 ,0»  A t ] .1 

M PL  A TE  = A » PNOM= 1 

pyp  I ( 1 ) = .o0073A6»  . » »3.37£>E7*  1 . 12PE7»  » A.3  75E7  » ,1  . 125E7 

PLATEU  A t a . 7 5 

T PLATE!  1)=1 *2»3»A»»1 

TPLATE(2)=2»3,6»5» *1 

TPLATF ( 3 ) =A,S  .8 ,7*  ♦ 1 

IPLATE (A) =S»S»9»8»  » 1 

T F P P = 1 » N F M a 1 * N T = 1 000  * NOT  = 1 0»DT=.0001 
md ( 3,3 ) =1 . , 1 . 

A'LP(  1 ) =A  »-?S. 

T( 1 ♦ 1 ) = 0 • * T ( 1»2)  = .00?*T ( 1 » 3 ) = • 0 1 » T ( 1 » A ) = 1 • 

F(  1 »1  ) =0. »F ( 1 »?)  = 1 • »F(  1*3)=0.»F ( 1»A)=0. 

X K K = 7 

TfOMP ( 3 » 5 1 =3 
KKK  S = 2 

I F ( 2 * 2 ) = 1 . 7 , 3 , A , S , 6 , 7 , 8 , o , ] o , 1 1 , 1 2 
TMS»'M=1.?.3»A,5»6 


5.2.5  Input  Data  - Case  3 - Static  Uniform  Pressure 


headfk=*stat ic  plate  with  reinforcing  ribs 

NUMNP=9,NFLTYP=2»LL=1 

r x ( i i = o*i ,1 *0*1*1 » 

T X ( 2)=n»o«o*0*l *1  * 

TX(3i=o»o,o*l»l*l * 
tx ( 4 ) =n,i * i ,n,n» 1 * 

T x { s i =o ,o ,o  ,o *o » 1 * 

TX(M=0*n*0,l»O*l  * 

T X ( 7 ) = 1 *1*1  ,0*0*1  * 

IX ( 8 ) = 1 *0, 1 *0*0* 1 » 

IX(9)=] *0,1,1 *0*1 , 

XYZT ( 1 ) =2o. 

XYZT ( 2 ) =2o. , 10*  * 

XV7T ( 3 ) =2o. ,20. , 

XY7T(4)=7o., 

XYZTIS )=lo.*10. * 

YV7T ( 6 ) =1 o. *20. * 

XYZT ( 7 1 =, 

YY7T(81=,10., 

XY7TI91 =*20. , 

NRFAM  = 8 ,rmtpc=2 ,RNMPC= 1 * 

RMpri  ] 1=7. r 7*. 3 37,  .000  77 /.A 
R FPC( 1 ) =.798, , ,.0395,.38?*.38i 
PFPC(?)=l.,:'Q,*,.07o*.76A*.766 

noAM  ( 1 ) =1  ,2  .4  * 1 * 1 * * 

PEAM(3)=3*A*4*1*1 »**,»**3 
p F A v (81=2*4*7*1  ,2*»**»*»3 

PFAM(7)=4,5,7»1*2*»*»»»*1 
roam(81=5*A,7,1 ,2* 

MPLATF=4»PN0M=1 

PMPI  ( 1 1 = .0  00  734  6,  , , *3.37SF7*  1 . 125F.7*  ,3. 375F7  , , 1 . 122E7  , 

T D L A T E (1  l = 1 *2  ,S  ,4*  ,1  , 1 

Of.  A TO  ( l 4 ) = . 2 S 

TPL ATE ( 1 1 = 1 *2  *5 ,4* ,1  * 1 

1 PLATE ( 3) =4  *3  *8 , 7,  , 1 , 1 

IDL ATL ( 4 1 =S  *6  *9  *8  * , 1 

CLMD ( 3 *1 1 = * ,-2b. , 


5.2.6  Input  Data  - Case  4 - Dynamic  - Uniform  Pressure 


hfader=*:>ynamic  analysis  of  a simply-supported  plate  with  2 
NUMNP  = Q »nfltyp=2 ,NF  = 1 n ,NDYN=2» 

T X ( 1 )=O»l,l,0»l»l»  XYZT(1)=20.» 

TX<2)=0»o,n,0»]»l»  XYZT(2)=2O.,10.» 

TX(3)=0»0,0»1»1»1»  XYZT  ( 3 ) = 20  » » 20  • , 

ix  ( 4 ) = 0 ,1 ,1  ,0  ,0,1  , XYZT(4)=lo., 

TX(5)=n»o,o»o»o»i»  xyzt(6)=io.»io., 

IX(6)=0»o»0*l»o,l,  XYZTI6 ) =lo.  ,20. , 

TX(  7) =1 ,1 , 1 ,0 ,0 , 1 » X YZ  T ( 7 ) = 

TX(B)=l,o,l,o*n,l,  XYZT(8)=0.»10. » 

T X ( 9 ) = 1 » o , i , 1 » o » 1 » XY7T(O)=0.,2o., 

MPFAM  = 4 ,c»MMpr=  1 »rnppc  = 1 , 
n.MPr  ( 1 1=3. F7, .333 ,.000 7346 
n^PCI 1 ) = • 795 , » » • 0395 , •383*. 383 
PFAM(1 1=1  ,2*4*1 *1 
RFAVf  2 ) =2 ,0  *4*  1 *1 
OPAM(3)=3, 6*4*1  ,1 
PFAM(4)=6,9,4»1*1 
NPLATE=4*PNDM= 1 

PMPI ( 1 ) =.oon7346, » , ,3.376E7,1 . 125E7* ,o.37SE7» , 1 .125E7 
PLATO ( 1 4 i = . 2 5 
T PL  A TF ( 1 1=1  ,2,6,4,  ,1 
T PL  A T F ( 71=2,3,6,5,  ,1 
T PL ATF ( 31=4, 6, 8, 7, ,1 
T PL  A TF (41=5, 6,9*8, *1 

TFPR=1 , NFN=1 ,NT=1 00O,N0T=l0,DT=.O0Ol 

M P ( 1 , 3 ) = 1 ,,-62.5 

A'P(?,3)=1,,-125. 

NP( 3,3  1=1  , ,-62 .6 
N P ( 4,3 ) =1 , ,-125. 

NP( 6 ,3 ) =1  , ,-250. 

NP(6,3)=1*,-125. 

MP( 7 , 3 ) = 1 , ,-62*6 
MP( R ,3 ) =1  , ,-126. 

MO( P,3 ) =1 , ,-62.6 
A'L°  ( 1 ) =4, 1 . 

T(1,1>=0.,T(1,2)=.C02»T(1,3)=.01*T(1,4)=1.0 
o( I ,1  I =0.  , F ( I , 21  = 1.0, F( 1 , 3 ) =0 • , F ( 1,41=0., 

/ r K = 2 

I COMP ( 3 , 5 1 = 3 , 

FKK  S=2 

T S <2, 2)=1, 2,3* 4, 5,6,7, 8, 9, 10, 11, 12, 

TS (6, 2)=1,2*3,4, 6, 6, 


BHBMMmMBB 


5.2.7  Results 


Case  1 - Static  - Concentrated  Load 


Maximum  Deflection  (at  center) 

Maximum  Bending  Moment 
(at  center  in  beam  element) 


2.695x10  in 


2.347x10  lb-in 


Case  3 - Static  - Distributed  Load 


Maximum  Deflection  (at  center) 

Maximum  Bending  Moment  (at  center 
in  beam  element) 


2.526x10  lb-in 


Plots  of  the  deflection-time  histories  for  the 
dynamic  loads  are  given  in  Figures  56  and  57. 


DISPLRCEMENT 
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DYNAMIC  ANALYSIS  OF  R SIMPLY-SUPPORTED  PLATE  WITH 
NORMALIZED  DISPLRCEMENT  RESPONSE  RT  NODE  3 


Figure  57.  \tr  itory  of  Plate  Center  - Uniform  Load 


5.3  Square  Plate  with  Six  Reinforcing  Ribs 
5.3.1  Load  Functions 

A simply  supported  square  plate  with  six  reinforcing 
ribs  (Figure  58)  is  subjected  to  the  same  loads  as  problem  5.1. 


Material:  Steel  Members:  40  in  x 40  in  x h in  plate 

2 in  x 2 in  x 40  in  box  beams 
with  wall  thickness  of  h in 

Figure  58.  Square  Plate  with  Six  Reinforcing  Ribs 


5.3.2  Finite  Element  Model 


Again,  only  one-quarter  of  the  plate  is  necessary  for 
this  analysis.  A group  of  four  thin  plate  elements  and  eight  3-D  beam 
elements  is  used  with  the  location  of  nodes  and  the  elements  shown  in 
Figure  59. 


Figure  59.  Element  Assignment  - Six-Ribbed  Plate 


• f 


5.3.3  Input  Data  - Case  1 - Static  - Point  Load 


HFADt.  R = *STATIC  PLATE  aITH  REINFORCING  K I OS 
NUVMP  = 9 »NFLTYP=2»LL=1 
TX(  1 ) = o » 1 , 1 . o j ] > 1 , 

TX( 2 ) =0»0.0,o, 1 , 1 , 

TV { 0 ) =0,0.0 , 1 , 1 , 1 , 

T x ( A > = n * l ,1 ,0*0,1, 

T X ( ci)=o, o,0, 0,0,1  , 

T X ( M =0»o,0  » 1 ,0 , 1 , 

TX( 7)=1  ,1  ,1  ,0,0,1 , 

T X ( 8 ) = 1 ,0,1 ,0,0,1, 

T X ( o ) =1  ,o,  1 , 1 ,o  , 1 , 

XY 7 T( 1 ) = ?n. 

X Y7  T ( 2 ) =2'-'.  , 10.  , 

XY7T(  3 5 =2".  , 20  • * 

XYZT(4)=lo.  , 

XY7T  ( * ) =1  o.  , 1 0.  , 

XYZT(6)=lo.,20., 

X Y7  T ( 7)  =, 

X Y7  T ( 8 ) = , 1 0 . , 

XV7  T ( 9 ) = , 20.  , 

NRF  AM  = # ,RMFP(-  = 2 ,RNMP0=1  , 

BMPC<1)=3.F7,.333,. 000 7346 
RFPC( 1 > = .793,  , , ,039b, .383, .383 
RFPC(2)=1.59», , • 079 , • 7(36 , • 766 

pFAY( 1 ) =1 ,2  ,4 , 1 ,1  ,,,,,,,  1 
RFAM(3)=3. 6,4,1  ,]  ,,,,,,,0 
pfam(5)=2,5,7,1  »2»»»»»»»3 
□ CAM( 7) = 4,3 ,7, 1 ,2  » , , , * , * 1 
PFAU(8)=3,6,7,1,2» 

MPL A TE  = 4 , PMHM= 1 

DMPI  ( 1 >=.00073  4 8,,  ,,-3.3  75b7,l.l2bE7,,3.3  75E7,,1.12t»E7, 
PFLML=  1 . , 

PL A TF ( 1 4 > = .28, -2*. 

TPLATF ( 1 ) =1  ,2,5,4  , , 1 , 1 
T PLATE ( 3 ) =4,5 ,8,7,  , 1 , 1 
T PL  ATE ( 4 ) =5 ,6 ,9 ,8  , , 1 
TPLATE ( 4 ) =5, 6, 9, 8,  ,1 , 

FI  V(  i ) = l.. 
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5.3.4  Input  Data  - Case  2 - Dynamic  - Point  Load 


HFADEK  = *DYNAMIC  ANALYSIS  of  a SIMPLY  SUPPORTED  PLATE  VV  REIN  KlbS 
NUVNP=9 »NFLTYP=2*NF=1o ,NDYN=7» 

TV( T )=0,1  ,T  ,0*1*1*  XY7  T ( 1 1=20.  * 

TX(7>=o*o,o,o*l,l,  XYZT(2)=20.*10., 


TX(3)=0»0,0*1»1*1»  XYZT(3)=20.*20., 

TX(4)=0*l*l»O»0»l»  X YZ  T ( 4 ) = 1 0 • » 

ix(6)=o»n,o»o*o»i»  xyzt(5)=io«,io«* 

TX(6)=0,n,o, 1,0*1*  XY7T(6)=l0.*7o., 

TX ( 7) =1  *1 .1  *0*0*1  * X YZ  T ( 7 ) = 

TX{B)=1*o*l, o*o»l*  XYZT ( 8 ) =0. » 1 o. , 

IX(9)=1»0,1,1*0»1»  XYZT(0)=0.»70., 

mrFAM=8»3MFPC=2»RNMPC=1 


RMPf ( 1 1 =3.F7, . 333* .0007346 
RPPr(  1 )=.795, , , .0393,  .383*.  383 
ROPC ( 2 1 =1 . SO, , , . 070**766  *.766 
R F A M ( 1 1=1, 7*4*1  *1*»*»*»*1 

0 r A M ( 3 ) = 3 » 6 » 4 , 1 » 1 *******  3 
PF AM ( 5 ) =2 , 5 , 7 , 1 , 2 * * * * * * * 3 
P F A m ( 7)=4, 6*7*1 *2*  * , * * * * 1 
R P A M (81=6*6*7*1 *2 

NPL  ATE  = 4 » PNDM  = 1 

PMP I ( 1 ) =.n00  7346, , , ,3.37bE7, 1. 12bE7* ,3.375E7,  ,1 . 12&E7, 
PLATE ( 1 ,3  *4 ) = . 26 
T PLATE ( 1 ) =1  ,2, 5 *4, ,1,1 
1PLATF(3)=4, 5*8*7* *1*1 
TPL ATL (4 ) =6 ,6,0,8  , , 1 

1FPR=1  *1*  rA'  = l , NT=  1 000  » NOT  = 10,DT=. 000 1 
MD(  3,3 ) =1 , , 1 
NLP (1  ) =4,-26  . 

T(1»U=0.,T(1»2)  = .002»T(1»3)=.01»T(1*4)=1. 
p(l  ,1  )=0.,F(  1,2)=1  .0,F(1,3)=0.,F(] » 4 ) = 0 • 

<KK  =2 

1 COMP (1  91  = 1,2*3*4,5,6 
r K K S = 2 

IS( 2*1 1=1 *2*3, 4, 6, 6, 7*8*9*10, 11, 12 
16(2,21=1,7, 3, 4, 6, 6, 7, 8, 9, m, u, 12 
I6(2»3)  = l, 2, 3, 4, 6, 6*7*8, p,10, 11, 12 
T6(?,4l=l  ,7,3,6,6,6,7,8,0,10,11  ,12 
16(7,61=1,2*3,4,5*6*7*8*0,10,11,12 
T 6 (?»6)=1, 7,3, 4 *5, 6*7, 8*9*10*11*12 
TS(7»  7 ) = 1 , 2*3*4, 5*6*7*8.9,10, 11*12 
IS (7*81=1, 7, 3, 4*5*6*7*8,0*10, 11, 12 
16(6*1  1=1  ,7*3, 4, 5*6 
16(6,71=1  .7*3, 4, 6, 6 
16(6*3  1=1  ,7*3*4, 5*6 
t6(6*41=1,?,3*4,5*6 
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5.3.5  Input  Data  - Case  3 - Static  Uniform  Pressure 


- 


t 


HEAD FR=* STATIC  PLATE  WITH  REINFORCING  RIBS 

numnp=q,nfltyp=?,ll=] 

T X ( 1 1 =0  » 1 » 1 » 0 » 1 » 1 » 

t x ( 2 1 *o  «o,o ,o , ] Hi 
IX ( 3) =0 ,0,0 , 1 , 1 , 1 , 
r x ( /» ) = o » i iiiOiAiii 
TX(  5 )=0,O.0,0,0,l  , 

I X ( 6 1 = O , n,0, 1,0*1, 

T X ( 7) =1  ,1 ,1 ,0,0,1  , 

T X ( 8 ) = 1 ,0,1  ,0,0,1, 

IX ( 9) =1  ,0,1  ,1,0*1  * 

XY7T(  1 1=  >o. 

XY7T ( 2 1 =7o.  ,10., 

XY7T  ( 3 ) =?0.  ,20. » 

XY7T ( ft ) = 1 o , , 

XY7T ( 5 ) =1 o.  , 10.  , 

XY7  T ( 6 ) =1 o.  ,20.  , 

XY7T  ( 7) =, 

XY7.T  ( 8 ) = ,1  0.  * 

XY7T(9)  = *?0.  , 

MnFAN=A»nNrPC=l  , RNMPC=1  * 

PMPC(  1 )=3.F7,.333 ,.0007346 
RPPr(  1 >=.798,  , ,.0395, ,3 83, .383 
PPAM(l)=l,?,ft»] ,1, 
n F A M ( 21=2,3,4,1  ,1  , 

PFAM (31=3,6,4,1,1, 

Q F A M ( 4 1 =6,9,4 , 1 , 1 , 

MPLATF=ft,DNDM= 1 , 

pwpI(l)=.r)0073ftfc,,,,3»37bL7,1.122L7,,3.37bt7,,l.]?5t7 
T PL ATF(  1 1=1  ,2,5,4,  ,1 
T PL  ATE (21=3,2,4,5,  ,1 
IPL ATE (31=4,5,8,7,,! 

I PL ATE (4) =8,6,9 ,8,  ,1 

DtT|  V[  = 1 . , 

DLATrn  ft  i = .25,-25. 

CLV( 7 ) = T . , 
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5.3.6  Input  Data  - Case  4 - Dynamic  - Uniform  Pressure 


HEADER=*DYNAMIC  ANALYSIS  OF  A SIMPLY  SUPPORTED  PLATE  .v  I T H REIN 
NUMNP=9  »NELTYP  = 2»NF=10  »NDYN=2» 

IX ( 1 ) =9  , l ,1  ,0, 1 * 1 » XYZ  T ( 1 1 =20  » , 

1X(?)=0»9,0»0»1»1»  XYZT ( 2 ) =20* , 1 0 » , 

TX(3)=0,9,0,1»1»1»  XYZT(3)=20,»20., 

!X(4)=o,l.l »0>n*lt  XYZT(4)=10., 

ix(5)=o»o,o,n»o*i»  XYZT(6)=io.»io.» 

!X(6)=0»o,n,l,0»l»  XYZT(6)=l0»*20«* 

TX(7)=l.],l.n^,1»  XYZT  ( 7 ) = 

lX(8)=l,o.i,o,o,l,  XYZT ( 8 ) =0. * in.  , 

!X{9)=1»0,1,1,0»1»  XYZT  (91=0,  ,20.  , 

MPFAM  = fl ,pmpdc=2 ,RMPPC=  1 , 

RVPr<l)=3.F7*. 3333*. 0007*46 

RFPC ( 1 1 =. 796,0. »0. , .0395 , .383, . 383 

R F P (*  ( 2 1 = 1 .59,0.  ,0. , .079, ,76 6, .766 

RFAM(  1)=1,2«4,1,1,,,****1 

REAM(3)=3,6,4,1,1,,,,,,,3 

REAM  (5)=2,5,7,1  ,2  ,3 

REAM { 7 ) =4, 5 , 7, 1 1 

REAM (81=5,6,7,1  ,2, 

NPLATE=4»PNDM= 1 , 

PMP I ( 1 1 = .0097346, ,,,3.375K7,1.125E7,,3.375E7,,1.125E7, 

9 1 ATFI  ] ,3, 4 1 = .25, 

T 9|  A T F ( 1 1=1  ,7,6,4,,!  ,1, 

T pl  A T F ( 3 1 = 4 , 6 » 8 » 7 , , 1 .1  , 

TDLATE(4)=6,6,9,8,,1  , 

TEPR=1  ,MFa)=1  ,NT  = 1 onn,NOT  = 10,DT  = .0nnl  , 

MP( 1 ,3)=1 , ,-62.6, 

NP( 2 ,31=1  . ,-l 25. , 

WP( 3,31=1 , ,-62.5, 

NP( 4,3 ) =1 , ,-125. , 

NP( 5 ,3 1 =1 , ,-259. , 

NP(6,3)=1 , ,-125. , 

MP( 7,31=1, ,-62.5, 

MP( 8 ,3 1 =1 , ,-125. , 

mo (9, 3 1=1 , ,-67.5, 

MLP ( 1 ) =4 , 1 . 

T(1  ,1  )=0.,T(] ,?)  = .00?,T( 1 ,3)=.01  ,T(1,4)=1. 

F ( 1 ,1  1=9., F(1  , 7 ) = 1 • 0 , F ( 1 ,31=0. »F( 1,4)=0. 
k"  X K - 2 , 

T COMP { 1 91=1,2,3,4,5,6 
KKK  S=2  » 

15(2,11=1,2,3,4,5,6,7,8,9,19,11,12, 

TM?,2)=1,2»3»4, 5*6,7,8,9,10, 11, 12, 

TS( 2,31  = 1, ?, 3, 4, 5, 6, 7, 8, 9, 19, 11, 12, 

16(2,41=1, 2, 3, 4, 5,6, 7, 8, 9, io, 11, 12, 

TS(?,5)=1,7, 3, 4, 6, 6, 7, 8, 9, in, 11, 17, 

1-6(7,61=1  ,7,3,4,5,6,7,8,9,10,11  ,12, 

T 5(7,71=1, ?, •»,/,, 5, 6, 7, 8, 9, lo,H, 12, 
16(7,81=1,2,3,4,5,6,7,8,9,10,11,12, 

16(6.,  11=1,7,3,4,5,6 
16(6,21=1,7,3,4,5,6 
16(6,31  = 1 ,7,3, 4 , 5 , 6 
I6(6,4)=l,?,3,4,5,6 


5.3.7  Results 
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Case  1 - Static  - Concentrated  Load 
Maximum  Deflection  (at  center) 


-0. 1454x10” 2 in 


Maximum  Bendinq  Moment  (at 
center  in  beam  element) 


= 1 . 637xl02  lb-in 


Case  3 - Static  - Distributed  Load 

Maximum  Deflection  (at  center)  = -0.19624  in 


Maximum  Bending  Moment  (at  node 
2,  or  5 in  beam  element) 


= 1. 585x1 0*  lb-in 


The  displacement  responses  at  the  center  node  for  the 
two  dynamic  forcing  functions  are  shown  in  Figures  60  and  61. 
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DYNAMIC  ANALYSIS  OF  A SIMPLY  SUPPORTED  PLATE  W RE  I 
NORMALIZED  DISPLACEMENT  RESPONSE  AT  NODE  3 


TIME  (SECOND) 


COMPONENT  1 


NORMALIZED  FACT HR 


COMPONENT  2 
COMPONENT  3 


Figure  60.  Nodal  Displacement  - Plate  Center 
Point  Load 


DISPLACEMENT 


DYNAMIC  ANRI  YSIS  OF  A SIMPLY  5UPP0RTFD  PLATE  WITH 
NORMALIZED  DISPLACFMFN1  RESPONSE  AT  NODE  3 


NORMALIZED  FACTOR:: 


TIMF 
.35F  01 


(SECOND) 


COMPONENT 

1 

COMPONENT 

2 

component" 

3 

Figure  61.  Nodal  Di splacement  Plate.  Center  - 
Uniform  Load 
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5.4  Thick  Rectangular  Plate  Analyses  Utilizing  Thick  Shell  Elements 


5.4.1  Load  Functions 

A simply  supported  rectangular  plate  is  subjected  to 
the  following  loads: 

Case  1.  Static  point  load  at  the  center  o^  the  top  face,  and 

Case  2.  Dynamic  point  load  at  the  center  of  the  top  face 
with  the  time  history  shown  in  Figure  62. 

Case  3.  Dynamic  uniform  pressure  on  the  top  face  with  the 
time  history  shown  in  Figure  62. 


60  in 


60  in 


40  in 

Point  Load 


Y 40  in 

X Distributed  Load 


Figure  62.  Thick  Rectangular  Plate  - Static  and  Dynamic  Loading 

101 


-i. 


5.4.2 


Finite  Element  Model 


Due  to  symmetry  one-quarter  of  the  plate  is  adequate 
for  this  analysis.  A group  of  four  thick  shell  elements  is  used 
with  the  location  of  nodes  and  elements  shown  in  Figure  63. 


z 


Figure  63.  Element  Assignment  - Thick  Plate 


102 


#-  c> 

— H r— H 

—I 

•> 

-H 

•« 

— • 

* 

r-4 

H H 

— t 

r-H 

•H 

i— 

D o 

r-H 

D 

r-H 

•>  #■ 

* 

•> 

•» 

.J 

o 

o 

r-H 

o 

X 

•*  •» 

#> 

•> 

UJ 

r-H  r-H 

o 

o 

-1 

-jJ 

J 

II 

a 

— 

j 

-* 

• 

#> 

r-H 

X) 

fNJ 

-3- 

* 

“O 

CO 

UJ 

o 

-4- 

O 

1— 

_l 

o 

_J 

4-> 

a_ 

#> 

C 

r— H 

o 

v 

Q_ 

u 

*— • 

O' 

• 

X 

m 

O 

1— 

•> 

X 

4-> 

'O 

<o 

UJ 

•> 

+-> 

h- 

f} 

(/) 

X 

f 

1 

0 

»• 

CL 

n 

r— 

CL 

-o 

0) 

t/) 

O 

LO 

-ot 

fO 

o 

> 

#> 

-J 

r-H 

1 

a 

II 

T 



X 

o 

(O 

wO 

CM 

X 

o 

#> 

e 

-*-> 

3 

OF 

AJ 

V. 

CL 

#> 

•> 

C 

# 

f\ 

CM 

i-h 

•— * 

H 

AJ 

II 

V/) 

II 

r 

» 

CO 

X 

J 

AJ 

Is- 

-J 

— ! 

l\J 

(' 

f. 

*d“ 

< 

*- 

II 

•> 

•» 

•> 

• 

r~ 

o 

C\i 

if) 

< 

II 

(\j 

AJ 

C<^ 

if 

a 

r- 

» 

•* 

U 

X 

X 

CP 

r— H 

•— « 

1— 

C 

— i 

fA 

if 

_j 

r— 

•> 

« 

♦ 

< 

Lu 

♦ II 

r- 

vC 

c 

K- 

r*  — 

r-H 

If 

* 

* (t~ 

#■ 

* 

* 

if 

vC 

f 

*-H 

c 

II 

if 

r » 

r-H 

A. 

r 

* 

II 

^ oc 

II 

#• 

# 

uj  a 

♦ 

vO 

c- 

. c 

i: 

r-  IX 

O 

r-H 

r-H 

r- 

< 

s 

w W 

w 

w 

w 

u 

X 

X X 

X 

X 

X 

X 

X T 

►— 

•— 

• • • if  • • 

*-  if  c • c c 

♦ r-  C\J  ► r—i  f\ 


••  • h C C » if  c c 

• C * «»H  f\  f-HOsJ 

•>  ♦lOiTiriririririr 


• • if  tf' 

r— < r—  • • 

*••>*•>  ••  • 

• • if  • • • • 

ir  c • c c • ir  c 

«—  (V  * c— « f\  if  r-  Aj 


• • • • • • •ccc\r\f\f\c\.c\cYf\:iriririririfirkririr 

CC  C CC  CCrffr(V(MC\rr\fVf\(\r-  r-  r-  r-  f>  r-  r-  r~ 

^ r p nr  m r,  ii  ii  ii  ii  ii  n n ii  ii  it  ii  n n n n n n n n ii 

II  II  II  II  II  II  ^ ^ ^ ^ 

ii  - scr-occc^fN.u'  v^f^cca  c a <r 

— i<tcrvDr^ao^r-f-HHr-H^Hr-r'f\jf\i(\j  <\.  f\j  A*  <\j  cm  co  n n 

w -w  w _ w • w ^ • w-^rww^w^-ww 

ifl-KKKH-KH-h-Kfr-KK-HKKKI-KKI-HHKHK^-K 
^KNNNNNNNNNNNNN  NNNNNNNNNKKNN 
>XX>>>>>>>  >>  > > X > > > X 

»-xyxyxxxyxxxxyx  x>xxxxxxxx>x> 


9 1 7 T 02 1 ( 3 ) = 49 , 4 8 » 2 7 » 2 8 » 3 8 
'-,1T08(4)=4  7,4S,2  4»?6,55,5>3»32»34 
^°TP]6(4)=46»36»2b»37»S>4»41*33»4 
r>  1 7 T 0 2 1 (4)  =S0*49,28»29»39 
c lmd  ( *>  * n = * .-in. » 


5.4.4  Input  Data  - Case  2 - Dynamic  - Point  Load 
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Input  Data  - Case  3 - Dynamic  - Uniform  Pressure  (Continued) 


Input  Data  - Case  3 - Dynamic  - Uniform  Pressure  (Concluded) 
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5.4.6  Results 

Case  1 - Static  - Concentrated  Load 

- 0.031335  in  (SAP  IV) 

Maximum  Deflection: 

- 0.036768  in  (Thin  Plate  Theory) 

The  displacement  response  curves  at  the  centroid  for 
the  two  dynamic  forcing  functions  are  shown  in  Figures  64  and  65. 


DYNAMIC  ANALYSIS  OF  SIMPLY  SUPPORTED  THICK  PLATE 
NORMALIZED  DISPLACEMENT  RESPONSE  AT  NODE  8 


TIME  (SECOND) 

COMPONENT  3 

NORMALIZED  FACT  OR : 

.G86E-01 

Figure  64. 

Nodal  Response  at  Plate  Center 
Point  Load 

6.  AIRCRAFT  SHELTER  COMPONENTS 


6.1  Static  and  Dynamic  Analyses  of  24-foot  Radius  Arch 

6.1.1  Introduction 

A 24- foot  radius  arch  was  chosen  to  represent  a segment 
of  a typical  aircraft  shelter  and  was  composed  of  reinforced  concrete 
that  was  equivalent  to  the  reinforced  concrete  covered  doubly  corrugated 
steel  panel  of  a hardened  aircraft  shelter.  The  cross-section  of  the 
shelter  that  was  modeled  is  the  same  as  that  shown  in  problem  4.5. 

6.1.2  Six,  21 -node  thick  shell  elements  were  used  to  model 
the  arch.  A typical  thick  shell  element  with  the  location  of  the  nodes 
and  the  natural  coordinate  system  is  shown  in  Figure  66. 


Figure  66.  21-Node  Element 


The  assignment  of  nodes  and  elements  for  the  model  is  shown  in 
Figures  67  and  68.  The  node  point  coordinates  were  input  in  terms  of 
cylindrical  coordinates  and  converted  to  the  x,  y,  z coordinate  system 
shown  in  the  figure  by  SAP  IV,  which  uses  the  following  relationships: 

x = R(sin0),  y = y,  z = R (cose) 


For  the  arch,  an  equivalent  thickness  was  chosen  to  match  the 
flexural  rigidity  in  the  circumferential  direction  of  a hardened 
aircraft  shelter.  Using  the  relation  I =(b/l3  h3,  with  an  I of  36,189 
in  and  a width  of  24  in  results  in  an 

hof  nnamEi . 26.25  in 

24  in 

6.1.3  Static  Analysis  of  24- foot  Radius  Arch 

For  the  static  analysis  the  arch  was  loaded  with  a 
uniform  pressure  of  36.7  psi.  This  pressure  approximated  the  peak  free- 
field  overpressure  of  the  Mixed  Company  event  at  a range  of  600  feet. 
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XYZT  ( *e.  )=31  »4R 

XYZT( 6? )=3l^.2S,A8 
XYZT  < 6« ) =R14.?5  *68 


Input  Data  - 24-Foot  Radius  Arch  - Static  (Continued) 

XYZT( 751=314.  25*48.  » 18, » 

XYZT(81)=314. 25*48. *0.* 

XYZT (4 1=30] .125*48. » 180. » 

XYZT ( 17)=3ni. 125*48. *144. * 


5MrLLl  ( ? 4 ) »*  *1  * * 
SHTLL2 ( 1 4 ) = *»»*  1 
5 1T08( 11=14*1*3*1 


6.1.5  Results 

As  an  example  of  the  results  for  the  static  load, 
consider  the  vertical  and  lateral  displacements  of  nodes  28  and  54 
which  are  symmetrical  about  the  crown. 


i 


NODE 

VERT.  DISPL. 

LAT.  DISPL. 

28 

-0.0644  in 

+0.0113  in 

54 

-0.0644  in 

-0.0113  in 

6.1.6 

Analysis  of  a 24- 

Foot  Radius  Arch  Under  a Dynamic 

For  the  dynamic  analysis  an  attempt  was  made  to  utilize 
the  data  from  the  Mixed  Company  event  to  construct  a simple  but  realistic 
loading  function.  Figure  69  shows  the  pressure  waves  that  were  chosen 
to  model  the  external  pressure  acting  on  the  arch  at  various  locations. 

As  was  done  in  problems  4.4  and  4.5  the  pressure  was  replaced  by 
equivalent  concentrated  forces  acting  on  22  nodes. 


P(psi) 


Blastward 


Time  (M  sec) 


Figure  69.  Applied  Load  - Arch 
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6.1.7  Dynamic  Results 


1 


The  first  15  eigenvalues  of  the  24-foot  radius  arch 
were  extracted  and  the  forced  dynamic  response  obtained  using  modal 
superposition.  The  first  five  natural  frequencies  are  wi  = 35-55 
Rad/Sec,  to2  = 77.33  Rad/Sec,  u)3  = 109.6  Rad/Sec,  ^4  = 156.9  Rad/Sec 
and  to5  = 229.7  Rad/Sec. 

The  nodal  displacements  corresponding  to  the  first 
natural  frequency  of  vibration  are  for  selected  nodes: 


NODE 

X 

Z 

27 

+0.0096 

+0.0019 

28 

0 

0 

29 

-0.0096 

-0.0019 

AO 

+0.0102 

0 

41 

0 

0 

42 

-0.0102 

0 

53 

+0.0096 

-0.0019 

54 

0 

0 

55 

-0.0096 

+0.0019 

These  nodes  are  located  symmetrically  about  the  crown 
of  the  arch.  The  displacements  listed  above  indicate  that  the  first 
natural  frequency  of  35.55  Rad/Sec  corresponds  to  a torsional  mode  which 
would  not  be  present  in  an  actual  shell. 

6.2  Dynamic  Analysis  of  a Hardened  Aircraft  Shelter  with  No 
End  Walls 


6.2.1  Introduction 

A hardened  aircraft  shelter  modeled  after  the  shelter 
placed  600  feet  from  ground  zero  in  the  Mixed  Company  Event  was  analyzed 
by  the  SAP  IV  computer  program.  The  72-foot  long,  24-foot  internal 
radius  shelter  had  the  same  cross-section  as  the  arch  described  in 
problem  6.1. 

6.2.2  Finite  Element  Model 

The  shelter  was  modeled  with  15  node  thick  shell 
elements  similar  to  those  used  in  problem  6.1.  The  assignment  of  nodes 
and  elements  for  the  model  is  shown  in  Figure  70.  It  should  be  noted 
that  symmetry  was  invoked  so  that  one-half  of  the  shelter  was  modeled 
with  15  elements,  each  of  which  was  12  feet  wide.  Cylindrical  coordi- 
nates were  used  to  input  the  nodal  coordinates. 
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4 ! 


Figure  70.  Element  Assignment  - Aircraft 
Shelter  With  No  Endwalls 


6.2.3  Dynamic  Analysis 

For  the  dynamic  analysis  the  shelter  was  struck  by 
the  same  pressure  waves  as  those  used  for  the  arch  in  problem  6.1. 
The  pressure  wave  was  replaced  by  equivalent  concentrated  forces 
acting  on  36  nodes. 
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6.2.4  Input  Data  - 24-Foot  Radius  by  72-Foot  Hardened  Shelter 
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6.2.5  Results 

The  first  20  eigenvalues  and  eigenvectors  were  de- 
termined and  used  in  the  nodal  superposition.  The  first  five  natural 
frequencies  are  =81.1  Rad/Sec,  w2  = 185.9  Rad/Sec,  u>3  = 240.8  Rad/ 
Sec,  oji*  = 242.2  Rad/Sec,  and  u)5  = 324.7  Rad/Sec. 

For  nodal  points  on  the  24-foot  radius  arch  and  the 
24  x 72  foot  shelter  located  at  similar  locations  near  the  crown 
(Figure  71),  it  can  be  seen  that  the  response  of  the  arch  is  quite 
similar  to  the  response  of  the  shelter,  when  subjected  to  the  pressure 
load  of  problem  6.1. 


Shelter 


Figure  71.  Comparable  Nodes  on  Arch  and  Shelter 
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6.3  Dynamic  Analysis  of  a 24-Foot  Radius  by  72-Foot  Long  Hardened 
Shelter  with  Endwalls 

6.3.1  Introduction 

Fifteen-inch  concrete  endwalls  were  attached  to  the  ends 
of  the  hardened  shelter  described  in  problem  6.1  and  the  shelter  was 
reanalyzed  using  SAP  IV. 

6.3.2  Finite  Element  Model 

An  endwall , composed  of  10-plate  elements,  was  attached 
to  the  shelter  model  described  in  problem  6.1.  Figure  73  shows  the 
assignment  of  nodes  and  elements  for  the  endwall.  The  endwall  was  placed 
at  y=0  to  be  consistent  with  shelter  that  had  been  previously  analyzed. 
Symmetry  was  again  invoiced. 
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Figure  73.  Element  Assignment  Shelter  Endwall 


6.3.3  Input  Data  - Hardened  Shelter  with  Endwall 
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Input  Data  - Hardened  Shelter  with  Endwal'l  (Concluded) 
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The  shelter  with  endwalls  was  struck  by  the  same  pressure  wave  as 
described  in  problem  4.4  and  6.2.  No  pressure  has  been  applied  to  the 
endwall  so  that  a comparison  might  be  made. 

The  first  five  natural  frequencies  are  n>i  = 126.2  Rad/Sec,  ^2= 

203.0  Rad/Sec,  oj3  = 236.5  Rad/Sec,  = 330.4  Rad/Sec,  and  w5  = 349.9 
Rad/Sec.  The  slightly  higher  natural  frequencies  confirm  that  the  shelter 
with  endwalls  is  indeed  slightly  stiffer  than  the  open-ended  shelter. 

The  eigenvalues  and  eigenvectors  were  used  in  the  subsequent  mode 
superposition  analysis. 

The  maximum  displacements  due  to  the  pressure  waves  at  a point 
24  feet  from  the  end  wall  (y  = 144  in)  for  the  two  shelters  that  were 
analyzed  are: 


24  x 72 

ft  Shelter 

Node 

Di spl acement 

Max  Value 

Time  at  Max 

50 

AX 

0.3875 

0.0280 

50 

AZ 

0.6308 

0.0260 

51 

AX 

0.1534 

0.0140 

51 

AZ 

0.5733 

0.0260 

52 

AX 

0.3511 

0.0640 

52 

AZ 

0.6106 

0.0280 

24  x 72 

ft  Shelter 

With  Endwall 

at  y = 432 

in 

Node 

Displacement 

Max  Value 

Time  at  Max 

50 

AX 

0.2323 

0.0240 

50 

AZ 

0.3110 

0.0220 

51 

AX 

0.1279 

0.0120 

51 

tz 

0.2728 

0.0220 

52 

AX 

0.2165 

0.0160 

52 

AZ 

0.2974 

0.0200 

It  can  be  seen 

that  when  the 

shelter  is 

struck  by  a dynami 

the  endwall  has  a significant  stiffening  effect. 

Figures  74  and  75  show  two  views  of  the  element  layout  for  the 
hardened  shelter  with  endwalls.  These  plots  were  produced  using  the 
mesh  plot  routine  of  Appendix  C. 
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6.4  Application  of  SAP  IV  to  a Prowed  Door 


6.4.1  Loading 


A prow-shaped  rolling  TAB-VEE  door  is  subjected  to  the 

following  loads: 

Case  1 - Static  uniform  pressure 

Case  2 - Dynamic  uniform  pressure  with  the  time  history  shown  in 


Pressure  (Psi) 


0.00  0.0182 


Loading  Function 


Figure  76.  Prow-shaped  Rolling  TAB-VEE  Door 


I 


6.4.2  Finite  Element  Model 

The  prow-shaped  TAB-VEE  door  consists  of  three  parts: 
collar  bulkhead,  armor  door,  and  floor  frame.  The  finite  element  model  for 
this  structure,  shown  in  the  following  diagrams,  includes  three  groups  of 
elements:  thick  shell,  thin  plate,  and  beam  elements.  The  thick  shell 
and  the  quadrilateral  thin  plate  elements  are  used  to  represent  the  bulk- 
head (Figure  77)  and  the  armor  door  (Figure  78),  respectively.  The  beam 
elements  are  used  to  model  the  reinforcing  ribs  and  the  floor  frame  (Figure 
79).  Also,  because  of  symmetry  only  half  of  the  door  has  to  be  modeled  with 
symmetry  conditions  imposed  to  the  center  plane  of  the  door. 

The  dynamic  calculation  was  performed  in  two  parts.  The 
first  set  of  data  for  the  dynamic  calculation  generates  the  eigenvalues  of 
the  armored  door  and  writes  a restart  tape.  The  second  set  of  data  causes 
SAP  IV  to  read  the  restart  tape  and  generate  the  model  solution  using  the 
calculated  eigenvalues. 


Numeral  - Nodal  number 

Q - Thick  shell  element 
number 


37^  64// 


71  // 


Figure  77.  Bulkhead 
143 


5^56 

Via  \ 

45  44  43U 

1 Q \ 

£147 1 48  x 

46  52  42  I 


Numeral 

A 


Nodal  number 
Beam  element  number 


Figure  79.  Floor  Frame 
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6.4.3  Input  Data  - Static  Uniform  Pressure  on  One  S^de  of  Door 
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Input  Data  - Static  Uniform  Pressure  on  One  Side  of  Door  (Continued) 
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Input  Data  - Eigenvalue  Calculation  and  Restart  Creation  (Continued) 
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I PLATE  < 3 ) = ?6»1 6 *1 7*27 


Input  Data  - Eigenvalue  Calculation  and  Restart  Creation  (Concluded) 
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6.4.5  Input  Data  - Dynamic  Response  Using  Restart  Tape 


HEAD£R=*REbT4RT  RtSPONbfc  H I b TOR  Y ANALYSIS  UF  AN  ARRiOK  LJOUk* 
NUMNP=72 »NElTYP  = 3 , NF  = 24  *NDYN  = -2 
NFN  = 1 *NT  = 500*NOT=4,r>T=.0001 
MP( 15*1 ) =1 , ,-10889. 
mp( 15  *3 ) = 1 ,,-14008.2 
MP( 1 7,1 )=1 , ,-21778. 

MP ( 1 7 * 3 ) = 1 , ,-28015.9 

MP{  1 8 ,1  ) =1  ,,-21 778 . 

MP ( 1 8 , 8 ) =1  , ,-280  16 . 5 
NP< 21  *1  )=1  , ,-21778. 

NP( 21 , 3 ) = 1 , ,-28016.5 
NP<  23*11=1  , ,-21778. 

\'P(  23*3  1 =1  ,,-28016.8 
MD( 2 6*11  = 1 , *-10889. 

MP ( 29 ,9 1 =1  ,,-20379.8 
NP( 28,1  )=t  , ,-l  9802 . 5 
M l*  ( 2 6 * 2 1 = 1 ,,  — 14581.5 
MP( 26*91=1 , ,-25497.9 
*'P(  27»i  )=i  ,,-9870  0.1 
NP( 27,? 1=1 , ,-24170.9 
NP( ?7, 9 ) =1 ,, -4 756ft. 6 
MP{ 2 8*1  1=1  , ,-31160.5 
NP< 28*2  1 = 1 , ,-16109.6 
NP( 28*31=1 , ,-40110.0 
N P ( 2 9 , l )=?  , ,-28665  . 

M°( ?9,7)=1 , ,-11072. 6 
NP( 29  *9 ) =1  ,,-9ft88o,8 

MP ( 90*1 1=1 , * -2 9?99 . 7 
NP(  90»? ) = 1 ,,-12110.1 
MP ( 30*3  1 = 1 ,,-37711  . 
m p ( 31*11=1, ,-l  5755.4 
MP( 31  *21  = 1 ,,-7R39.4 
NP( 31 ,3 ) = l , ,-74739.5 
NP( 32*11=1, , -17826. y 
N P ( 32  *2  1 = 1 * ,-28703.0 
K|P ( 72,3  1 = 1 , ,-22978.7 
NP( 93,11=1 , ,-30023.8 
NP(93,2)=1 ,,-48341.8 
MP( 33,3  )=1  , ,-3R70l  . 

N p < 34,1  1 =i  , ,-] 3 7ft 6 . ? 

*)P(  n,?)=1  , ,-30213.7 
MP(  34,3  1 = 1 ,,’-2  4 188.1 
MP ( 36  * 1 ) = 1 , ,- 1 376  4 . 

MP(  35*21=1  * ,-22149./. 

MR( 56*3 )=1 , ,-17728.9 
NP  ( 3ft  , 1 ) = 1 ♦ 1 01 7ft . 2 
NP( 36*2  1=1  * ,-16384.7 
NP( 36*31=1 **-13117. 


Ntt#. 


Input  Data  - Dynamic  Response  Using  Restart  Tape  (Concluded) 


mp  ( 7 7 

» 1 1 =i  * 

,-891 7.5 

NiD  ( ■>  7 

♦ 7 ) = 1 ♦ 

.-14351 .5 

M P ( 07 

» 8 ) = 1 ♦ 

♦-23129.° 

MD  ( o p 

» 7 1=1  * 

♦ -15012. 1 

»'D  ( 3 8 

>7 ) =1  » 

♦-26 170.° 

M O ( OP 

♦ 7 ) =1  , 

,-40980.7 

MD  ( "3  0 

* l )=i  , 

♦-9787.6 

MP  ( 30 

♦ 7 ) =i  ♦ 

♦-1 6106.6 

MP  ( 3 g 

♦ ° 1 = 1 ♦ 

♦-4 9139. 

► r>  ( 4 n 

♦ 1 ) = 1 ♦ 

♦ -•6877. 

MD(  40 

♦ 7 1 = 1 * 

♦-1 1072.6 

MP  ( 40 

♦ 7 ) =1  , 

♦043077.4 

MD  ( Lx  1 

♦ 1 > = 1 ♦ 

.-7521.3 

MD  ( 41 

*71=1  * 

♦-1711°. 1 

MP  ( 4] 

*7  ) = 1 ♦ 

-33541.5 

MP  ( 4? 

>3 ) =1  * 

♦-6769.4 

M D ( AO 

• o \ = , , 

,-14566  .') 

MD  ( P 0 

*7  ) = 1 , 

♦-1 6771.7 

MP(  58 

*7  ) =1  , 

♦ -26884 .4 

'io(51 

»7 )=1 » 

♦-73815.9 

NP(  62 

*7 ) =1  , 

♦-12774.7 

•NIP  ( 67 

>7  ♦ 

♦-7S079.3 

MP  ( 70 

’7  ) =1  , 

♦ — 41483.7 

Ml.P  ( 1 

1=40  . 

T ( 1 .1 

> =0. * T ( 1 *2 1 = .001  * f 

c ( 1 ♦ ! 

l =0. ♦FI 1 >7 ) = ] . *F(  1 

fc'K’K  = ?*lSP=7 
T Comp ( ?6^i 

♦3? .77. 78) =1 

risr  S = 

2 * 15P8 

= 7 

tS(  ?♦ 

16  ) =1  ♦ 

? *3*4-*5  **> 

T S ( ?* 

IP  ) = 1 ♦ 

7 ♦ 7 ♦ 4 ♦ 5 ♦ 6 

T S ( ? ♦ 

21 ) =i  » 

7^7^4^5^6 

T S ( 4 » 

1 ) =1  ♦ 7 

♦ 7 ♦ 4 ♦ 6 * 6 

T S ( 6 * 

2 ) = 1 ♦ ? 

♦3^4^5»fe 

1 S ( 6 » 

3 1=1*2 

♦ 3 ♦ 4 ♦ 5 ♦ 6 

7 ) = . Cl 82*1  ( l = 1 . 

=o. »f { 1 ) =o. 


i 


I. 


7.  SUMMARY 

These  problems  have  served  to  provide  some  guidance  with  regard  to 
the  accuracy  that  may  be  expected  with  SAP  IV.  For  example,  the  static 
response  of  beams  and  plates  agrees  quite  well  with  classical  theory  with 
very  few  beam  and  plate  elements  respectively.  Similarly,  the  use  of  a 
small  number  of  thick  shell  elements  with  very  large  aspect  ratios  appeared 
to  be  adequate  for  thick  plates,  arches  and  shells. 

For  those  cases  where  other  analytical  predictions  or  experimental 
data  were  available,  the  results  of  SAP  IV  were  reasonably  close  when 
account  is  taken  of  the  rather  large  elements,  different  boundary  con- 
ditions and  linear  assumptions  that  were  used  throughout.  There  is  no 
question  that  SAP  IV  can  be  extremely  useful  in  assisting  an  engineer  to 
optimize  an  aircraft  shelter  or  shelter  system  design  with  respect  to 
protection  or  cost.  This  latter  conclusion  is  based  on  the  analysis  of 
the  shelter  and  the  prowed  door  that  represent  fairly  complicated  but 
realistic  components  of  a shelter  complex. 

These  problems  have  also  illustrated  the  value  of  the  free-format 
input  program  and  the  graphics  package.  For  each  problem,  each  input 
card  contains  nontrivial  data  that  is  easily  constructed  and  understood. 

The  flexibility  of  the  plot  package  is  indicated  with  the  various  views 
of  the  finite  element  models  as  well  as  with  the  time  histories  of  several 
parameters  that  are  significant  to  the  engineer. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  initial  decision  of  choosing  SAP  IV  appears  to  have  been  a good 
one.  The  problems  associated  with  installing  and  using  the  program  were 
much  fewer  than  might  be  expected  for  a large  code.  Several  problems  with 
classical  solutions  were  run  to  gain  familiarity  with  the  code  and  the 
numerical  results  were  accurate.  Furthermore,  the  solution  routines, 
especially  the  modal  technique  for  dynamic  problems,  were  efficient  and 
easy  to  use. 

Although  the  basic  code  was  quite  satisfactory,  it  was  evident  that 
very  little  effort  had  been  made  to  make  the  input  and  output  data  conven- 
ient for  the  user.  The  input  data  followed  a rather  rigid  format  with  the 
requirement  of  blank  spaces  and  blank  cards  very  common  and  somewhat  con- 
fusing at  best.  There  was  no  method  to  check  the  input  geometric  data 
except  by  repeated  scrutiny,  which  is  not  a very  viable  method  especially 
for  large  problems.  The  output,  in  the  form  of  printed  data  and  plots  from 
a printer,  was  cumbersome  and  rather  unsuitable  for  the  analysis  of  compli- 
cated structures.  Furthermore,  it  was  evident  that  the  code  had  been  con- 
structed primarily  with  metal  structures  in  mind,  since  there  was  no 
indication  of  property  identification  for  concrete  or  reinforced  concrete 
structures. 

Since  the  primary  objective  of  this  study  was  a code  that  could  be  used 
in  the  dynamic  analysis  of  concrete  and  reinforced  concrete  structures,  it 
was  necessary  to  derive  a method  for  computing  equivalent  linear  bending 
stiffnesses.  For  any  realistic  dynamic  environment  a given  structural 
member  will  experience  both  positive  and  negative  moments.  Since  many 
members  used  in  aircraft  shelter  systems  are  not  symmetrical  with  respect 
to  positive  and  negative  bending,  it  was  necessary  to  extend  an  existing 
method  to  this  case.  The  equivalent  bending  stiffness  that  was  proposed 
was  an  average  of  the  stiffnesses  for  the  cracked  and  uncracked  sections 
for  both  positive  and  negative  bending.  The  results  appear  to  be  reasonable 
based  on  a limited  amount  of  other  analytical  data. 

The  free- format  input  program  that  was  developed  simplifies  the  input 
of  data  for  the  average  user.  This  program  takes  data  in  a simplified  or 
free-format  form  and  converts  it  to  the  format  required  by  the  SAP  IV  code. 
Thus,  the  user  has  the  option  of  using  the  free-format  form  or  of  bypassing 
this  program  and  submitting  data  in  the  SAP  IV  format.  With  the  free- 
format  input  program,  cards  can  be  input  in  any  order  and  only  nonzero  data 
is  required.  When  certain  data,  such  as  boundary  conditions  are  repetitive, 
or  nodes  are  equally  spaced,  the  program  generates  much  of  the  data 
for  the  user.  The  result  is  that  the  free-format  input  consists  of  the 
minimal  amount  of  data  that  is  required  for  defining  a problem.  Not  only 
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is  this  extremely  convenient  to  the  user,  but  the  chance  for  error,  such 
as  a wrong  number  or  a correct  number  in  a wrong  column,  is  significantly 
reduced. 

Even  with  the  free-format  input  there  is  still  the  possibility  of 
errors  in  node  or  element  definition  that  are  extremely  difficult  to 
detect.  The  mesh  plot  program  was  devised  to  assist  in  this  important 
aspect  of  data  checking.  Furthermore,  a graphical  display  of  elements 
can  easily  show  elements  that  have  a bad  aspect  ratio  and  can  give  a 
general  indication  concerning  the  degree  of  refinement  that  is  almost 
impossible  with  printed  data. 

The  mesh  plot  program  has  several  options  that  were  incorporated 
to  assist  the  user.  The  viewing  point  is  specified  and  can  be  changed 
from  plot  to  plot.  Some  of  the  options  that  are  available  include  the 
f ol 1 owi ng : 

1.  A plot  of  all  nodes  and  elements  with  labels 

2.  An  unlabeled  plot  of  all  nodes  and  elements 

3.  Labeled  or  unlabeled  plots  of  just  the  beam,  plate  or  thick 
shell  elements 

4.  Labeled  or  unlabeled  plots  of  specified  subsets  of  beam,  plate 
or  thick  shell  elements. 

The  regular  output  of  SAP  IV  for  dynamic  problems  consists  of 
specified  numerical  data  at  each  time  step  and  printer  plots  of  the  same 
information.  Neither  of  these  forms  is  generally  adequate  for  efficiently 
studying  the  response  of  a structure  or  for  incorporating  such  data  in 
a report.  Thus,  a time  history  program  was  constructed  which  maintained 
essentially  the  existing  control  parameters  as  those  used  for  the  printer 
plots,  but  which  allowed  the  use  of  plotters  that  are  usually  available  at 
any  large  computer  installation. 

A review  of  the  output  parameters  for  a typical  dynamic  analysis  re- 
vealed that  time  history  plots  could  be  conveniently  combined  in  groups 
of  three.  According  to  the  particular  element  and  node  model,  these  groups 
were: 

1.  Displacements  at  a node 

2.  Rotations  at  a node 

3.  Bending  moments  and  torque  for  a beam  element 

4.  Shear  resultants  and  axial  force  for  a beam  element 
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5.  Membrane  resultants  for  a thin  plate  element 

6.  Moment  resultants  for  a thin  plate  element 

7.  Normal  stress  components  for  a thick  shell  element 

8.  Shear  stress  components  for  a thick  shell  element 

Within  any  such  group,  one,  two  or  all  three  histories  can  be  requested. 
The  largest  absolute  value  of  the  requested  functions  over  the  specified 
time  period  is  used  as  a normalizing  factor  and  this  factor  together 
with  identifying  symbols  are  printed  with  each  graph. 

Although  these  additions  to  SAP  IV  have  greatly  enhanced  its 
suitability  for  the  analysis  of  aircraft  shelter  structures,  there  still 
remain  several  items  which,  if  addressed  successfully,  could  improve  the 
usefulness  of  SAP  IV  even  further.  These  are: 
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force  vector  is  computed  and  stored  in  core  for  all  time  steps  prior  to 
the  integration  of  the  modal  equations.  Consequently,  the  core  require- 
ments of  a problem  can  vary  quite  significantly  depending  upon  the  number 
of  time  steps  requested  for  solution  and  can,  on  occasion,  cause  a 
problem  to  exceed  storage  limits.  It  appears  that  this  difficulty  can 
be  corrected  by  an  alteration  to  SAP  IV  in  which  the  generalized  forcing 
function  vectors  are  computed  and  stored  outside  of  core  to  be  accessed 
when  needed.  The  required  vectors  can  then  be  brought  into  core  either 
sequentially  or  in  blocks  as  the  nodal  integration  proceeds. 

2.  Inclusion  of  a capability  of  calculating  and  plotting  of  principal 
stresses,  strains  and  directions.  This  feature  would  be  very  helpful  in 
determining  the  possibility  of  failure. 

3.  Additional  graphical  display  showing  both  deformed  and  undeformed 
states  of  the  structure. 

4.  Inclusion  of  advanced  band  with  minimization  techniques  which 
would  decrease  the  computer  time  necessary  for  a particular  problem. 

5.  Resolving  the  1 initiation  of  a maximum  of  eight  displacement 
components  for  output  plotting. 


/ 
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APPENDIX  A 

FREE  FORMAT  INPUT  PROGRAM 


The  Free  Format  Input  Program  (FFIP)  produces  an  input  data  file 
acceptable  to  the  SAP  IV  program  through  extremely  flexible  input  in- 
structions and  is  available  for  three  element  types  - beam,  thin  plate, 
and  thick  shell  (21  nodes). 

The  two  basic  parts  in  FFIP  are  key  words  and  entries.  These  two 
parts  are  always  linked  with  an  equal  sign.  A key  word  must  appear  be- 
fore an  equal  sign,  and  the  entry  must  appear  to  the  right  of  the  equal 
sign.  Key  words  may  be  nonsubscripted,  single  subscripted,  or  multiple 
subscripted.  For  each  single  subscripted  and  multiple  subscripted  key 
word,  a set  of  parentheses  must  be  included  before  an  equal  sign.  The 
entries  may  be  single  or  multiple  and  may  be  in  floating  point,  fixed 
point,  or  alpha-numeric.  This  information  is  specifically  stated  for 
each  key  word.  Variables,  values  or  strings  of  characters  may  be  placed 
anywhere  on  a card  as  long  as  three  elements  are  intact  in  the  following 
order:  key  word,  equal  sign,  and  entry.  More  than  one  key  word  can  be 
placed  on  one  card,  again,  as  long  as  a set  of  three  elements  as  described 
earlier  remains  intact  and  is  separated  from  another  set  by  a comma  or 
a blank.  Furthermore,  when  two  or  more  key  words  are  punched  on  one 
card,  the  number  of  entries  required  for  a previous  key  word  must  be 
completely  satisfied.  All  key  words  listed  in  this  appendix  are  recog- 
nized by  FFIP,  and  any  misspelled  key  words  will  not  be  recognized. 
Diagnostic  error  messages  will  be  printed  for  misspelled  key  words,  and 
the  job  will  be  aborted. 

Further  requirements  for  input  include: 

- There  must  not  be  any  blanks  within  an  exponential  field. 

- An  integer  or  floating  point  number  may  be  terminated  by  a comma. 

Consecutive  commas  may  be  used  to  indicate  missing  fields,  but 
care  should  be  taken: 

7 . . 05  means  the  second  number  is  missing  (default=0)  but 

70., |tf5  means  the  second  and  third  numbers  are  missing. 

- A * + 1 in  column  80  causes  the  next  line  to  be  read  and  the  free 

format  scan  continues.  (Usage  of  this  feature  should  be  spared  to 
preserve  the  flexibility  of  the  data  deck.) 
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FFIP  KEY  WORDS  AND  ENTRIES 


Explanations: 

Keyword  and  type 
of  allocation 


No.  of  entry 

Type  of  entry 
Usage 


Ref. 


key  word  and  type  of  allocation  which  may  be 
NS,  SS  or  MS  where: 

NS  - nonscripted 

SS  - single  subscripted 

MS  - multiple  subscripted 

number  of  values  to  be  entered  to  the  right  of 
the  equal  sign 

self-explanatory 

the  first  line  is  reserved  for  description  or 
usage.  Another  line  derailed  by  two  dashes 
gives  an  example  of  how  input  should  appear 
on  card. 

more  specific  information  on  usage  is  referenced 
to  the  SAP  IV  manual  by  section  and  page  number. 


U5HDE 


REF. 


ICYWORD  HND  NO.  TTFE 

TYPE  OF  OF  CF 

ALLOCATION  ENTRY  ENTRY 


ALPHA 

NS 

1 

REAL 

.DAMP INC  FACTOR  FOR  A FORCED  DYNAMIC 
.RESPONSE  ANALYSIS 
. — ALPHA~ENTRY 

RT 

S3 

1 

PEAL 

.ARRIVAL  TIME 

AT(El)=E2 

BEAM  SS  12  INTEGER 


BELFX  NS  >4  REPL 


BELFY  NS 


4 REAL 


BELFZ  N5 


4 fEFL 


WHERE 

ElrflRRIVRL  TINE  NUMBER  (INTEGER! 

• LE.  NOT 
E2=ENTRT 
BEAM  DATA 

— BE AM (El I =E2,E3. E4.E5. E6, E7. E3.  E9.E10.E11 . 
E12.E13 
WHERE 

El -BERM  ELEFENT  NUMBER  .LE.  FBERM 
E2=N0OE  NUMEER  I 
E3=NBDE  NUMBER  J 
E4=N0DE  NUMEER  K 
ES=MATERIPL  PROPERTY  NUMBER 
ESrELEHENT  PROPERTY  NUMBER 
E7;FIXEQ  END  FORCE  10  FOR  ELEMENT  LOAD 
CR3E  R 

E8=FIXED  END  FORCE  ID  FOR  ELEMENT  LORD 
CASE  8 

E9:FIXED  END  FORCE  ID  FOR  ELEMENT  LOAD 
CASE  c 

E10=FIXED  ENT  FORCE  ID  FOR  ELEFCNT  LORD 
CASE  D 

EllrENO  RELEASE  CODE  RT  NODE  I 
E12=END  RELEASE  CODE  AT  NODE  J 
E13rDPTI0NAL  PARAMETER  K FOR  AUTOMATIC 
GENEFFTIDN 

ELEMENT  LOAD  FACTOR  - MULTIPLIER  OF  GRAVITY 
LORD  IN  THE  +X  DIRECTION 

— BELFXrEl,  E2.E3.  E4 
WHERE 

EUELEMENT  LOAD  CASE  A 
E2;ELEMENT  LOAD  CASE  R 
E3=ELEMENT  LORD  CASE  C 
E4:ELEHENT  LOAD  CASE  D 
ELEMENT  LORD  FACTOR  - MULTIPLIER  OF  GRAVITY 
LOAD  IN  THE  +Y  DIRECTION 
— BELFY =E  1 , E2 . E3.  EU 
WHERE 

ElrELEMENT  LOAD  CASE  A 
E2=ELEMENT  LORD  CASE  B 
E3=ELEHENT  LOAD  CASE  C 
EU=ELEMENT  LOAD  CASE  D 
ELEMENT  LORD  FACTOR  - MULTIPLIER  OF  GRRVITT 
LOAD  IN  THE  +Z  DIRECTION 
--8ELFZrE2.E2.E3. E4 
WHERE 

E1=ELEWENT  LOAD  CASE  A 
E2=ELEHENT  U3A0  CASE  B 
E3rELEMENT  LOAD  CASE  C 
E4=ELEMENT  LORD  CASE  0 


VII. 8 
VII. 12 

IV. 2.3 


IV. 2. 2 

IV. 2. 2 

IV. 2. 2 
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KEYWORD  AND  NO.  TYPE 

TYPE  QF  OF  CF 

ALLOCATION  ENTRY  ENTRY 


BEPC  SS  6 FEAL 


BETH 

N5 

1 

FEFL 

. DAMPING  FHCTOfi  FOR  ft  FORCED  DTNRMIC 
.RESPONSE  ANALYSIS 
. — BETA=ENTRY 

.VII. 8 

BFEFI 

SS 

6 

REAL 

.BEAM  FIXED  END  FORCES 
. —BFEFI  (El)  =£2.  E3.E4,  E5.E6.  E7 
. WHERE 

.IV. 2. 2 
* 

■ 

BFEFJ  SS  6 PEAL 


BEAM  ELEMENT  PFiCPERTY 

-BEPC  ( El)  =E2 . E3,  EH . E5 , E6 , E? 

WHERE 

E ^GEOMETRIC  PROPERTY  NUMBER  (INTEGER) 
.LE.  Bf£PC 
E2=AXIRL  FAEA 

E3:SHEAR  FRER  ASSOCIATED  WITH  SHEAR 
FORCES  IN  LOCAL  2-DIRECTION 
E4=3HERR  FfiEfl  R330CIRTED  NITh  SHEAR 
FORCES  IN  LOCAL  3-DIRECTION 
E5=TDRSI0NAL  INERT  I ft 
ESrFLEXURFL  INERTIA  ABOUT  LOCAL  2-AXIS 
E7=FLEXURAL  INERTIA  ABOUT  LOCAL  3-AXIS 


E INFIXED- END  FORCE  NUMBER  (INTEGER) 

.LE.  BhfEFS 

E2=FIXED-END  FORCE  IN  LOCAL  1-DIPECTI0N 

RT  NODE  I 

E3:FIXED-END  FORCE  IN  LOCAL  2-DIPECTION 

AT  NODE  I 

EH^FI  XED-END  FORCE  IN  LOCRL  3-DIf=ECTlON 


AT  NODE  I 

E5=PIXED  END  MONENT  ABOUT  LOCAL 

1 - DIRECT ION  AT  NODE  I 

E6=FIXED  END  MOMENT  ABOUT  LOCAL 

2- DIPECTION  RT  NODE  ( 

E7:FIXED  END  MOMENT  ABOUT  LOCAL 

3- DIRECTION  AT  NODE  I 
BEAM  FIXED  END  FORCES 
—BFEFJ  (El)  =£2.  E3.EN,  E5.E6.  E7 

WHERE 

ElnFIXEO-ENO  FORCE  NUMBER  (INTEGER) 

.LE.  BAFEF3 

E2;FIXED-END  FORCE  IN  LOCAL  1-DIRECTION 

AT  NODE  J 

E3=FIXED-END  FORCE  IN  LOCAL  2-DIRECTION 

AT  NODE  J 

E4=FIXED-END  FORCE  IN  LOCRL  3-DIFECTION 

AT  NODE  J 

E5:FIXED-END  MONENT  ABOUT  LOCAL 

1 - DIRECT ION  AT  NODE  J 

E6=FIXED-END  MOMENT  ABOUT  LOCAL 

2- 01 FECT  ION  RT  NODE  J 

E7=FIXED-END  MOMENT  ABOUT  LOCAL 

3- DIFECTION  AT  NODE  J 


IV.2.1 


IV. 2. 3 


I 


KETNQRD  AND  NO.  TYPE 

TYPE  OF  OF  OF 

ALLOCATION  ENTRY  ENTRY 


BMPC 


SS 


REAL 


BNEPC 

NS 

1 

INTEGER. NUMBER  OF  BEAM  ELEMENT  PROPERTY  SETS 
. — BNEFC=ENTRY 

.IV. 2.1 

ENFEF5 

NS 

1 

INTEGER. NUHBER  OF  FIXED  END  FORCE  SETS 
. — BNFEFS=ENTRY 

IlV.2.1 

BNNFC 

NS 

I 

INTEGER. NUNBER  OF  BERM  MATER I RL  PROPERTY  SETS 
. — BNMPC=ENTRY 

1 IV. 2.1 

CLNO 

MS 

6 

REAL  .CONCENTRATED  LOAO/MASS  DflTfl 

. — CLMD  ( El . E2)  :£3.  E4 , E5 . E6 , E7 . E8 
. WHERE 

!v,  i 

USAGE 


BEAM  MATERIAL  PROPEATY 
— BMPCIE1)=E2.E3.E4,E5 
WHERE 

E IMMATERIAL  ID  NUMBER  (INTEGER) 
.LE.  BKMPC 
E2=YBUNGTS  MODULUS 
E3rP0ISSONfS  RATIO 
EU=MASS  DENSITY 
ESrHEIGTH  DEN3TIT 


ElrNDDAL  FOINT  NUHBER  (INTEGER) 

.LE.  NLMKP 

E2M3TRUCTURE  LOAD  CASE  NUMBER 
(INTEGEFi)  .LE.  LL 
GE.l.  STATIC  RNRLYSIS 
EQ.O.  DYNRHIC  ANALYSIS 

E3=X-DI RECT ION  FORCE  OR  TRANSLATIONAL 
MR35  CCEFFICIENT 

E4=Y-0IRECTI0N  FORCE  OR  TRANSLATIONAL 

MRS5  CCEFFICIENT 

E5:Z-DIRECTI0N  FORCE  OR  TRANLS AT IONAL 

MASS  CCEFFICIENT 

E5:X-RXI5  MOMENT  DR  ROTATIONFL  IfERTIH 

E7=Y-AXIS  MOMENT  OR  ROTATIONAL  IKERTIA 

EOzZ-flYIS  MCMENT  OR  ROTATIONFL  IlERTIfi 


O0FQ 

NS 

1 

REAL 

.CUT-OFF  FREQUENCY 
. — COFQrENTRY 

.VII. 3 

J 

CT 

SS 

1 

ALPHA 

. SYMBOL  DESCRIBING  COORDINATE  SYSTEM  FOR 

illl.l 

THIS  NODE 
EQ.  , (BLANK) 

EQ.C, 

— CTCEl i =E2 

WHERE 

E1=NOOAL  POINT 
.LE.  NUMfE 
E2:ENTRY 


CARTESIAN  IX.Y.Z) 

CYLINDRICAL  (R,Y.  THETA) 


NUHBER  (INTEGER) 


DAMP 

NS 

1 

FEFL 

. DAMPING  FRCTOR  (NDYN=2) 

. — OAMPrENTRY 

.VII.B 

DT 

NS 

1 

PEAL 

.SOLUTION  TIME  STEP 
. — OT^ENTRT 

I Vt I , 8 

ELH 

SS 

4 

FEFL 

.ELEMENT  LOAD  MULTIPLIERS 
.—  ELM  (El)  =E2.  ES.E4.E5 

!vi.i 

WHERE 

EI=LDAD  CASE  NUHBER  (INTEGER)  LE.  LL 
E2=MULTIPLIER  FOR  ELEMENT  LOAD  CASE  A 
E3rMULTIPLIER  FOR  ELEMENT  LOAD  CASE  B 
EmMULTIFLIER  FOR  ELEMENT  LDFD  CASE  C 
E5:MULTIPLIER  FOR  ELEMENT  LOAD  CASE  0 


REF. 


IV. 2.1 
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KEYWORD  RND  NO.  TYPE  . U5RGE 

TYPE  OF  0F  I . 

ALLOCATION  ENTRY  ENTRY  . 


PEAL  .TINE  FUNCTION  DEFINITION  FOR  RESPONSE  .VII. Id 

.HISTORY  ANALYSIS 
F£El.E2)iE3 
. WHERE 

ElrTIME  FUNCTION  NUMBER  (INTEGER) 

.LE.  NFN 

E2rDEFINITIGN  POINT  NUMBER  [INTEGER) 

.LE.  NLP (NFN) 

E3=FUNCTION  VALUE 

PEFL  .FACTOR  FDR  X-DIRECTION  INPUT  FOR  .VI I. 23 

. RESPONSE  SPECTRUM  ANALYSIS 
. — fX-ENTRT 

REAL  .FACTOR  FOR  Y-DIFECT ION  INPUT  FOR  . V [ I . 23 

.RESPONSE  SPECTRUM  RNHLT5I5 

. —FYrENTRY 

PEFL  .FACTOR  FDR  Z-DIPECTION  INPUT  FOR  .VI I. 23 

.RESPONSE  SPECTRUM  ANALYSIS 

. — FZ=ENTRY 

ALFPifl  . HERDING  INFORMATION  .1.1 

. — HEADER:*52CHARACTERS  INCLUDING  BLANK* 

FLFHfl  .LABEL  INFORMATION  FOR  RESPONSE  MISTORT  .VCI.13 
.ANALYSIS 
.— HEO(El)=*E2* 

. WHERE 

E1=FUNCTIQNAL  DEFINNITION  PD1NI  NUNBtH  . 

. (INYEDER)  .LE.NFN 

E2:5B  CHARACTERS  INCLUDING  BLANK 

ALPHA  .HERDING  INFORMATION  USED  TO  LRBEL  THE  .VC 1 . 

.SPECTRUM  TRBLE 

.— HEDSP=*5B  CHARACTERS  INCLUDING  BLANK ■ 

INYEGEfl. DISPLACEMENT  CUTFUT  ORTH  .VII. IB 

. --3C0MP  (El)  =E2,  E3.EU,  E5.E6,  E7 

. WHERE 

ElrNDORL  POINT  NUMBER  ( INTEGER) 

. .LE.  NUMNP 

E2=0ISPLACEPENT  COMPONENT,  REQUEST  1 
£2:01  SPLPCENENT  COMPONENT.  REQLEST  2 
. E4:DI 3PLACENENT  COMPONENT,  REQUEST  S 
E5:DISPLACENENT  COMPONENT.  REQUEST  U 
E6=DI SPLACENENT  COMPONENT,  REQUEST  5 
E7:0ISPLACENENT  COMPONENT,  REQUEST  6 

INTEGER. FLAG  FOR  PRINTING  INTERMEDIATE  MATRICES  .VI I. 3 
. — IFPRrENTRT 

INTEGER. FLAG  FOR  PERFORMING  THE  STURM  SEQUENCE  CHECK. VI 1. 3 

. — IFSS=ENTRY 

INTEGER. STANDARD  INTEGRATION  ORDER  FOR  TPE  NATURAL  .IV. 8. 5 
. IR.SI  DIRECTIONS 
. — )NTR5=ENTRT 

INTEGER. STANDARD  INTEGRATION  ORDER  FOR  TPE  NATURAL  .IV. 8. 5 
. IT)-OIRECTION 
. — )NTT:ENTRY 


PEADER  NS  6 


ICOMP  MS  6 


I 


I I 


MEYWORD  RND  ND.  TTFE 

TYPE  OF  OF  CF 

ALLOCATION  ENTRY  ENTRY 


IPLATE  SS 


7 INTEGER 


12  INTEGER 


KEGB  NS 


KKK  NS 


1 INTEGER 


1 INTEGER 


PLATE  DATA 

--1PLATE ( El ) =E2 . E3 , E4 . E5 . E6 , E7 , E8 

WHERE 

E1=PLATE  ELEMENT  NUNBER  .LE.  NPLRTE 

E2=NDDE  I 

E3=N0DE  J 

E4=NDDE  K 

E5=NODE  L 

EStNODE  0 

E7:MATERIPL  IDENTIFICATION  NUMBER 
EB=ELEHENT  DfiTfl  GENERATOR 
ELEMENT  STRESS  COMPONENT  OUTPUT 
— 1SIEUE2)  =E3,EU.E5,E8.E7,E8,  EB.ElO.Ell ,E12 
E13.E11 
WHERE 

EUELEMENT  TYPE 
E2;ELEHENT  NUMBER 
.LE.  NEE  AM.  IF  E1=2 
.LE.  NPLRTE,  IF  El  =6 
.LE.  NS0L21.  IF  E1=S 
E3=STRE3S  COMPONENT  NUMBER.  REQUEST  1 
E4:5TRE3S  COMPONENT  NUMBER,  REQUEST  2 
E5=STRE3S  COMPONENT  NUMBER,  REQUEST  3 
ES=STRESS  COMPONENT  NUMBER,  REQUEST  ^ 
E7=STRESS  COMPONENT  NUMBER.  REQUEST  5 
E3=5TRE35  COMPONENT  NUMBER,  EEQUE5T  5 
E9rSTRESS  COMPONENT  NUMBER.  REQUEST  7 
E10=STRE3S  COMPONENT  NUMBER.  REQUEST  B 
El l “STRESS  COMPONENT  NUMBER,  REQUEST  9 
E12=BTRESS  COMPONENT  NUMBER,  REQUEST  10 
E13r5TRE55  COMPONENT  NUMBER.  REQUEST  II 
El  (^STRESS  COMPONENT  NUMBER,  REQUEST  12 


I3P 

NS 

1 

INTEGER. PLOT  SPACING  INDICATOR  FOR 
.OUTPUT 
. — ]SP=ENTRT 

DISPLACEMENT 

.VII, 15 

ISPS 

NS 

1 

INTEGER. PLOT  SPADING  INDICATOR  FOR 
.—  JSPSzENTRY 

STRESS  OUTPUT 

ivtl.17 

1ST 

NS 

1 

INTEGER. INPUT  SPECTRUM  TTPE 
. — 1ST =ENTRT 

I VI I . 25 

IX 

S3 

6 

INTEGER. BOUNDARY  CONDITION  DATA 

. — 1XCEI ) =E2.E3.  B4.E5. E6.E7 

1 1 1 1 . 1 

WHERE 

E1=ND0RL  POINT  NUMBER  .LE.  NIMNF 
E2:X-TRPNSLATI0N  BOUNDARY  CONDITION  CODE 
E3=Y-TRPN3LAT ION  BOUNDARY  CONDITION  CODE 
EH:Z-TRflNSLATI0N  BOUNDARY  CONDITION  CODE 
ES=X-RB7ATICN  BOUNDART  CONDITION  CODE 
E6=T-RDTRTICN  BOUNDART  CONDITION  CODE 
E7:Z-ROTATI0N  BOUNDART  CONDITION  CODE 
NUMBER  OF  DEGREES  OF  FREEDOM  (EQUATION) 

PER  SLOCK  DF  'STORAGE 

— KEQB=ENTRY 

OUTPUT  TYPE  INDICATOR  (DISPLACEMENT) 

EQ. 1 . PRINT  HISTORIES  AND  MAXIMA 
EQ.2,  PRINTEFi  FLOT  HISTORIES  AND  RECOVERT  OF 
MAXIMA 

EQ.3,  RECOVERY  0F  MAXIMA  ONLY 

— KKKrENTRT 


IV. 6. 2 


VII. 17 


VII. 15 


168 


METWDRD 

RNO 

NO. 

TTFE 

. U5RGE 

. REF. 

TYPE 

OF 

OF 

CF 

ALLOCATION 

ENTRY 

ENTRY 

KKKS 

NS 

1 

INTEGER 

! OUTPUT  TYPE  INDICATOR  (STRESS) 

. EQ . 1 , PRINT  HISTORIES  AND  MAXIMA 

IVII. 17 

.ED. 2.  PRINTER  PLOT  HI  STORIES  AND  RECOVERY 

DF. 

MAXIMA 

.EQ.3.  RECOVERY  CF  MAXIMA  ONLY 

LL 

NS 

1 

INTEGER 

.NUMBER  OF  STRUCTURAL  LORD  CASES 
. — LL=ENTFiT 

.11.1 

MASSON 

SS 

1 

BEAL 

.HASS  DENSITY  CF  THE  MATERIAL  USED  TO  COMPUTE. IV. 8. E 

.THE  MASS  MATRIX  IN  A DYNAMIC  RNRLT3I 3 
-HASSON  ( El  )=E2 

. WHERE 

E IMMATERIAL  IDENTIFICATION  NLMEER 

(INTEGER)  .LE.  NUMHAT 

E2MENTRT 

MHTDES 

S3 

4 

RLPHfl 

.SHELL  MATERIAL  DESCRIPTION 
.--HATDE3(E1)=»E2* 

.IV. 8. B 

. HHEFiE 

E1=MATERIPL  IDENTIFICATION  NUMBER 

(INTEGER)  .LE.  NUMHAT 

E2;3E  CHAFfiCTERS  INCLUDING  BLANKS 

MATEY 

MS 

7 

FEFL 

. 3HELL  MATERIAL  PROPERTY  DATA 
. — HATEVCE1, E2)=E3.E4, E5.E6, E7. EB.E9 

.IV. 3. 7 

. WHERE 

F1-MBTFRTO  TnFNTlFirATMN  NIMFFR 

. LE.  NUHMAT 

E2-TEMFEfTFTTLFIE  NUHBER  I INTEGER) 

.LE.  NTP (NUMHAT) 

. E3=TEMPERATIJRE 

E4:0RTHOTR0PIC  ELASTIC  MODULO  (Ell) 
£5=ORTH0TFOF'IC  El  TSTIC  MODULO  (E22) 
E6:0RTH0TR0FIC  ELASTIC  MODULO  (E33) 
E7=PBISS0WS  RATIO 

E8=PDI330NTS  RATIO 

E9=PDISS0NfS  RATIO 

MATGfl 

MS 

6 

REAL 

.SHELL  MATERIAL  PROPERTY  DATA 
. — NATGA(E1,E2)=E3.E4,E5»E6.E7.  EB 

.IV. 8. 7 

. WHERE 

E IMMATERIAL  IDENTIFICATION  NLMEER 

(INTEGER)  .LE.  NUMHAT 

E2=TEMPERRTURE  NUMBER  (INTEGER) 

.LE.  NTP  (NUMHAT) 

E3=SHEAR  MODULO  (G12) 

E4-5HERR  MODULO  IG13) 

E5:SHEAR  MODULO  (&23) 

E6=CBEFFICIENT  OF  THERMAL  EXPANSION 
E7=CDEFFICIENT  OF  THERMAL  EXPANSION 
E8=C0EFFICIENT  OF  THERMAL  EXPANSION 

NFIXNDD 

NS 

1 

INTEGER 

.SHELL  MAXIMUM  NUMBER  OF  NODE 5 USED  TO 
.DESCRIBES  ANT  ONE  ELEHENT 
. — NAXNOD=ENTRY 

.IV. 3.1 

MAXTP 

NS 

1 

INTEGER 

.SHELL  MAXIMUM  NUMBER  OF  TEMPERATURE  POINTS 
.USED  IN  THE  TABLE  FOR  ANY  MATERIAL 
. — NAXTPrENTRY 

.IV. 8.1 

MQDEX 

NS 

1 

INTEGER 

. PROGRAM  EXECUTION  MODE 

- 1 1 . 1 

. — NODEX=ENTRT 

• 
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HH 


__ 


USAGE 


WORD  FIND  ND.  TTFE 
TYPE  OF  OF  CF 
ALLOCATION  ENTRY  ENTRY 


REF. 


NAD 

NS 

1 

INTEGER. TOTAL  NUMBER  OF  VECTORS  TO  RE  USED  IN  A .11.1 

. SUBSPACE  ITERATION  SOLUTION  FDFi  EIGENVALUES/. 

. VECTORS 
. — NAD-ENTR'f 

NAT 

NS 

1 

INTEGER. NUMBER  OF  DIFFEFENT  ARRIVAL  TIMES  FOR  THE 
.FORCING  FUNCTIONS 
. — NAT=ENTRT' 

.VI I. 8 

NRTX 

NS 

1 

INTEGER. GROUND  MOTION  CONTROL  DATA  - ARRIVAL  TIME 
• NUMBER,  X-DIRECTION 
. — NATX=ENTRY 

ivii.n 

NRTY 

NS 

1 

INTEGER. GROUND  MOTION  CONTROL  DATA  - ARRIVAL  TIME 
.NUMBER,  T -DIRECTION 
. — NATY=ENTRY 

ivii.n 

NATZ 

N3 

1 

INTEGER. GROUND  MOTION  CONTROL  DRTR  - RRRIVFL  TIME 
.NUMBER,  Z-DIRECTIGN 
. — NRTZ=ENTRY 

'.Vll.ll 

NEEAN 

NS 

1 

INTEGER. NUMBER  OF  BEAM  ELEMENT 
. — NBEAH=ENTRY 

! IV. 2.1 

NDLS 

NS 

1 

INTEGER. NUMBER  OF  DIFFEFENT  DISTRIBUTED  LOAD 
. — NDLSrENTRY 

ilV.8.1 

FOTN 

NS 

1 

INTEGER. ANALYSIS  TTFE  CODE 
. — NDTN=ENTRY 

ill.  1 

FELTYP 

NS 

1 

INTEGER. NUMBER  OF  ELEMENT  GROUPS 
.—NELTTF -ENTRY 

ill.  1 

NF 

NS 

1 

INTEGER. NUMBER  OF  FREQUENCIES  TO  BE  FOUND  IN  THE 
.EIGENVALUE  SOLUTION 
.—  NF=ENTRY 

ilt.i 

NFN 

N3 

1 

INTEGER. NUMBER  OF  DIFFEFENT  TIME  FUNCTIONS 
. — NBfcENTRr 

ivil.7 

NFNX 

NS 

1 

INTEGER. GROUND  MOTION  CONTROL  DATA  - TIME  FUNCTION 
.NUNBER  DESCRIBING  THE  GROUND  RCCELERHTION 
.IN  THE  X-DI  FEET  ION 
. — NFNX=ENTRY 

ivn.ii 

• 

» 

■ 

FFNY 

NS 

1 

INTEGER. GROUND  MOTION  CONTROL  DATA  - TIME  FUNCTION 
.NUMBER  DESCRIBING  THE  GROUND  ACCELERATION 
. IN  THE  Y-DIRECTION 
. — NFN Yr ENTRY 

.VII. 11 

nfnz 

N3 

1 

INTEGER. GROUND  MOTION  CONTROL  DATA  - TIME  FUNCTION 
.NUMBER  DESCRIBING  THE  GROUND  ACCELERATION 
.IN  THE  H-DIFECTION 

i VI I. 11 

NFS 

NS 

1 

INTEGER. NUMBER  OF  STORTING  ITERATION  VECTORS  TO  BE 
.READ  FROM  TSFE  10 
.—NFC  ENTRY 

iVfI.3 

N3N 

NS 

1 

INTEGER. GROUND  MOTION  INDICATOR 
. — NG«=ENrrir 

i VI 1.8 

NITEH 

NS 

1 

INTEGEH. MAXIMUM  NUMBER  OF  ITERATIONS  ALLOWED  TO 
.REACH  THE  CONVERGENCE  TOLERANCE 
. — MITEMrLJTFIY 

ivtl.3 

fLP 

3S 

2 

INTEGER. NUNBER  OF  FUNCTION  DEFINITION  POINTS 
.—  NLP(E1)=E2.E3 

i VI I. 13 

WHERE 

Ei=TIME  FUNCTION  NUNBEfi 
. LE.  NFN 

E2;NUMBER  OF  FUNCTION  DEFINITION  POINTS  . 
EBrSCALE  FfiCTGR  TO  BE  RFFLIED  TO  F(T) 
VALUES 


■ 
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I 


' • 

p 

HETWCIRO  AND 
TYPE  OF 
ALLOCATION 

NO. 

OF 

ENTRY 

TYPE  . 

CF 

ENTRY  . 

USAGE 

. REF. 

NOPSET 

NS 

1 

integer! 

SHELL  NUMBER  OF  SETS  OF  DATA  REQUESTING 

I IV. 8. I 

" ; ' 

STRESS  OUTPUT  AT  VARIOUS  ELEMENT  LOCATIONS 

— NQPSETr ENTRY 

NORTHO 

NS 

1 

INTEGER. 

SHELL  NUMBER  OF  DIFFERENT  SETS  OF  MATERIAL 

.IV. 0.1 

! 

AXIS  ORIENTATION  DATA 

i 1 

— N0RTH0=ENTRY 

1 j 

NOT 

NS 

1 

INTEGER. 

OUTPUT  PRINT  INTERVAL  FOR  STRESSES. 

.VII. 8 

I 1 

DISPLACEMENTS 

— NOTzENTRT 

NF 

MS 

3 

MIXED  . 

T I HE- VARYING  LOAD  DATA 

.Vtl.10 

— NP(EL.E2)=E3,EW.E5 

WHERE 

ElrNODAL  POINT  WHERE  THE  LOAD 

COMPONENT  (FORCE  OR  MOMENT!  IS 

[ 1 

APPLIED  (INTEGER)  , LE.  NUPNP 

1 1 

E2=DEGREE  OF  FREEDOM  (INTEGER) 

1 i 

.EE.  1 .PND.  .LE.  6 

I 1 

E3=TIME  FUNCTION  NUMBER  (INTEGER) 

1 1 

.LE.  NFN 

1 1 

EU;ARRIVAL  TINE  NUMBER  (INTEGER) 

ES:3CRLRFi  MULTIPLIER  FOR  THE 

1 1 

TIME  FUNCTION  IREAL) 

■ » 

NPLRTE 

NS 

1 

INTEGER. 

NUMBER  OF  PLATE  ELEMENT 

.IV. 6.1 

— NPLATEiENTRY 

: 

NF-TS 

NS 

1 

INTEGER. 

NUMBER  OF  DEFINITION  POINTS  IN  TPE  SPECTRUH 

. V ( I . 2U 

TRSLE 

} 

— MPTS=ENTRY 

N3DL21 

NS 

1 

INTEGER. 

NUMBER  OE  SPELL  ELENENTS 

.iv.e.i 

— NS0L2I=ENTRY 

NT 

NS 

1 

INTEGER. 

TOTOL  NUMBER  OF  SOLUTION  TIME  SREPS 

. V [ 1 . 8 

— NT-ENTRT 

1 w 

WTP 

ss 

1 

INTEGER. 

NUMBER  OF  DIFFEFENT  TEMPERATURES  AT 

.IV. 8. 6 

i t 

MHJCH  PROPERTIES  ARE  GIVEN 

— NTP(E1)=E2 

1 

WHERE 

\ r 

ElrMATERIAL  IDENTIFICATION  NUMBER 

(INTEGER)  .LE.  NUMMAT 

E2=NUM6ER  OF  DIFFERENT  TEMFERRTUFE3 

• LE.  MAKTP 

NUMMHT 

NS 

1 

INTEGER. 

NUMBER  OF  DIFFEFENT  MATERIALS  (SHELL) 

.IV. 6.1 

— NUMMAT GENTRY 

f f 

NUHNP 

NS 

1 

INTEGER. 

NUNBER  OF  NOD PL  POINTS 

.K.l 

— NUMNF-ENTRT 

FELML 

NS 

4 

FEFL  . 

PLATE  ELEMENT  LOAD  MULTIPLIER 

.IV. 6.1 

— PELML=E1.E2,E3,E4 

WHERE 

i 

E1=DI STRIBUTEO  LATERAL  LOAD  MULTIPLIER 

FOR  LOFC  CRSE  R 

E2rDISTRIBJTED  LATERAL  LOAD  MULTIPLIER 

\ £ 

FOR  LORD  CASE  B 

V 

E3-DISTRIEUTED  LATERAL  LOAD  MULTIPLIER 

1 § 

FOR  U3PO  CASE  C 

p 

E4:DISTRIBUTED  LATERAL  LOAD  MULTIPLIER 

V £• 

FOR  LOAD  CASE  D 

METWORO  RND 
TYPE  OF 
ALLOCATION 


PELM7  NS 


PELMX  NS 


PELMY  NS 


PELMI  NS 


PLATE  SS 


NO.  TTFE 
QF  EF 
ENTRY  ENTRY 


4 FEflL 


U5HGE 


PLATE  ELEMENT  LOAD  MULTIPLIER 
— PELMT=E1,  E2,E3. E4 
WHERE 


REF. 


IV. 6. 1 


4 FEflL 


E2=TEMERR7UFE  MULTIPLIER  FOR  LORD  CASE  B 
E3=TEMERATUFE  MULTIPLIER  FOR  LOAD  CASE  C 
E4=TEMERRTUFE  MULTIPLIER  FOR  LORD  CASE  D 
PLATE  ELEMENT  LOAD  MULTIPLIER 
— PELMX=E1.E2.E3,  EU 
WHERE 

El=X-DIRECTION  ACCELERATION  FOR  LORO 
CR3E  R 


IV. 6.1 


E2:X-DIRECTI0N  ACCELERATION  FOR  LOAD 

CASE  B 

! E3:X-DIRECTI0N  RCCELERRTION  FOR  LORD 

case  c 

! E4=X-0IFECT ION  RCCELERRTION  ^OR  LORD 

Case  o 

4 REAL  ! PLATE  ELEMENT  LOAD  MULTIPLIER 

PELMT=E1,E2.E3,E4 

WHERE 

! El:T-DIRECTION  RCCELERRTION  FOR  LORD 

CASE  A 

! E2=Y-DIRECTION  RCCELERRTION  FOR  LORD 

CASE  B 

! E3=Y-DIFECTI0N  ACCELERATION  FOR  LOAD 

CR5E  C 

E4:Y-DIFECTI0N  ACCELERATION  FOR  LOAD 

CAGE  D 

4 REAL  .PLATE  ELEMENT  LCRQ  MULTIPLIER 


IV. 6.1 


IV. 6.1 


— PELMZ=E1,E2,E3,E4 


4 REAL 


WHERE 

E1=Z- DIRECT ION  ACCELERATION  FOR  LOAD 

CR3E  R 

E2rZ-DIRECTI0N  ACCELERATION  FOR  LOAD 

CRSE  B 

E3:Z-DIFECTI0N  ACCELERATION  FOR  U3AD 

CASE  C 

E4=Z-DIF£CTI0N  ACCELERATION  FOR  U3RO 
CASE  D 
PLATE  DATA 

— PLATE  (El)  =E2.  E3 . E4 . E5 

WHERE 

EHFLHTE  ELEMENT  NUMBER  I INTEGER) 

. LE.  ff  LATE 

E2=ELEMENT  THICKNESS 

E3=DISTRIBUTED  LATERAL  LOAD 

EUrMEAN  TEMPERATURE  VARIATION  T FROM  THE 
REFEFENCE  LEVEL  ) UNDEFOFMED  FOSITIO 

E5=MEAN  TEMPERATURE  GRADIENT  ACROSS  THE 

SHELL  THICKNESS 


IV. 6. 2 
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NETWDRD  RND 

NO. 

TTFE 

. U3RGE 

. REP. 

TYPE  OP 

OF 

CF 

ALLOCATION 

ENTRY  ENTITY 

PMPI  SS 

ID 

REAL 

! PLATE  MATERIAL  PROPERTY  INFORMATION 

iiv.e.i 

— PMPI (El ) =E2.E3. E4.E5.E6iE?, E8, E9.E10.E11 
WHERE 

El'iMRTERIFL  IDENTIFICATION  NUMEER 


(INTEGER)  .LE.  PNDH 

E2=MRS5  DENSITY 

E3:THERMRL  EXPANSION  COEFFICIENT  (X) 
EU=THEfirtAL  EXPANSION  COEFFICIENT  (Y) 
E5-THERMAL  EXPANSION  COEFFICIENT  (2) 
E6=ELASTICITT  ELENNT  CXX 
E7=ELA3TICITY  ELENNT  CXT 
ESrEL ASTI CITY  ELENNT  CXS 
E9=ELASTICITY  ELENNT  CYT 


PNOM 

N3 

1 

INTEGER. NUNBER  OF  PLATE  MATERIAL 
PNDMrENTRY 

PROPERTY  SET 

.IV. 6.1 

R0TL 

NS 

i 

REAL 

•CONVERGENCE  TOLERANCE  FOR  THE  HIGHEST 
.REQUESTED  EIGENVALUE 
. — RQTLrENTRY 

Ivi I . 3 

S03LDF 

33 

7 

FEAL 

.DISTRIBUTED  3LRFACE  LOAD 
.IF  LT  . EQ.  1 

DATA 

ilV.8.9 

E10-ELR5TICITY  ELEMENT  CT5 
El  INELASTICITY  ELEMENT  CXT 


— SDSLDF (El  I =E2  < E3  > E4 » E5 

WHERE 

E1=L0A0  SET  ID  NUMBER  I INTEGER) 

.LE.  NCL5 

E2:PRESSUFE  AT  FACE  NODE  N1 
E3=PfiE53UFE  AT  FACE  NODE  N2 
E4-PRESSURE  AT  FACE  NODE  N3 
E5=PRE£SUFE  PTR  FACE  NODE  N4 
IF  LT  .EQ.  2 

— SDSLDF ( El ) -E2 . E3 , E4 . E5 . E6 , E7 , E8 

HHERE 

EULDAD  SET  ID  NUMBER  (INTEGER) 

.LE.  NCLS 

E2:  WEIGHT  DENSITY  OF  THE  FLUID 
ESrX-OROINATE  OF  POINT  S IN  THE  FREE 
SURFACE  CF  THE  FLUID 
E4=Y -ORDINATE  OF  POINT  S IN  THE  FREE 
SURFACE  OF  THE  FLUID 
E5:Z-0ROINATE  OF  POINT  S IN  THE  FREE 
SURFACE  OF  THE  FLUID 
E5:X-0RDIMTTE  OF  fl  POINT  N DN  THE 
NDRHFL  TO  THE  FLUID  SURFRCE 
E7=Y-0RDINATE  OF  fl  POINT  N ON  THE 
NDRHFL  TO  THE  FLUID  SURFRCE 
E3=Z-ORDI NATE  OF  fl  POINT  N DN  THE 
NDRHFL  TO  THE  FLUID  SURFRCE 
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KEYWORD  AND 
TYPE  OF 
ALLOCATION 


SOSLDI  SS 


5ELCNF  N5 


SELCHT  NS 


SELCHX  NS 


SELCHT  NS 


ND. 

TYPE 

. U5RGE 

. REr. 

OF 

CF 

. 

ENTRY 

ENTRY 

• 

• 

2 


INTEGEA 


DISTRIBUTED  SLFlFRCE  DATA 
— SDSLDI (E1)nE2.E3 
WHERE 

EInLDAD  SET  ID  NUNBER  t INTEGEFI) 

.LE.  NOLS 

E2=ELEMENT  FACE  NUMBER  ON  WHICH  THIS 
DISTRIBUTED  LOAD  IS  ACTING 
E3=LOHD  TYPE  CODE 


IV. 8. B 


4 PER- 


IL FEPL 


.EQ.  1.  PRESCRIBED  NORMAL  PFE33URE 

INTEfvIS  IT  IE3 

.EQ.  2.  HTDROSTAT I CALLY  VARYING 
PRESSURE  FIELD 
.EQ.  0.  DEFAULT  SET  TO  1 
SHELL  ELEMENT  LDFD  CASE  MULTIPLIERS 
— SELMCPnE1,E2.E3.E4 
WHERE 

E INFRACTION  OF  PRESSURE  LOADS  TO  BE 
APPLIED  IN  ELEHENT  LOAD  CASE  A 
E2nFRACTI0N  OF  PRESSURE  LOADS  TO  BE 
APPLIED  IN  ELEHENT  LOAD  CASE  B 
ESnFRRCTICN  OF  PRESSURE  LOADS  TO  BE 
APPLIED  IN  ELEHENT  LOAD  CASE  C 
E4=FR ACTION  OF  PAE33URE  LOADS  TO  BE 
APPLIED  IN  ELEHENT  LOAD  CASE  D 
SHELL  ELEMENT  LOAD  CASE  MULTIPLIERS 


IV. 5. 15 


IV. 8. IS 


4 REFL 


—5ELCMT:E1.E2.E3.EIJ 

WHERE 

E INFRACTION  OF  THERHAL  LORDS  TO  BE 
APPLIED  IN  ELEHENT  LOAD  CASE  A 
E2=FRACTIOJ  OF  THERHAL  LORDS  TO  EE 
RFFLIED  IN  ELEHENT  LORD  CR5E  B 
E3-FRACTI  ON  OF  THERMAL  LORDS  TO  BE 
APPLIED  IN  ELEHENT  LOAD  CASE  C 
E4nFR ACTION  OF  THERHAL  LORDS  TO  EE 
APPLIED  IN  ELEHENT  LOAD  CASE  D 
SHELL  ELEMENT  LOAD  CASE  MULTIPLIER 
— SELCHX;E1.E2.E3.EU 


IV. 8. 13 


U REAL 


HHERE 

EInFRACTICN  OF  X-DIRECTION  GRAVITY 
APPLIED  IN  ELEHENT  LOAD  CP6E  A 
E2nFRHCTICN  OF  X-DIRECTION  GfHVITY 
APPLIED  IN  ELEHENT  LOAD  CASE  B 
E3NFRHCTIOJ  OF  X-DIRECTION  GRAVITY 
APPLIED  IN  ELEHENT  LOAD  CASE  C 
E4=FRACTIOI  OF  X-DIRECTION  GRAVITY 
APPLIED  IN  ELEHENT  LOAD  CASE 
SHELL  ELEMENT  LORD  CASE  MULTIPLIER 


TO  BE 
TO  BE 
TD  BE 
TQ  BE 


IV. 8. 13 


—3ELCMTNE1.E2.E3.E4 

HHERE 

E1=FRPCTICN  OF  Y-DIRECTION  GRAVITY  TO  BE 
APPLIED  IN  ELEHENT  LOAD  CASE  A 
E2=FRACTI0N  OF  Y-DIRECTION  GRAVITY  TO  BE 
APPLIED  IN  ELEHENT  LOAD  CASE  B 
E3=FRACTIDN  OF  Y-DIRECTION  GRAVITY  TQ  BE 
APPLIED  IN  ELEHENT  LOAD  CASE  C 
E4NFRACTICN  OF  Y— DIRECTION  GRAVITY'  TQ  BE 
APPLIED  IN  ELEHENT  LOAD  CASE  D 
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HETWDRD  RND 
TYPE  OF 
ALLOCATION 

NO. 

OF 

ENTRY 

TYPE 

CF 

ENTRY 

. USAGE 

. REF. 

SELCHZ 

NS 

4 

REAL 

[SHELL  ELEMENT  LOAD  CASE  MULTIPLIER 

IlV.8.13 

SELCMZ=E1,E2.E3.EU 

. WHERE 

E INFRACTION  OF  2-DIRECTION  GRAVITY  TO 

EE. 

APPLIED  IN  ELEMENT  LOAD  CRSE  A 

E2=FRACTICM  OF  Z-DIRECTION  GRAVITY  TO 

BE. 

APPLIED  IN  ELEMENT  LOAD  CRSE  B 

. E3=FR ACTION  OF  Z-DIRECTION  GRAVITY  TO 

BE. 

. AFPLIED  IN  ELENENT  LOAD  CASE  C 

EUrFRACTICN  OF  2-DIRECTION  GRAVITY  TO 

BE. 

APPLIED  IN  ELEMENT  LOAD  CASE  D 

SFTR 

NS 

1 

FEAL 

.SCALE  FACTOR  USED  TO  ADJUST  THE  DISPmCEHENT.VII.21 

. IOR  ACCELERATION)  ORDINATES  IN  THE  SPECTRUM  . 

.7RELE 

SFTRzENTRY 

SFTR3 

NS 

1 

FEFL 

.SCALE  FACTOR  LB  ED  TO  ADJUST  THE 

.VI I. 21 

.DISPLACEMENT  (OR  ACCELERATION)  ORDINATES 

.IN  THE  SPECTRUM  TABLE 

SFTRSrEN 

. — SFTRSrENTRY 

SHELL  1 

S3 

5 

INTEGER. SHELL  DATA 

.IV. 8. 16 

. — SHELL 1 (El ):E2,E3. EU.E5.  E6 

SHELL2  SS  8 INTEGER 


WHERE 

EUELEMENT  NUMBER  .LE.  NS0L21 

E2=NUMBER  OF  NODES  TO  BE  USED  IN 

DESCRIBING-  TME  ELEMENT  5 DISPLACEMENT 

FIELD 

E3=NUMBER  OF  NODES  TO  BE  USED  IN  THE 
DESCRIPTION  OF  ELEMENT  GEOMETRY 
E4=MATERIFL  IDENTIFICATION  NLMBER 
. LE.  NJftIRT 

E5=IDENTIFICRTI(3N  NUMBER  OF  IRE  MATERIAL 
AXIS  ORIENTATION  SET  .LE.  NORTHO 
E6= IDENTIFICATION  NUMBER  OF  THE  STRESS 
OUTPUT  SET  .LE.  NOPSET 
SHELL  DATA 

— SHELL2 1 El  I rE2 . E3 . EU . E5 . EG . E7 . E8. E9 
WHERE 

EliELEMENT  NUHBER 

E2=N0OE  NUMBER  INCREMENT  FORM  ELEMENT 

DRTR  GENERATION 

E3= INTEGRATION  ORDER  FOR  NATURAL 
COORDINATE  IR. 5)  DIRECTIONS 
EU: INTEGRATION  ORDER  FOR  NATURAL 
COORDINATE  IT)  DIRECTION 
E5=FLAG  INDICATING  THAT  THE  STIFFNESS 
AND  MASS  MATRICES  FOR  THIS  ELEMENT 
ARE  THE  ERNE  R5  THOSE  FOR  THE 
PRECEDING  ELEMENT 

E6=PRE59JFE  SET  FOR  EUEMENT  LOAD  CASE  fl 
E7:PRES3JRE  SET  FOR  ELEMENT  LOAD  CASE  B 
E8=PRESSURE  SET  FOR  ELEMENT  LOAD  CASE  C 
E9=PRESSUFE  SET  FOR  ELEMENT  LOAD  CASE  D 


IV. 8. 16 
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PETWORO  RND 
TYPE  OF 
ALLOCATION 


SMflOS  SS 


SSQRLS  SS 


3SPEC  SS 


snoa  ss 


S17TB21  SS 


S9T016  SS 


*•  - • ^ — 


NO.  TTFE 
OF  CF 
ENTRY  ENTRY 


3 INTECEH 


7 INTEGER 


1 FEFL 


8 INTECEH 


5 INTEGER 


8 INTEGER 


U5RGE 


SHELL  MATERIAL  AXES  ORIENTATION  SET 
— SMftOS  C E 1 ) =E2.  E3 , EM 
WHERE 

El: IDENTIFICATION  NUMBER  . LE.NQRTHO 
E2:N0DE  NUMBER  FOR  POINT  I 
E3=NBDE  NIMEER  FOR  POINT  J 
E4=N0DE  NUMBER  FOR  POINT  K 
SHELL  STRESS  OUTPUT  REQUEST  L0CRTI0N  SETS 
— S30RL3 ( El ) =E2 , E3 . E4 . E5 , E6 , E7 . E3 
WHERE 

E1=3ET  NUPBER  .LE.  N0P3ET 
E2:L0CRTICM  NUMBER  OF  OUTPUT  POINT  1 
E3=LDCRTICN  NUMBER  OF  OUTPUT  POINT  2 
E4:LBCRTICN  NUMBER  OF  OUTPUT  FOINT  3 
E5:L0CRTICN  NUMBER  OF  OUTPUT  POINT  4 
ES=L8CRTI ON  NUMBER  OF  OUTPUT  POINT  5 
E7:L0CRTICN  NUMBER  OF  OUTPUT  POINT  6 
E8=LBCBTI0N  NUMBER  OF  OUTPUT  POINT  7 
VALUE  OF  DISPLACEMENT  (OR  ACCELERATION) 
FOR  RESPONSE  SPECTRUM  ANALYSIS 
— 3SPEC(E1)=E2 
WHERE 

E1=DEFINITI0N  POINT  NUMBER  .LE.  FPTS 
E2=ENTRY 
SHELL  DATA 

— 5 1 TOB I E 1 ) =E2.  E3  ,E4 . E5 . E6 . E7,  EB . ES 

WHERE 

E1=ELEHENT  NUNBER  .LE.  N3BL21 
E2=N0OE  I NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 


2 

3 

4 


6 

7 


E3=NBDE 
E4:N0DE 
E5=N0DE 
E8=ND0E 
E7:NBDE 
E8=NBDE 

E9:N0DE  8 NUMBER 

SHELL  DATA 

--3 1 7T021  ( E 1 ) =E2.  E3 . Et . E5  < E6 

WHERE 

E1=ELEHENT  NUMBER  .LE.  NSBL21 
E2=N0DE  17  NUMBER 
E3=NBDE  18  NUMBER 
E4:NDDE  19  NUNBER 
E5:N0DE  20  NUMBER 
E8=NDDE  21  NUMBER 
SHELL  DATA 

— SQTOietEl )=E2.E3. EU.ES. E6, E7. E3. ES 


REF. 


IV. 8. B 


IV. 8. 13 


VII.24 


IV. 8. 17 


IV. 8. 17 


IV. 8. 17 


WHERE 

EUELEHENT  NUMBER  .LE.  NS0L21 

E2=NBDE  9 NUMBER 

E3=N0DE  10  NUNBER 

E4=NBDE  11  NUMBER 

E5=N0OE  12  NUMBER 

ESiNBDE  13  NUMBER 

E7:NDDE  14  NUNBER 

E8:N0OE  15  NJNBER 

E9=N0DE  16  NUMBER 
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APPENDIX  B 
TIME  HISTORY  PLOT 


The  time  history  plot  package  is  now  an  integral  part  of  the  SAP  IV 
program.  Response  history  analyses  of  displacements  and  stresses  are 
produced  for  three  element  types  - beams,  thin  plates,  and  thick  shell 
(21  nodes).  Analyses  of  other  element  types  can  be  obtained  but  only  a 
generalized  labeling  is  available. 

The  plot  package  has  been  added  to  the  SAP  IV  program  with  minimal 
modifications  to  SAP  IV. 

The  input  requirement  for  time  history  plot  output  is  identical  to 
the  one  for  the  response  history  analysis  output  request  as  given  in  the 
SAP  IV  manual,  and  no  other  input  is  required.  However,  the  sequence 
of  output  for  the  time  history  plot  is  different  from  the  arrangement  of 
output  for  response  history  analysis  output  as  indicated  in  the  SAP  IV 
manual.  A maximum  of  three  curves  will  be  plotted  within  one  frame  for 
each  node  or  element  member. 
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APPENDIX  C 


MESH  PROGRAM 


PURPOSE 

Program  MESH  produces  undeflected  mesh  plots  from  data  on  TAPE  8 written 
by  the  SAP  IV  program.  Using  different  input  instructions,  the  entire 
element  model,  one  part  of  the  model,  or  a combination  of  different  parts 
of  the  model  may  be  plotted  on  one  frame.  Options  of  labeling  all  or  selected 
nodes  and/or  elements  are  available.  A unique  line  pattern  represents  each 
type  of  element— that  is,  a solid  line  for  a beam  element,  a series  of  long 
dashed  lines  for  a thin  plate  element  and  a series  of  short  dashed  lines  for 
a thick  shell  element. 

INPUT  INSTRUCTIONS 

Card  one  is  required  for  the  successful  generation  of  one  or  more  plots. 
Two  additional  cards  are  required  for  each  plot  to  be  generated. 

Card  One: 


Col.  1-10  Vertical  length  of  a plot  in  inches.  Eight  inches 
is  used,  allowing  1 inch  for  the  legend  at  the 
bottom  and  a 7-inch  square  for  the  graphics.  The 
positioning  and  sizing  of  all  character  symbols 
are  fixed. 

Col.  11-20  Scale  factor.  This  allows  the  user  to  reduce  or 

expand  frame  size  with  all  proportions  of  the  entire 
layout  kept  intact.  It  is  defaulted  to  1 if  set 
blank. 


Card  Two: 

This  represents  the  first  input  card  of  each  plot  to  be  generated 
and  is  in  free  format.  Keyword  parameters  available  in  the  list  below 
must  be  used. 

X x-coordinate  in  cartesian  system  of  viewer  position 

Y y-coordinate  in  cartesian  system  of  viewer  position 

Z z-coordinate  in  cartesian  system  of  viewer  position 

A All  element  members 
B Beam  element  members 


P Plate  element  members 

S Shell  element  members 

AN  All  element  members  with  nodes  labeled 

BN  Beam  element  members  with  nodes  labeled 

PN  Plate  element  members  with  nodes  labeled 

BN  Shell  element  members  with  nodes  labeled 

AE  All  element  members  with  element  labeled 

BE  Beam  element  member  with  element  labeled 

DE  Plate  element  members  with  element  labeled 

SE  Shell  element  members  with  element  labeled 

ANE  All  element  members  with  nodes  and  elements  labeled 

BNE  Beam  element  members  with  nodes  and  element  labeled 

SNE  Shell  element  member  with  nodes  and  element  labeled 

The  order  and  position  of  the  keyword  parameter  is  immaterial  and  is 
column  independent  on  card.  Each  keyword  parameter  must  be  followed  by 
a pair  of  parentheses,  even  though  an  entry  is  not  required.  Keyword 
parameter  X,  Y,  and  Z require  the  entry,  if  any,  to  be  in  floating 
point  format--that  is,  it  must  have  a decimal  point.  Entries  for 
other  parameters  are  to  be  in  fixed  point  format  or  integer.  Except  for 
parameters  X,  Y,  and  Z,  multiple  entries  are  allowed  and  must  be  separated 
by  a comma  with  no  blank  imbedded.  The  user  can  generate  entries  auto- 
matically between  two  given  entries  that  are  separated  by  one  blank 
character.  Samples  are  given  below. 

Card  Three: 

This  represents  the  second  input  card  of  each  plot  and  is  used  for 
repositioning  of  the  plot  pen  for  next  plot  or  at  the  end  of  the  plotting 
series. 

Col.  1-10  +X  - direction  in  inches 

Col.  11-20  +Y  - direction  in  inches 
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Graph  size  (vertical  length)  is  8 inches.  Three  plots  are  to  be 
generated  with  the  following  input  data.  The  first  plot  has  all  elements 
plotted  with  the  viewer  position  at  (1,  0,0)  in  cartesian  space.  The 
second  plot  uses  a viewer  position  of  (1,  1,  1)  with  all  elements  plotted 
and  beam  elements  No.  1,  3,  and  5 through  10  with  nodes  labeled.  The 
third  plot  has  the  same  viewer  position  as  the  second  plot.  In  the  third 
plot,  beam  elements  1 through  10  have  nodes  labeled,  beam  elements  8 through 
15  have  elements  labeled,  beam  elements  16  through  20  have  both  nodes  and 
elements  labeled,  and  shell  elements  1 through  10  have  nodes  labeled. 

Card  1 e3. ii 

Card  2 X(l.)  A(  ) 

Card  3 7lO.11  is0.2i 


Card  4 X(l.)  Y(l.)  Z(l.)  A(  ) BN(1 ,3,5*10) 

Card  5 10.  0. 

7 11  iru 

Card  6 X(l.)  Y(l.)  Z(l.)  BN(1*10)  BE(3015)  BNE(16*20)  SN(1*10) 

Card  7 7IO.11  is0.2i 

Card  8 7-8-9 


APPENDIX  0 


CONTROL  CARD  SAMPLES 


CREATE  UPDATE  FILE  FOR  FREE  FORMAT  INPUT  PROGRAM, 
SAP  IV  WITH  TIME  HISTORY  PLOT  PACKAGE,  OR  MESH 
PROGRAM,  AND  SAVE  ON  PERMANENT  FILE. 

JOB  CARD  with  account  number 
REQUEST  (NEWPL,*PF) 

UPDATE  (F,N) 

CATALOG  (NEWPL,  FFIPNEWPL,  ID=) 

(OR  CATALOG,  NEWPL,  SAP4THNEWPL,  ID=  .) 

(OR  CATALOG,  NEWPL,  MESHNEWPL,ID=  .) 

7-3-9 

*DECK  FFIP 
(OR  *DECK  SAP4TH) 

(OR  *DECK  MESH) 

-FORTRAN  SOURCE  PROGRAM 
6-7-8-9 


COMPILE,  CREATE  BINARY  OBJECT,  AND  SAVE  ON 
PERMANENT  FILE. 

JOB  CARD  with  account  number 
ATTACH  (OLDPL,  FFIPNEWPL, I D= 

(OR  ATTACH,  OLDPL,  SAP4THNEWPL,ID=  .) 

(OR  ATTACH,  OLDPL,  MESHNEWPL, ID=  .) 

UPDATE (F) 

REQUEST  (LGO,*PF) 

FTN,A,I=COMPILE. 

CATALOG, LGO.FFIPLGO, ID  = 

(OR  CATALOG,  LGO,  SAP4THLG0,  ID=  .)  j 

7-8-9 

UPDATE  DIRECTIVE 
6-7-8-9 

1 

| 

EXECUTE  FREE  FORMAT  INPUT  PROGRAM,  SAP  IV  WITH  DATA 
MODEX=l , AND  MESH  PROGRAM. 

JOB  CARD  WITH  ACCOUNT  NUMBER 
ATTACH  FFIP ,FFIPLGO,ID= 

LDSET,PRESET=ZERO 
FFIP. 

RETURN, FFIP. 

REWIND  (TAPE!) 
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ATTACH,  PLOT,  PLOTLIBRARY,  ID=CA. 
ATTACH,SAP4,  SAP4THLG0,ID  . 

LDSET,  PRESET=ZERO,  LIB=PLOT 
SAP4 (TAPE! ) 

RETURN, TAPE! ,SAP4 
REWIND(TAPE8) 

ATTACH, MESH, MESHLGO,ID= 

LDSET,  PRESET=ZERO,LIB=PLOT 
MESH. 

RETURN, MESH, TAPE8. 

REWIND, TAPE! . 

ATTACH, CON VLGO, I D=CA. 

CONVLGO. 

ROUTE ,TAPE2 ,DC=PR,FC=PL ,EL=A6 ,IC=ASCII 
7-8-9 

FFIP  INPUT  DATA 
7-8-9 

MESH  INPUT  DATA 
6- 7-8- 9 


m 

